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MEETING OF 1957 JANUARY 11 
Professor Sir Harold Jeffreys, President, in the Chair 


‘The President announced that the Council had awarded the Gold Medal 
of the Society to Professor Albrecht Unsold for his work on the physics of 
stellar atmospheres. 

The President announced that the Council had awarded the Jackson-Gwiit 
Medal and Gift to Mr R. P. de Kock for his systematic observations of variable 
stars over a period of more than twenty years. 


tv trier Wace bas 


The election by the Council of che following Fellows was duly confirmed :— 

Colin Robert Edwards, c/o Engineering Division, Public Works Department, 
Perth, W. Australia (proposed by R. W. Boggis); 

*Kenneth Hugh Fea, 4 Manor Drive, Hinchley Wood, Esher, Surrey (proposed 
by C. W. Allen); 

Roger Francis Griffin, 12 Fir ‘Tree Road, Banstead, Surrey (proposed by 
F. M. Holborn); and 

*Philip James Message, Gonville and Caius College, Cambridge (proposed 
by H. Jeffreys). 


The election by the Council of the following Junior Members was duly 
confirmed :— 
David Francis Carter, 2 Peterborough Gardens, Ilford, Essex (proposed by 
G. W. Wakeford); and 
Peter Calvin Kendall, 28 Grove Hill, South Woodford, London, E.18 
(proposed by V. C. A. Ferraro). 


Seventy-three presents were announced as having been received since the 
Jast meeting, including :— 
Lund Observatory: (1) Tables for the conversion of galactic into equatorial 
coordinates and for the direction cosines in the equatorial system; and 
(2) Table for the reduction of velocities to the local standard of rest 
(presented by Lund Observatory). 


* Transferred from Junior Membership. 
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INTERFEROMETRIC STUDY OF PROFILES OF 
FAINT FRAUNHOFER LINES 


Z. Suemoto 
(Communicated by the Director of the € ambridge Observatories) 
(Received 1957 February 4) 


Summary 

Profiles of well isolated Fraunhofer lines in the red were measured 
interferometrically by crossing a Fabry-Perot interferometer (with dielectric 
coating) with the 9 m spectrograph of the Cambridge Observatories. The 
average effective resolving power was as high as 6 x 10°. The wave-length 
distribution of the selective absorption coefficient was derived from the 
measured profiles, using Vitense’s Model II. For every point of the disk 
faint lines give broader Doppler widths than do stronger lines. It is suggested 
that there is a circulation in the solar atmosphere, vertical upward motion in 
granules being converted into horizontal motion at about 70-4, and finally 
changing into downward motions back to the convection layer. 





1. Introduction.—From the measurement of profiles of faint Fraunhofer lines 
in the blue and violet regions Allen (1) found that the true widths of the faintest 
lines, as deduced from extrapolation of the measured widths, increase as we go 
from centre to limb on the Sun’s disk. He interpreted this as due to an anisotropy 
of turbulence in the reversing layer, horizontal velocities being greater than vertical 
velocities. However, as is well known, even the faintest lines which can be 
measured with reasonable accuracy are not free from saturation, in the sense that 
the intensity absorbed at a given point of the line profile does not increase in pro- 
portion to the absorption coefficient. In other words, even the faintest lines are 
wider than true absorption coefficient widths by an appreciable amount which 
varies according to the position on the disk as well as from line to line. If we try 
to extrapolate from the measured widths to the widths of lines of infinitesimal 
intensity the results are therefore very uncertain. 

In the present work the profiles of well-isolated faint lines in the red were 
measured at several positions on the disk, using a Fabry-Perot interferometer 
before the slit of the gm spectrograph. ‘The profiles of the absorption coefficient 
were deduced from the measured profiles, using one of the modern models of the 
solar atmosphere. ‘The results show that at any one point on the disk the Doppler 
widths decrease with increasing line intensity, suggesting that there is a decrease 
of turbulence upwards and that at about r=0-4 the horizontal velocities con- 
siderably exceed the vertical. From the velocity fields at various depths it seems 
that the vertical upward motions at about r=1 are deflected to the horizontal 
near t=0-4, and are finally converted into downward motions back to the con- 
vection layer. 

2. Observations.—Most faint lines, say of Rowland intensity less than 1, are 
blended appreciably, to such an extent as to make accurate measurement of their 
profiles difficult. Only in the red could numbers of well isolated lines be found. 
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Measurements were made in the solar tunnel of the Cambridge Observatories, 
by crossing a Fabry-Perot interferometer with the gm spectrograph. The 
optical arrangement is shown in Fig. 1. Since the interferometer was used 
before the spectrograph, we obtained parabolic channels similar to those of 
Shane (2), instead of circular channels. But the slit width was wider than 
Shane’s and only a little narrower than inter-channel distances (Plate 1(b)). As 
was pointed out by Treanor (3), this arrangement has the disadvantage that the 
high resolving power of the interferometer can be lost by the finite resolving 
power of the spectrograph in the direction along the slit. It will be shown later, 
however, that in spite of some loss of resolution we still reached an average 
effective resolving power of 6 x 10°. 





Fic. 1.—Optical arrangement. Full lines: optical path for solar parabolic channels. Broken 
lines: path for Pointolite lamp used with 0-08 mm spacer. 

The solar image, of 16cm diameter, was focused about 2 m before the spectro- 
graph slit S, on another auxiliary slit S,. The sunlight was then collimated 
and re-focused on to the spectrograph slit by two lenses L, and L,, between 
which the Fabry—Perot interferometer I was situated. ‘These two lenses, both 
of 50cm focal length, were situated in such a way as to act as field lenses as well. 
‘The image of the telescope lens O is first focused on to the Fabry—Perot inter- 
ferometer by L, and then the second lens L, projects this image on to the spectro- 
graph collimator lens C. As the distance of the second lens L, from the slit S, 
is small compared with the spectrograph focal length, it acts approximately as 
a field lens, such as would normally be put just in front of the slit. This modified 
method worked quite satisfactorily and we were quite free of difficulties from dust 
or bubbles on the field lens. ‘The lens L, also serves as a “ ring lens” to project 
the interference fringes from infinity on to the spectrograph slit. The spectro- 
graph is of the auto-collimation type, with an 18cm collimator lens of gm focal 
length. A high quality grating, lent to the Cambridge Observatories by the 
Mount Wilson Observatory for eclipse work, was used in the second order. 

Dr A. H. Jarrett kindly coated the interferometer plates especially for this 
investigation. They had a high efficiency 7-layer dielectric coating, giving a 
reflection coefficient of 94 per cent in the red, with very small absorption. A 
7°5mm quartz spacer was used and ideally, with perfect plate figures, the 
resolving power would have been about 1-2 x 108. The light from ghost images 
of the first and second kind (4), was prevented from entering the slit by rotating 
each of the interferometer plates. 
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Because of the high efficiency of both grating and interferometer the exposure 
times were only about 4o sec at the centre of the disk. About thirty standardized 
photographs were obtained during the summer of 1953, covering five positions 
on the disk, in two spectral regions. For standardization, each plate (Ilford 
Rapid Process Panchromatic) was cut into halves, one of which was used for 
photographing the parabolic channels, the other for the standardization spectrum. 
The latter was obtained by replacing the interferometer by a neutral filter and 
putting a Hilger step wedge of rhodium over the slit, the exposure time varving 
by not more than a factor 2 between interferometer and wedge photographs. 
The two halves of the plate were brush developed together in Azol. 

3. Measurement.—Most of the microphotometer tracings were made in the 
‘Tokyo Astronomical Observatory by the self-registering microphotometer 
described elsewhere (5). Each tracing was made along a parabolic channel with 
the microphotometer slit held parallel to the spectrograph dispersion. The 
tracing, therefore, is very similar to that of a normal spectrum, except that one 
has to steer on a slightly curved path, although holding the slit orientation fixed. 
‘The microphotometer magnification was 20 x , and the projected dimensions of 
its slit o-o1 x 0-05 mm, the slit length being considerably less than the average 
width, 0-19 mm, of one channel. ‘The slit was steered along this curved spectrum 
during the tracing, watching the plate through the microscope and adjusting the 
Y-axis slow motion by hand. ‘Tracings were made in this way using three or 
four successive orders for each line on each plate. Since for each position on 
the disk there were three plates taken on different days, there were more than 
nine tracings for each line. The nine profiles were fitted together to give a final 
average profile. 

The step wedge was calibrated at a sufficient number of wave-lengths by the 
National Physical Laboratory. In constructing characteristic curves, small 
corrections for light scattered in the spectrograph were applied to the wedge 
transmission. ‘The measured intensity of this scattered light for 18mm slit 
length was 2-1 per cent. ‘The effective slit length for the calibration spectra was 
about 34mm, after making allowance for the opacity of the steps, so that when 
viewed through this spectrograph the apparent transmission of each step exceeded 
the true transmission by 0-4 per cent. Although very small, this correction 
becomes comparatively important for the steps of small transmission: e.g. 
for the last step the real transmission was 6-4 per cent, the corrected trans- 
mission 6:8 per cent. Following this procedure, a more satisfactory fit was 
obtained between sections of the characteristic curve having different average 
densities. 

4. Corrections.—There are three kinds of correction to be considered in 
spectra of this sort, arising respectively from uniform scattering, horizontal 
scattering from neighbouring parabolic channels, and the finite resolving power 
along the parabola. 

The first of these was estimated as follows: The measured scattering in the 
red for 18 mm slit length was about 2-1 per cent, as already mentioned. ‘The 
parabolic channels were photographed with the same slit length. But there 
are dark bands between the channels, and since the average channel width was 
o"19mm, with inter-channel distance 0-30 mm, the 2-1 per cent is reduced by a 


factor o-19/0-30. The correction for uniform scattering is therefore on the 
average 1-3 per cent. 
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The horizontal scattering from neighbouring channels arises from the grating 
ghosts and from the outer wings of the instrumental profile. The ghost in- 
tensities were estimated by using the yellow mercury lines; the sum of the first 
four ghosts in the segond order was o°5 per cent. In the ghost spectrum there 
was always another image of the same colour as the primary line, in appearance 
like a ghost line but obliquely located. This image was attributed to double 
reflection at the two parallel interfaces of the components of the doublet spectro- 
graph lens. The light going through the spectrograph passes twice through this 
layer so that the image should be the superposition of two, each having an 
intensity about 5-10°? x §-10°*= 2-5 x 10°* of the primary line, in good agreement 
with observation. ‘This false image is a little out of focus because of the lens 
action of the air between the two inner faces, and is slightly displaced because of 
the wedge shape of the air lens. It must be emphasized that the false image has 
nothing to do with the white light reflections at the outer surfaces of the lens. 
A spurious line of this sort can conveniently be included among the ghosts. Th« 
resulting ghost correction was finally estimated to be about 1 per cent. 

The instrumental profile of the spectrograph was measured with specia! 
attention to the wings, by using A 4047 and A5461 of a water-cooled Hg 198 isotope 
lamp for the core, and the red lines of an ordinary Ne discharge tube for the 
wings. All these measurements were consistent within the accuracy needed, 
so that only the Ne profile was used in the final estimate for the effect of the 
neighbouring parabolae. 

‘The instrumental profile thus found was integrated numerically over the 
width of the channel, 0-19 mm, giving a kind of channel profile. ‘Three of these 
were plotted side by side at the channel spacing, 0-30 mm. At the centre of the 
central channel the neighbouring channel on the violet side contributes an 
intensity 1-8 per cent, while that on the red side gives o-8 per cent. This 
asymmetry is in the instrumental profile, and appears in all the Hg and Ne lines 
measured. Adding the two contributions we have 2-6 per cent false light at the 
channel centre. However, the amount varies a little across the channel and the 
average correction for the scattering from neighbouring channels, over the 0-05 mm 
width used for the microphotometry tracing, was estimated to be 2-7 per cent. 

It is of interest to note that the minimum intensity between the channels can 
be understood in the same way. ‘The superposed contributions from two 
channels, one from each side, give 25 per cent for the o-19mm channel width, 
whereas the measured minimum intensity was 27 per cent. Similarly, with the 
parabolic channels obtained with slit width o-o5 mm, the minimum intensity 
between channels was about 6 per cent, whereas the calculated intensity was 
5 per cent. 

The vertical instrumental profile can be divided into three parts: the vertical 
profile of the spectrograph, the resolving power of the Fabry—Perot interferometer, 
and the resolving power of the fringe imaging lens, including any de-focusing of 
the fringes on the slit. 

The vertical profile of the spectrograph was measured by using a horizontal 
slit in place of the usual vertical slit. ‘This was supplemented by another 
measurement, using a thin wire across the vertical slit. The measured vertical 
protile had a width 0-065 mm on the plate, agreeing with the theoretical value, 
0-069 mm, for A6300, with a square diaphragm 82 mm high over the g m spectro- 
graph lens. Converted into wavelengths this width corresponds to 3 mA, 7mA 
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and 11 mA respectively at the first, second and third of the interferometer fringes. 
We shall be dealing with solar line widths about 70 mA, and in comparison these 
widths should be negligible. However, as described later, they do affect the 
narrower iodine lines by which de-focusing of the fringe system on to the spectro- 
graph slit was checked. 

Various considerations show that the resolving power of the Fabry—Perot 
interferometer was near 10°, but no reliable method could be devised to measure 
the actual resolving power given with the 7-5mm spacer. It is well known, 
however, that the resolving power of a Fabry-Perot interferometer depends on 
the reflectivity of the coating, the surface quality of the two plates, and the 
adjustment. In this work the first two were measured together by using a 0°08 mm 
spacer in place of 7-5 mm, and the last factor was checked by an iodine absorption 
tube. Owing to the possibility of variations of reflectivity or plate quality, the 
same area of the interferometer had to be examined as was used for the solar 
spectrum. ‘Thus the same set-up was used, changing only the spacer, and using 
a Pointolite lamp for source instead of the Sun, in order to avoid absorption 
lines. The lamp (Fig. 1), put near the focus of the condenser lens O, illumi- 
nated the secondary mirror M in just the same way as sunlight from the coelostat 
mirror. After reflection from the secondary mirror, the light was focused on 
the auxiliary slit S, and thence followed the same path as light from the Sun, so 
that exactly the same area of the interferometer was used. What one now gets 
on a spectrogram is parabolic channels, essentially similar to those from the 
thick spacer, but now with a very large inter-channel separation, 28 mm (Plate 1(¢)). 
If one slowly narrows down the spectrogram slit the width of the channels at first 
decreases too, but finally one reaches a limit beyond which the channel width 1s 
more or less independent of slit width. ‘The channel profile with the slit width 
below this limit is the profile of the interferometer, as seen through the spectro- 
graph. ‘The ratio of the inter-channel distance to the half-width of this profile 
equals the “effective number of reflections’’, which in turn is directly related to 
the reflectivity of the coating. : 

The measured interferometer profile was fitted to the formula 


i 

J /Ja= 1/ (34 = ee 4 

Ut Rp a 
where ¢ is the phase angle, using two values for R, the reflectivity, g4 per cent 
and g2 per cent, the corresponding effective number of reflections being 49 and 36 
respectively (6). In general the first value gives the better fit, but if one is 
concerned only with the core of the profile the lower reflectivity is better. This 
discrepancy is presumably associated with the surface quality, which does not 
appear in the formula. Were the coated surfaces not plane the fringe system 
formed by one part of the interferometer aperture would not be the same as that 
formed by other parts. The fringe system from the whole aperture would there- 
fore be an aggregate of profiles. whose centres arc not exactly coincident. The 
lower reflectivity, 92 per cent, must therefore be interpreted as a kind of effective 
reflectivity, taking into account the surface quality, while the higher value, 94 per 
cent, more nearly represents the reflectivity of the coating itself. We are 

interested in the resolving power and should therefore use the g2 per cent. 

With the 7-5 mm spacer a reflectivity of g2 per cent means a resolving power 
88x 10°. In wavelength units this is 7mA, which is so small as not to affect 
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the profiles of the solar lines. _ It is to be noted in this connection that the intensity 
between orders was so low that it could be neglected altogether, especially since 
we are interested only in faint lines whose residual intensities are higher than 
50 per cent. 

The focusing of the fringe system on to the spectrograph slit was checked 
before or after the solar observations, using the same set-up to photograph 
through an iodine absorption tube (Plate 1 (c)). ‘The iodine lines are supposed 
to be very sharp, and any broadening found from these photographs is to be 
attributed to the finite resolving power of the spectrograph along the slit. This 
etfect is negligible near the centre of the fringe system, but becomes gradually 
more important as we go towards the outer fringes, for the dispersion diminishes 
outwards, from an infinite value at the centre. ‘The measured iodine line profiles, 
expressed in wavelength units, are very narrow for the first fringe, but become 
shallower and wider as we go to the outer fringes, the first four fringes giving 17, 
21, 22 and 23mA respectively. ‘The widths to be attributed to the finite vertical 
resolving power of the spectrograph are 3, 7 11 and 18 mA respectively, so that 
the corrected widths become 17, 20, 19 and 14mA respectively. ‘Che tendency 
for the widths to increase outwards appears now to have been eliminated, and this 
has been taken to mean that the vertical instrumental profile of the spectrograph 
affects the definition of the fringe system more than de-focusing on to the spectro- 
graph slit. 

The iodine line width, after correction for the finite resolving power of the 
interferometer, was 14mA. ‘This seems a little too large, since iodine lines are 
frequently used to test resolving powers of the order of 500000 (7). It is to be 
noted, however, that to get measurable absorption lines in the red the temperature 
of the iodine had to be raised as high as 4odeg.C. It is possible also that there 
may have been some maladjustment of the 7-5 mm spacer but, if so, we are 
certain that it was so small as to give a contribution to the solar line widths much 
less than, say, 1omA., 

Summarizing the discussion, the corrections to be applied are as follows: 
1°3 per cent from uniformly scattered light, 1 per cent from grating ghosts plus 
reflection at the inner faces of the doublet spectrograph lens, and 2-6 per cent 
from light scattered from neighbouring channels. As these corrections are 
more or less uniform over a wide spectral range, they could be simply added 
together and regarded as the scattering. ‘The sum, 4:9 per cent, was therefore 
applied to the measured profiles of nine Fraunhofer lines. Other corrections 
depending on resolving power were neglected altogether, since the average 
effective resolving power along the parabolic channels was as high as 6 x 10°, i.e. 
gave a line width about 1omA only. 

5. Profiles, equivalent widths and half half-widths.—'Vhe corrected profiles are 
given in Table I. ‘The lines are arranged and numbered in order of decreasing 
rate of diminution of » with depth (cf. Fig. 2). The table gives the profiles in 
terms of line depths in per cent, instead of residual intensities. Equivalent 
widths and observed half half-widths are in ‘l'able II, both in mA. Explanations 
of the last two columns in this table will be found in section 8. Line depths and 
equivalent widths at the disk centre are compared with mean values from Allen 
(8), in Table ITI. The agreement in equivalent widths appears to be very 
satisfactory. Allen’s line depths, however, are systematically a little larger 
than ours, suggesting that his measures have been slightly over-corrected. 
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TasBie | 
Line depths in per cent in the lines at five positions on the disk 
No. 1. 6252 V RI: —1 LEP : 0:29 eV 
cos @ AA=o 10 20 40 60 to 100 120 140mA 
1°O 12°1 11°6 10°7 8:6 °7 3°0 1°4 o's o"4 
od 1°9 11'S I1°‘2 9g'0 6-0 3°6 1°7 o'6 o'r 
o°6 12°9 12°6 11°6 9g°2 6°4 3°6 I*4 o'9 o'2 
o°4 15°0 14°7 13°6 10°6 ye 4°8 2°5 '°3 o'7 
0°3 4 14°2 13°4 IIo 73) 571 2°6 1°4 o’7 
No. 2. 6330 Cr RI: 1 LEP: o-94eV 
1°O 27°5 26°5 23°4 15'S S8°5 4°2 1°9 o°7 o'4 
o'8 24°3 24°0 22°3 16°3 9°9 $°7 2°7 I'l O'5 
o°6 24°5 23°8 21°8 '7*3 i? 6°7 3°5 1°6 0°6 
O'4 27°2 26°6 24°5 19°3 13°4 Oe i 3°9 I°4 O°3 
0°3 26°4 25°9 24°4 19'I 130 7°6 4°2 2°0 o's 
No. 3. 6393 Fe RI: 3 LEP: 2:27eV 
1°O 16°9 16°3 14°3 100 5°9 *3 1d 0°7 o°2 
o's 16-0 15°5 I4°1 9°6 5°7 3°2 2°0 a3 O*4 
0°6 17°O 16°5 I5‘1 11°5 79 4°5 2°0 O'7 o°2 
O'4 17°0 16°4 4°93 11'S 7°06 4°4 2°! o'9 o'2 
0°3 17°4 17°0 16°0 Fe; g'l ef 9°13 1°6 O°4 
No. 4. 6229 Fe RI: 1 LEP: 2°83 eV 
1'O 35°60 34°5 30°2 19°7 10°7 5°3 1°7 O*4 o'o 
o's 32°85 32°2 29°1 20'S 12°4 5°5 23 o'7 o°2 
o°6 40°97 30°S 28°4 ai°2 13°9 79 3°4 | are) O'5 
O'4 32°8 32°0 29°58 23°1 15°3 9°6 4°0 hej 0°4 
o°3 320 31° 28°8 22°3 15°6 10'O 4'°8 2°6 o's 
No.5 6227 Fe mee LEP: 3:87 eV 
1°O 25°5 24°7 21°4 14°38 S*4 E 1°5 o'6 o°2 
o's 25°0 24°2 22°0 15°4 S°9 4°60 1°6 o°7 o°2 
o°6 24°5 24°0 22°2 16°53 1170 6:1 2°9 o'9 o'o 
O*4 24°6 24°0 22°4 17°6 11°9 7°0 3°6 1S 06 
o'3 23°5 23°3 219 16°8 I1°5 6-9 3°8 r'S oo 
No. 6. 6224 Ni | LEP: 4:09 eV 
ro 25°5 24°5 21°7 151 8°38 4°4 I° °o'9 o'2 
o's 24°6 24°! m7 15°5 9°5 5°3 1°6 o°2 ore) 
0°6 23°3 22°5 20°9 16°1 10°9 5°6 2°4 I°l o°5 
o°*4 23°8 23°5 21°9 17°4 12°3 y he 3°8 1°8 0°3 
0°3 22°4 21°6 20°0 15°5 10°3 5°7 2°6 ore) c'o 
No. 7. 6369 Fe* RI: 0 LEP: 2:88 eV 
1°O I7°2 16°7 14°8 10°5 5°7 2°9 os o'4 o'r 
o's 15°2 14°9 13°9 I1‘o 6°6 ae 1°5 03 0'0 
0°6 14°8 14°6 13°6 10°2 6°3 3°4 I'l 0-3 o'r 
o"4 14°0 13°6 12°8 10°2 6°5 3°4 m7 oF 03 
0°3 1470 13° 12°5 9°8 6°7 4°3 2°2 er O'4 
No. 8. 6238 Fe Ri:-2 LEP: 3:87 eV 
oe A IS ce eS A PT, 
fe dae 5 22457 9°7 4°8 2°2 o'9 
o°6 30°7, 29°90 27'8~—s 22-t—Sssa5*8 10" 5:8 3"1 o"4 
o"4 28-9 28-4 26°9 21°8 16°1 9°6 5°5 26 1's 
0°3 24°3 23°9 22°90 18-6 13°7 8-0 3°0 1°2 0°5 
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TaB_Le [—continued 


No. 9. 6248 Fe RI: 2 LEP : 3°87 eV 
cos 8 AA=o 10 20 40 60 8o 100 120 140mA 
1:0 45'4 44°3 40°5 28°9 18°5 10°7 5°6 371 16 
o's 42°6 41°6 38-2 29°6 19°2 cas 6°4 3°5 iS 
o'6 39°85 39°2 37°0 28°8 20°6 12°6 6°5 3°4 15 
O'4 36°9 36°2 34°3 27°8 19°4 12°! 7°O 3°3 7 
O°3 33°0 32°2 30°3 24°5 17°6 10°9 5°6 2°4 o°6 


Tase II 
Centre limb variations of observed equivalent widths, observed and calculated half half-widths, 
and their effective optical depth for 6310 
W : equivalent width in mA. AA: observed half half-width in mA. AA°: corrected half 
half-width in mA derived from the profile of the absorption coc ficient. 7: effective optical 
depth for AA® at 4 6310 


No. 1. 6252 V RI: I No. 2. 6330 Cr mic 2 
cos 4 W AA AX W AA AX‘ T 
1° 14°7 59°0 55°5 27°38 44°6 40°3 0°70 
o's 15°3 64°5 58°8 28-4 53°8 48°0 o'58 
ob 16-0 62°6 55°8 30°1 58:0 §2°2 O'4I 
o'4 19°3 60°6 54°2 33°6 60°0 50°3 0°25 
o'*3 1g°S 65°6 59°6 33°0 59°5 srs O'17 
No. 3. 6393 Fe RI: o No. 4. 6229 Fe ae 
cos 6 WwW AA AX T W AX AX z 
1° 17°35 46:2 44°5 0°52 34°4 44°5 35°6 0°57 
o's 17°6 47°3 44°1 0°66 35°0 50°6 2°! 0°47 
°o°6 20°1 56°7 51°0 0°48 37°0 54°0 44°0 0°34 
o'4 20°3 57°8 50°2 O31 40°1 56°8 44°58 0°21 
o°3 22°9 61-0 55°0 o'21 40°3 59°0 45°4 o'18 
No. 5. 6227 Fe te No. 6. 6224 Ni me: 2 
cos 6 Ww AA AAS T Ww A AX r 
1°O 25°4 46°9 40°2 o'71 26°4 47°2 411 o'71 
o's 26°4 49°6 43°1 0°56 26-2 52°5 43°9 0°56 
06 28°9 §7°2 48:0 o°41 27°6 58:0 48°6 0°42 
o"4 31°0 59°4 48°8 0°26 30°6 60-4 51°3 0°27 
03 30°0 59°6 46°7 o'18 26°5 60°0 45°6 ors 
No. 7. 6369 Fet RI: o No. 8. 6238 Fe* mi: 2 
cos 4 WwW AA AX T WwW AA AAS T 
1°O 17°7 45°0 44°0 0°82 41°0 54°2 36°6 0°65 
o'8 18-0 55°6 51:2 0°67 40°6 57°8 37°4 0°48 
o°6 17°O 538°9 48°9 0°50 40°2 65°1 42°6 0°36 
O°4 17°O 60°00 52°3 0°33 39°2 66-2 45°6 0°24 
0°73 17°8 600 50°0 0°24 31°9 70°6 46°6 0°17 
No. 9. 6248 Fet Be: Ss 
cos @ WwW AA AAS T 
1° 52°5 g5°5 30°9 0°44 
o'8 52°2 55°9 33°7 0°36 
o'6 52:2 60°9 36°6 0:26 
O'4 49°6 63°4 38°8 o'18 


O73 4371 64°5 37°7 orrl 
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Comparison with Allen at the disk centre 


Line depths Equivalent widths 
Wave- Present Present 
length EF). RI work Allen work Allen 
6223°996 Ni I 25°5% 28%, 26-4 mA 26 mA 
6226-745 Fe I 25°5 29 25°4 26 
6229°240 Fe I 35°6 37 34°4 34 
6238-396 Fe* 2 35°3 34 41°0 44 
6247°569 Fe* 2 45°4 50 52°5 53 
6251°845 Vv —I 12°1 13 14°7 
6330°103 Cr 1 27°5 30 27°8 27 
6369-479 Fe? ° 17°2 14 17°7 7 
6392°545 Fe ° 16°9 13 17°8 14 


In Tables I and II most of the weaker lines show limb strengthening and the 
exceptions are lines from ionized atoms. Going from centre to limb, lines tend 
to become broader and shallower. 

We shall now try to interpret the observed centre limb variation on the basis 
of one of the modern model solar atmospheres, and then will derive the true 
distribution of absorption coefficient with wave-length from the observed line 
profiles using the chosen model. In order to come to any conclusions about 
turbulence we shall need to derive the variation in the selective absorption co- 
efficient over the line profile, because the measured profiles are all appreciably 
affected by saturation. In fact, a simple calculation on the Milne-Eddington 
model shows that, even for the faintest measured lines, with central line depths 
about 17 per cent, the correction for the half half-widths at the centre of the disk 
is of order of 10 per cent (cf. (g)). 

6. Equation of transfer and model solar atmosphere.—Let us divide the equation 
of transfer into two parts, one for the line and the other for the continuum. The 
equation of transfer for the line can be written in the following form: 

dl I : Nz 

= =— =r (2- & B) hvl + gab Nelly, 
where A, B’ and B are respectively the Einstein probabilities of spontaneous 
and stimulated emissions and absorptions, and p, is the normalized function 


giving the form of the line. The equation can be brought into the conventional 
form 





COs gle =I,-E, 
dr, 
if we put 
dr, cos 6= — K,dx, 
pe oS . N,B’ 

“— ab, (B- N, ) hw, 

; , . Ne 
E,= N,A/N, (2 -—_ b’) ° 

N, 


It is to be noted that if we introduce the excitation temperature T., by 


No _ 82 .—hwlk Te 
N, 5° 
orF 


rs) 


~ 
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E, is formally equivalent to the Planck function for 7. Hence 
_ 2h I nrg . 
E,- = (ec -1) PT.) 


Therefore, in any non-equilibrium state, we may write 


dl, se 
cos —- =1],— PT). (1) 
‘The equation of transfer for the continuum can be written 
dl 
cos 80 — =],-—P(T. 
cos 0S" =1,— PT.) (2) 


where 
dr cos 0= —xdx 
and 7), represents the temperature distribution as derived from limb darkening. 

We may superpose these two equations of transfer (1) and (2) by remembering 
that the total intensity /, suffers simultaneously two kinds of absorption « and x,. 
If we denote x,/« by , the resulting equation of transfer can now be written as 

dl, - _ 
cos 0-- =(1+)1,—nP(T,)—PAT,). (3) 

‘This equation cannot be solved unless the distributions of 7;, and 7, are given 
at all depths. As 7’. is normally given by limb darkening, the remaining problem 
to be considered in the theory of line formation concerns 7',. ‘The assumption 

T,, = T,, everywhere 
has usually been called local thermodynamic equilibrium, while the alternative 
assumption 
P(T,)=1 
has been called monochromatic radiative equilibrium. 

‘These assumptions have so far been very useful in many problems, but if the 
problem involves both the uppermost layers, including the chromosphere and 
the deeper reversing layer, either assumption could be made satisfactorily only 
by considering a very complicated stratified atmosphere. In fact, the chromo- 
sphere may make an important contribution, even in the case of weak lines, the 
etfect of which must be considered if we are interested in centre-limb variations. 
In the present state of our knowledge it seems that the best way would be to use 
an empirical distribution of 7\, and to avoid reference to either assumption. 
Although to solve the equation in this way may not be tantamount to solving the 
problem of line formation, it will lead to a sufficient check on the model solar 
atmosphere we wish to use in evaluating the saturation effect. 

In what follows I have used the Bohm—Vitense Model II (10), which as far as 
I can see represents the empirical ¢‘stribution of 7, through the upper layers, 
with minimum deviations from the continuous limb-darkening. Since the 
excitation temperatures naturally vary from line to line, this procedure is justi- 
fiable only from a statistical point of view. In the uppermost layers, for which 
limb darkening measures give us no data, we can equate 7\. to 7, without violating 
any observed limb darkening. In deeper layers, where we have no empirical 
determination of 7\,, the discrepancy between 7, and 7.. might in principle be 
considerable, depending on the degree of departure from thermodynamic 
equilibrium. But in this particular wavelength region the limb darkening 
coefficient happens to be about 1-5, and the difference between “absorption” 


averaged over all solid angles, 


v’ 
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. : r : — 
and “scattering” for weak lines is thus smaller than for other regions. Thus, 
irrespective of the mode of formation, we can safely put 7, = T. in deep layers. 
The equation of transfer (3) now becomes 


cos 0 =(1+)1,-(1+y)PAT,), (4) 
where 7, is the temperature distribution of B6hm—Vitense’s Model I. 

(4) is to be regarded as an empirical equation, rather than as referring to pure 
absorption, although formally there is no difference. P, will not be assumed to 
be linear in r. In the next section exact solutions will be given for equation (4) 
for general (7) and P,(r). 

7. Centre-limb variation of equivalent width.—From equation (4) the intensity 
I, is given by 


1,(0,0)=| P,e~ttre? d(r+7,) sec, 


- 


where 


Tt, = | ndr. 


If we introduce an auxiliary function 
x) 
*(7, 0)= | Pie "°° drsec, 
/ £ 


the intensity outside the line is ®*(0, @). 
The line absorption J, —J, is now transformed as follows: 


I,(0, 4) — 1,(0, 0) = *(0, 8) — | Poe tere? dreccd 


re 
¥ | Poe" ge-*9™ “drecc@ 
“U0 

ao “a 

=*(0, 6) + O*(z, 0) en7 ' - | @*(7, 0) —nsecO)e"8" " dr 
u “0 
*@ 

a | ne . 5 .e7e™* drsecé 
“0 


= | (O*(7,0)— Pier" Jnemte°" drsecd. 


Hence the line depth 7(6) can be written 
[,(0,9)—1,(0,8) _ (” 





(0)= Fh) =), Bt) 9-7 dr, (5) 
where 
g(7, 8) = Oe Gp (D*(7, 8) — Pylr)e *"] sec 6 (6) 


It is to be noted that this expression for the weighting function is identical with 
that derived by Unsdld for general P,(z), but for» <1, (g). It can easily be shown 
that the weight function (6) may also be written as 
o(r, 8) = I [ * dP (7) 
orl (0,0) .4 dr 
‘This last expression gives higher accuracy in numerical computation. It also 


has a definite physical meaning, namely that the weight function at 7 is proportional 
to the observable effective emissivity gradient below r. 





e-? *° “de sec 6 (7) 
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If P,(z) is represented by a linear function 1 + 897, (7) reduces to 


Bo —7 sec 0 
8(7, #)= 1 + Bycosé (3) 
If, moreover, » is small compared with unity, the line depth becomes 
Pa ee lintel 
r(@)= : 1+ ,cos0~ Ui] dr (9) 


This is the expression frequently used by Uns6ld, which is valid only for very 
small line depths. It is the correction factore~*’*’ in our formula (5) that 
measures the degree of saturation with which we are primarily concerned. 

For very large y, the line depth in a “linear”’ atmosphere, corresponding to 
equations (5) and (8), tends to a finite value 





aay . _ bo cos? 
W)> Ts Bocas” (10) 
In a general temperature distribution it becomes 
r (°aP AT) is sece 
1(9) > g(0, 0) cos = D*o a | ie dz (11) 


from equations (5) and (7). It is of interest to compare these two series of 
limiting line depths. ‘The comparison is made in Table IV. 


TABLE IV 
Comparison of limiting line depths 
cos 6 ize) o's 06 O'4 0°3 
Present work 0°88 0°85 o'82 o'81 0°76 
Linear approximation 060 0°55 0°47 0°38 O°31 


The Bohm-—Vitense model gives decidedly greater limiting depths and a much 
less marked decrease in line depth from centre to the limb than does the linear 
approximation. ‘These features would be inherent in every modern model of the 
solar atmosphere in which there is a very sharp decrease to a very low surface 
temperature near the upper boundary. ‘This is apparent in our formula (11) 
according to which the limiting line depth is proportional to the observed effective 
emissivity gradient in the reversing layer. Modern models of this kind there- 
fore give less saturation effects in the lines than an atmosphere based on the 
linear approximation. 

The weight functions (7) were evaluated by numerical integration after 
converting the temperature distribution into the emissivity distribution P,(r) 
for 46310; they are tabulated in Table V for cos @= 1-0 and 0-3 for a number of 
optical depths, and are compared with the weight function (8) for the linear 
approximation with B)=1°5. ‘The comparison shows that in the model we are 
using there is a greater contribution to the line depth from the upper layers than 
there is in the linear approximation, especially near the limb. This can also be 
understood from equation (7) as a consequence of the adopted sharp decrease in 
temperature towards the surface. 

7» was calculated for each line on the same model, using the Saha and Boltzmann 
formulae, partition functions being taken from Claas’ tables (11). Continuous 
absorption coefficients were read from Vitense’s diagrams (12). The resulting 
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TABLE V 


Comparison of weight functions 





Present work Linear approximation 

T6310 cos § =1°0 o°3 1°o o°3 
0°00 0°877 2°525 0-600 1°034 
0°04 772 1°Q17 576 0°905 
orl 57! 0824 543 742 
o'2 500 531 491 531 
O"4 385 320 402 273 
1°O 198 043 221 037 
1'5 0'120 0008 O'134 0°007 
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Fic. 2.—Relative ’s of nine lines. Numbers refer to the line number in Table 1. 


7's for the nine lines studied are shown in Fig. 2. ‘The numbers in the diagram 
refer to line numbers in Table I. The smaller the number the steeper the 
decrease of % with depth. For a number of lines with smaller numbers the »’s 
increase very rapidly at very small optical depths not shown in the diagram. 
The calculated centre-limb variations of central line depths from formula (5) 
are shown in Fig. 3 by broken lines and compared with observation (continuous 
lines). Numbers on each side of the diagram refer to line numbers in Table I. 
Those on the left refer to both observed and calculated curves, those on the right 
to observed curves only. On the whole the agreement with observation seems 
as good as one could hope. It is well established that lines of stronger limb 
strengthening are always associated with steeper 7 distributions. However, if. 
there is an anisotropy or depth-variation of turbulence in the reversing layer, as 
will be discussed later in detail, this comparison can be only qualitative, since the 
limb profiles may be shallower on account of larger turbulent velocities in the 





No. 1, 1957 profiles of faint Fraunhofer lines 15 


upper layers. Centre-limb variations for separate parts of the profile would then 
give no significant test of the adopted model. It is well known, however, that 
centre-limb variations of the equivalent widths are not affected by turbulence, 
provided this is of large scale. ‘Though we do not necessarily postulate that the 
turbulence is of large scale, the comparison of centre-limb variations of equivalent 
widths would be more significant. This is especially true because even the 
strongest lines we measured are still on the linear part of the curve of growth. 
The first step was to determine the wave-length distribution of the absorption 
coefficient for the centre of the disk, and with this then to calculate the expected 
profiles for other positions on the disk. Integration of the calculated profiles 
gave calculated equivalent widths, which were compared with observation. ‘The 
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Fic. 3.—Centre limb variations of line depths 
at the centre of lines. Continuous lines: observa- 
tions. Broken lines: calculations. Numbers on 
both sides of the diagram refer to line numbers. 


Fic. 4.—Centre limb variations of equi- 
valent widths. Continuous lines: obser- 
vations. Broken lines: calculations. Numbers 
on both sides of the diagram refer to line 
numbers. 
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results are in Fig. 4. Once more the numbers on the left refer to the observed 
and calculated values, whereas those on the right refer to observed values only. 
‘The agreement appears satisfactory except for lines of high excitation potential, 
especially those of ionized iron. For two or three of these lines the observed 
equivalent widths definitely exceed the calculated equivalent widths near the 
limb. This discrepancy may possibly be associated with an error in estimating 
the excitation from the ground state. The excitation was calculated by 
Boltzmann’s formula, using the temperature, empirically determined, either from 
the limb darkening or from the central intensities in the visible region. For excita- 
tion from the ground state, however, the effective radiation comes from the far 
ultra-violet, of which we are quite ignorant. In this wave-length region there might 
be an excess of radiation such as was suggested by B6hm (13). But we shall not 
examine this problem in full detail here, for the discrepancy under discussion 
appears not to affect the final conclusion of this work appreciably. We would 
like to believe that the adopted model is good enough to evaluate the saturation 
effect in the observed profiles. 

8. Centre-limb variations of Doppler widths.—We do not know exactly whether 
the turbulence we are going to deal with is predominantly of large scale or small 
scale. In this section we shall determine the true Doppler widths on both of 
these assumptions. 

We shall first derive the Doppler widths assuming small scale turbulence. 
In this case the distribution of the absorption coefficient over the profile can be 
derived directly from the line depths using the known relation between 7(@) and ». 
The absorption coefficient distributions were calculated and half half-widths 
were determined for each line at various positions on the disk. The results are 
in the fourth column of Table II]. The table shows that the widths of all lines 
except that of vanadium increase towards the limb, as was pointed out by Allen. 

The vanadium line 46251-845 is evidently broadened by hyperfine structure 
(14), and we cannot use it for discussing turbulence. The other lines, of Cr, Fe 
and Ni, may safely be assumed to be free of hyperfine structure effects, since 
these elements have even atomic numbers and their primary isotopes have even 
mass numbers (1§). The vanadium line is excluded from what follows. 

If one wishes to get average true widths for each cos @ from Table II, one finds 
values which are more scattered about the mean values than are the observed 
widths. However, if true widths are plotted against the corresponding equivalent 
widths, there is a beautifully linear relation. Accordingly in Fig. 5 we employ the 
more exact procedure of plotting true widths (represented by the larger symbols) 
egainst effective optical depths of line formation. These effective depths are 
calculated by means of the relation 


(7+ 70.5») sec 6 = I 


where 7,,;, means 7, at the part of the profile where 7 has half its central value 7). 
There seems to be no doubt that the true widths again lie on straight lines. 
Since the half half-width of therraal motion is about 23 mA, the rest of the width 
must be due to turbulence. We may conclude that the turbulence velocities in 
a given direction increase with optical depth, and that the velocity at r~0-4 is 
greater in the horizontal direction than vertically. 

At this point one might raise the objection that turbulent velocities apparently 
decreasing with height could have resulted from an overestimate on our part of 
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the saturation effects, as the stronger lines which need more correction for this 
effect proved to have narrower true widths. The answer is given by a comparison 
of 46238 with 46248, both lines belonging to the same multiplet. In spite of the 
tact that 46238 is the weaker of the two, it actually has a greater observed width. 
This would be quite incomprehensible from the point of view of saturation effect, 
unless the true width of \ 6238 is greater than that of 16248. Another example 
of this kind is given by comparing 6369 and A6229. Again the weaker 16369 
has a greater observed width than the stronger (6229. The tendency for the 
line of sight components of turbulence velocities to decrease with height may thus 
be said to be clearly visible in the observed widths themselves. If we take the 
brightness difference of the ascending and descending granules into account, 
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Fic. 5.—Corrected line widths as a function of the optical depth. Larger symbols refer to half 
half-widths of the absorption coefficient, whereas smaller symbols refer to half half-widths derived 
from 1/10-idths. 


the variation of velocities with height would be intensified, as we would naturally 
expect greater contrast in deeper layers. 

One might also suspect that the anisotropy found near t=0-4 is due to the 
brightness difference of the ascending and descending granules (16). In fact, 
if the descending granules had zero brightness, the vertical velocities would be 
measured at about half their real values, supposing no account to be taken of any 
velocity shift of the whole line. But it seems difficult to explain the strong 
anisotropy at 70-4 on this basis, and there seems no doubt that at about this 
level the horizontal motions are predominant. It is quite possible, however, 
that there are ascending and descending granules of diameters about 1ookm 
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with surface brightness ratio as high as 3:1, corresponding to a temperature 
difference of about 1000 deg., and if so the line widths are probably affected. ‘The 
velocity field of turbulence derived in the next section is to be understood as 
provisional. 

If the turbulence increases with depth the lines would be widened in the 
wings more than at the centre, since the wings are formed effectively at a greater 
depth. To check this the line widths were computed for a point where the 
absorption coefficient is 1/10th that at the line centre. ‘The damping parameter 
a can be taken to be near zero for these faint lines and then the ratio of 1/10-width 
to half-width is 1°82. The 1/10-widths were therefore reduced by this factor 
to give fictitious half-widths corresponding to the derived 1/10-widths. ‘The 
fictitious half half-widths were plotted in Fig. 5 (small symbols) against effective 
depths given by the relation 


(7 +79.,)secO=1 


where 7,.,, means 7, at the point of the profile where y»=0-1m. This diagram 
tends to confirm the distribution of velocities with depth suggested by the half 
half-widths. Were a taken to be about o-1 instead of zero, the widths would be 
reduced by about 3 per cent and the fit in the diagram would not be appreciably 
altered. But in view of the low accuracy of the 1/10-widths this argument is 
not very significant. 

On the assumption of large scale turbulence we again obtained about the 
same velocity distribution with depth as is shown in Fig. 5. ‘The widths were 
calculated in the following way. For simplicity, lines at the centre of the solar 
disk were assumed to be broadened by small scale turbulence only. This 
assumption is plausible in view of the fact that the line widths there are consistent 
with the small scale turbulence velocities derived from the curve of growth by 
Wright (17), and ten Bruggencate and Houtgast (18). From these profiles at 
the disk centre, profiles at other positions on the disk were computed using the 
previous calculations of the centre-limb variation of the absorption at each point 
of the line profile. The calculated profile widths were compared with those 
observed and the amount of blurring necessary to bring them into agreement 
gives the velocities of the large scale turbulence. However the variation of 
velocity with depth which was suggested by the half half-widths is not apparent 
in the observed line profiles. Velocities derived from the 1/10-widths are about 
the same as the velocities from the half half-widths of the same lines, in spite of 
the difference of average depth of the layers in which they originate. If we 
take damping into consideration the situation becomes worse. 

Although we do not attach great weight to the 1/10-widths, we are inclined 
to believe that the turbulence responsible for the line widths is predominantly 
of small scale. It does not necessarily follow, however, that the small scale 
turbulence suggested has to be much smaller in size than the scale height. The 
small scale turbulence we have in mind can be as big as the scale height and 
accordingly may have a free path comparable with this. This can be under- 
stood if one bears in mind that turbulence of this size also includes much smaller 
turbulence and that the greater part of the effects of “small scale” turbulence 
comes from the smaller turbulence elements (19). 

Wright (17), found 1-okm/sec for the velocity of the small scale turbulence 
from the curve of growth at the centre of the disk. The corresponding half 
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half-width after adding the thermal velocity is 28mA at 46300. ‘This is to be 
interpreted as the velocity at the depth where lines lying on the Doppler part of 
the curve of growth are formed. These lines have equivalent width about 130 mA. 
Since the equivalent width of our strongest line is 52 mA, the effective depth of 
the Wright velocity should be much higher in the photosphere than t=0-4. It 
seems probable from Fig. 5 that our vertical velocity is very near the Wright 
value there and our assumption that the turbulence is of small scale thus does not 
contradict Wright’s results. But it is evident that the larger value of small 
scale turbulence found by ten Bruggencate and Houtgast (18), viz. 1-7 km/sec, 
is inconsistent with our profiles, since the corresponding half half-width in Fig. 5 
is this time 35 mA. 

g. Variation of turbulence with depth in the reversing layer.—The velocity 
field of turbulence at various depths was derived from Fig. 5 after subtracting 
the thermal velocities for the corresponding temperatures, and this is shown in 
Fig. 6. In this diagram the most probable turbulence velocities €, are drawn 
separately, according to whether they are ascending or descending, and the 
optical depths for A5o010 are also shown. Strong asymmetry around 7~0°4 is 
again evident here. In deeper layers, where at present we have no means of 
measuring horizontal velocities, we do not know whether, for example at 7~1, 
the turbulence is predominantly vertical or horizontal, or isotropic. At any 
rate it looks as though the upward motion at r~1 of some of the granules escaping 
from the top of the convection layer into the radiative equilibrium zone are 
deflected into horizontal motions at t~o-4 and that these are finally converted 
into downward motion. As the speed of vertical motion at r~1 is about the 
same as of horizontal motion at 7~o-4, there should be no difficulty in satisfying 
continuity conditions. 

Allen’s anisotropy (1) of 1°74km/sec and 2-79 km/sec for the vertical and 
horizontal motions respectively now have to be interpreted as the velocities at 
7=1 and say 70-2 respectively, now that we find evidence for a strong variation 
of turbulence with optical depth within the reversing layer. Our velocities, 
however, appear to be somewhat lower than the corresponding Allen values, a 
difference which may be associated with differences in the way of correcting for 
saturation effects. 

Quite recently McMath et al. (20), have found increasing large scale turbulence 
towards the limb from the zigzag appearance of A 4626, Crl, and 5853, Ball. 
They found that €, increases from 0-6 km/sec at the centre of the disk to 0-8 km/sec 
at the limb. ‘This points to a variation of velocity of the kind found by Allen 
and supported by the present work. We do not know, however, whether the 
McMath results, which undoubtedly refer to much higher layers than ours, are 
directly related to our velocity field. 

The turbulence variation with optical depth proposed here is the opposite 
of that suggested by Su-Shu Huang (21), for F type stars and also the Sun. 
Su-Shu Huang found a strong increase of turbulence with height in F-type 
giants and a very weak increase in the Sun. But as far as the Sun is concerned 
his argument is not very convincing, since he based his conclusions on Allen’s 
profile measures. We do not believe that the turbulence/height relation need 
necessarily be the same in F giants and G dwarfs. In the F-type stars the con- 
vection layer must be located much higher than on the Sun, and may include the 
whole reversing layer, whereas in the Sun there is a layer in radiative equilibrium 
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Fic. 6.—Angular distributions of most probable velocities of turbulence at various depths. Con- 
tinuous arrows are observed velocities, whereas broken arrows are extrapolated values. 
above the convection layer, acting as a barrier to the granules ascending from 
below. What we see in the F-type stars is the increasing turbulence with height 
which is associated with decreasing density. But what we see in the Sun 
is only the tail end of convection prevailing beneath r~1. _‘It is understandable 
that the granular elements surging into the reversing layer are decelerated or 
deflected, or diminish in number as they ascend, giving a decreasing turbulence 

with height. 

Summing up, we may conclude that the vertical motions of the granule 
elements from the convection layer are gradually deflected and converted into 
horizontal motions as they penetrate the radiative equilibrium zone, and are 
finally converted into downward motions back to the convection layer, without 
much change in their original speed of about 2km/sec. ‘The size of the granules 
would be about the scale height, and the mean free path would be the path length 
of the circulation just described. No substantial part of the granulation energy 





No. 1, 1957 profiles of faint Fraunhofer lines 21 


is transferred to the reversing layer. The turbulence observed in the highest part 
of the reversing layer and in the chromosphere is to be regarded only as a secondary 
disturbance arising from a small fraction of the granulation energy. ‘This 
secondary turbulence will probably increase with height, due to the decreasing 
density. 

The author has great pleasure in acknowledging with thanks his indebtedness 
to Professor F. J. M. Stratton, President of Commission 38, for a grant per- 
mitting him to stay at the Cambridge Observatories, and to Professor R. O. 
Redman, for encouragement and for the facilities put at the author’s disposal, as 
well as for reading and revising the manuscript of this paper in detail. He is 
especially indebted to Dr H. von Kliiber, without whose active guidance and 
valuable discussions this work could not have been done. 
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Summary 


The interference spectrum of the brighter lines in the chromosphere of the 
uneclipsed Sun has been photographed by means of the Oxford 35 m solar 
telescope, used in conjunction with a Fabry-Perot etalon of high reflectivity 
and large range, and a low dispersion spectrograph. The very narrow 
spectroscope slit was set tangentially to the 32 cm solar image. Under these 
conditions, the emission lines of hydrogen, helium and calcium appear as a 
system of interference fringes, with a resolving power of nearly 10°. 

A method is developed for the quantitative photometric analysis of the 
fringe system, from which it is possible to determine chromospheric line- 
profiles as a function of height, and the chromospheric intensity gradient. 

This technique has been employed to determine the profile of the D, line 
ofhelium. It is found to be nearly independent of height over the range 1500- 
4250 km, and, after correction for the effect of fine structure and photo- 
metric broadening, corresponds closely to a Gaussian distribution of the 
form I= I, exp | —(8A/0°34)?]._ These results are interpreted as implying that 
there is little self-absorption, and that the atoms which give rise to the line are 
subject to turbulent velocities of the order of 17 km/sec. 





1. Introduction.—The application of interferometric techniques to solar 
studies, initiated by Fabry and Perot in 1902 (6), received an impetus some years 
ago, when it became possible to exploit the progress in the production of high 
quality silver and dielectric reflecting films and narrow band interference filters 
(12). 

In the field of photospheric interferometry, the observational achievements 
of the past few years have been of considerable significance. Investigations of 
the mean solar magnetic field by Thiessen (18) and by von Kliiber (25), the 
study of the Einstein shift by Adam (1), the fine structure analysis by Hindmarsh 
(7), and Kinman’s search for solar deuterium (9) are examples from a varied 
field of interferometric solar researches, each of which has made a definite con- 
tribution to outstanding problems in solar physics. 

Nevertheless, as the writer pointed out (20), the small range between successive 
orders makes the Fabry-Perot interferometer a somewhat difficult and limited 
instrument for use in the study of an absorption spectrum, and it is largely for 
this reason that it has been overtaken by the grating. In studying the emission 
spectra of the chromosphere, of solar prominences, and of the corona, these 
limitations are not encountered, and the great versatility of the etalon commends 
it in various ways. In it we have a compact spectroscope whose dispersion and 
resolving power can be chosen at will. With relatively high Juminous efficiency, 
it combines the possibility of simultaneous study of line structures and 
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displacements over extended areas (a fact made use of by Jarrett and 
von Kliiber (8) in their coronal interferometry), or of recording on a single plate 
high resolution interference spectra of the whole of the photographically 
accessible spectrum. ‘The potential value of such an instrument for the simul- 
taneous study of chromospheric line structures, or for using to the best 
advantage the rare and brief moments of total eclipse, is sufficiently evident. 

Attempts to realize this promise of interference spectroscopy in the study of 
solar emission spectra have in fact almost all been made under the difficult and 
chancy conditions of total eclipse. A number of early attempts from Sproul 
Observatory (26) to obtain coronal fringes failed, as the writer has pointed out 
(21), because of defective experimental design. ‘The success of Jarrett and 
von Kliiber’s coronal interferometry (8) in the eclipse of 1952 was due entirely 
to the highly controlled conditions of their observations. On the other hand, 
Briick and Jackson in their attempt (5) to secure an interference spectrum of the 
chromospheric lines during the eclipse of 1952 February 25 did not succeed in 
photographing fringes, since the line width, to which from their observations 
they were able to assign a minimum value of 1-2cm~! (0:24 A at 4471), exceeded 
the range of the etalon. 

Now the stronger chromospheric lines, and the lines in solar prominences, 
can be observed outside of eclipse under favourable conditions by setting the slit 
of the spectrograph across or tangentially to the solar limb so that it intersects 
the chromosphere and the overlying prominences. ‘The possibility of developing 
and testing an interferometric technique in the more tranquil atmosphere pre- 
vailing outside of eclipse is clearly worth exploring. ‘This fact is emphasized 
by the success of Harrison’s recently invented echelle grating in obtaining high 
order interference spectra of prominences outside of eclipse (13). ‘The fact that 
the chromospheric emission is confined to a narrow layer of less than 10” of arc 
immediately bordering on the solar limb naturally makes the problem of chromo- 
spheric interferometry observationally more difficult, since quite small excursions 
of the Sun’s limb, due to scintillation or guiding errors, and light scattered 
atmospherically or instrumentally from the photosphere, frequently produce an 
intolerably heavy scatter spectrum. 

In the summer of 1955, the new Oxford 35 m solar telescope (16) had been 
completed, but awaited the installation of its large spectrograph, the place of 
which was occupied by the Hilger medium glass prism spectrograph which had 
been used with this telescope to obtain material for Bishop’s study of the latitude 
variation of D, (2). The dispersion of this instrument is about 40 A/mm at 
5876A. ‘This arrangement provided the writer with an opportunity to attempt 
to develop a technique for photographing the interference spectrum of the 
chromosphere outside of eclipse and for analysing the fringe systems in order to 
obtain from them the line profiles as a function of height, and the chromospheric 
intensity gradient. 


2. Photography of the chromospheric interference spectrum 

2.1. General arrangement.—The writer has discussed elsewhere the factors 
determining the best arrangement of an etalon in conjuncuon with a spectrograph 
for solar studies (20). When absolute wave-lengths are not required, and when 
the image-forming properties of the spectroscope are limited only by slit width 
and plate grain, as in the Hilger instrument, the arrangement of the etalon, in a 
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collimated beam before the slit, and that in which it follows the formation of the 
spectrum, are equivalent. The choice of the former arrangement was necessitated 
by the strongly inclined focal plane of the Hilger Spectrograph. 

The general arrangement of the system is shown in Fig. 1. The solar image, 
32cm in diameter, was formed on a screen A containing a rectangular slot, which 
was re-imaged on the slit of the spectrograph by means of the field lens, collimator 
and ring-lens system BCD. The slot was so related to a circle ruled on the 
screen, of the same diameter as the solar image, that when the solar image was 
guided on the circle the slit was tangential to the re-imaged solar limb, with the 
base of the slit at the point of contact. The slit was about 2cm in height. Under 
these conditions, the full height of the slit corresponds to an elevation of 8” above 
the solar limb. ‘The guiding was controlled during exposures (of one to two 
minutes’ duration) by varying the orientation of a plane parallel guiding plate F. 
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Fic. 1.—Experimental arrangement for chromospheric interferometry. 


‘The etalon E was placed in the collimated beam. The Hilger plates, 2cm in 
diameter, had a wedge angle of 15’, and were so oriented that the two principal 
first order ghosts produced images of the slot in the focal plane of the slit, to 
either side of the principal image, so that no scattered light from these ghosts 
could enter the system. ‘The etalon was tilted a few degrees about a horizontal 
axis, so that the centre of the fringe system fell vertically below the slit. In 
4 way, several chromospheric fringes could be formed along the length of the 
slit. 

2.2. Some considerations of experimental designs.—The success of this kind 
of work depends greatly on the right choice of slit width, etalon spacer and 
reflectivity, and of the focal lengths of the field lens, collimator, and ring-lens 
system. 

2.2.1. Slt width.—If there were no scattered light from the solar limb, it 
would be possible to use a wide slit, since the only effect of so doing, in the case 
of an isolated emission line, would be to admit a longer arc of each fringe, without 
altering its intensity. In fact, however, the emission fringes will be accompanied 
by a continuous photospheric scatter spectrum. It is obvious that the intensity 
of a continuous spectrum is directly proportional to slit width. On this general 
ground one would use as narrow a slit as possible, the useful limit being set at 
0°03 mm by plate grain and light loss by slit diffraction. In fact, the scatter 
spectrum is not continuous, but takes the form of a channelled spectrum. How- 
ever, for all but the very narrowest slit widths, the channels overlap and obey 
the same intensity law as a continuous spectrum, and even when they do not the 
use of a narrow slit will, as Shane has shown (17), reduce scatter from adjacent 
orders. We conclude, therefore, that a very narrow slit is essential. 
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2.2.2 Etalon range.—It is known that the strong chromospheric lines are 


greatly broadened. Ha extends over sure than 2A, D, over more than 1A. 
The range of a Fabry—Perot etalon is giver |» 
pe =A?/2i (1) 


where A is the wave-length, ¢ the spacer length. We can therefore only use 
spacers of less than 1 mm for Ha and less than about 2mm for D,. There is one 
significant qualification on this limit. If the slit width is sufficiently small for 
adjacent parabolic channels to be separated, it remains possible to resolve adjacent 
orders without overlap. An interesting example of this is illustrated in Plate 2(c) 
in which 3 orders of Hf are clearly distinguished by their relation to the channel 
system, although they just exceed the etalon’s range. 

2.2.3. Resolving power.—The resolving power of an etalon of reflectivity r, in 
order n is given by 

A/dA = [2-98r"?/(1 —1)]n. (2) 
Having fixed t by considerations of range, we can increase the resolving power only 
by increasing r. However, since plate irregularities and defective parallelism 
may be expected to broaden the fringes to 1/30th of an order, the maximum 
resolution we can expect under working conditions with reflectivities exceeding 
0-92 is about 30m, where n is the order of interference, that is to say 100000 at 
6000 A for a 1mm spacer, a value that will not be increased by increasing the 
reflectivity. 

2.2.4. Dispersion.—In order that this resolving power of 10° may not be lost 
by plate granularity, the dispersion should be great enough for the resolving 
limit, 0-06 A at 6000 A, to exceed the dimensions of plate grain (say 0-03 mm), a 
condition that is well satisfied by a dispersion of 1 A/mm. Since the dispersion 
of the etalon is given by 

D=f,/Aa mm/A, (3) 
where f/f, is the focal length of the ring lens and a the tilt of the etalon in radians, 
we can obtain the desired dispersion by various combinations of a and f,. 

In order to preserve the solar image size on the one hand, and on the other to 
preserve its intensity, the focal lengths of C and hence the field lens B, should be 
more or less equal to that of D. But the focal length of the field lens determines 
the size of the image of the solar telescope objective, which is formed on the etalon, 
and which must be kept sufficiently small to ensure that there is no loss of beams by 
multiple reflection. A 50cm field lens already forms an image 0-7 cm in diameter 
on the etalon. Accepting this indirect upper limit to the focal length of D, we 
obtain the required dispersion for tilts not exceeding 5°. In this way the optical 
constants of the system are completely determined. 

2.3. Observations.—To satisfy the preceding requirements, the spectrograph 
slit was reduced in width to 0-02 mm, which, however, on account of the tilt of 
the focal plane, corresponds to a slit image width of 0-o4mm. _Etalon spacers of 
invar steel, 1 mm and o-5 mm in length, were employed. ‘Through the kindness 
of Dr Kuhn, silver coatings were deposited shortly before use, with mean measured 
reflectivities and tranmissions of 95-5 and 2-5 per cent respectively at 5500 A, 
giving a theoretical resolving power, for the 1 mm spacer, of 260000. On each 
plate of the chromosphere, control fringes were photographed using the 6438 A 
line from a Cd discharge. The actual resolving power deduced from micro- 
photometric measurements of the profiles of the Cd fringes, was 90 000, in good 
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agreement with the indications of the previous section. The lenses C and D were 
visual achromats of 45°7cem focal length. ‘The field lens B was a spectacle lens 
of 50cm focal length. A stop, only slightly larger than the image of the objective, 
was placed immediately in front of the etalon, to eliminate scattered light and to 
ensure that the Cadmium lamp, when interposed in the space between B and 
C, illuminated the same area of the etalon as the Sun. ‘The Cadmium lamp was 
also used to adjust the etalon for parallelism in situ and to orientate the etalon with 
respect to the slit. 

A plate normaily contained five exposures. ‘I'wo of these were of the chromo- 
sphere. ‘The third was an exposure made through a rhodium step wedge cali- 
brated previously in the D region by Bray (3), using the solar centre as source. A 
fourth exposure was made without the wedge, as a test for uniformity of illumina- 
tion along the slit, again using the solar centre. ‘The interferometer was left in 
position for these exposures. All of these exposures were made with the narrow 
slit. The final exposure, superimposed on the fourth, and with a wide slit, was 
the Cd lamp exposure already referred to. In addition to serving as an adjust- 
ment check, it provided a system of sharp fringes from which the centre of the 
chromospheric fringe system could be determined. Plate 3(a) shows a greatly 
enlarged portion of this fringe system, superimposed on the interference spectrum 
of the solar centre. Exposure times were equalized by means of neutral filters at A. 

During the summer of 1955 a successful series of interference photographs 
was taken with a separator of 1 mm (23). These showed fringe systems of 
Ha, Hf, and of the H and K lines. The most perfect fringe system however 
was that shown by the D, line of neutral helium, as was to be expected in view of 
the ratio of the line width to the etalon range. Plate 2(a) and (5) show enlarge- 
ments from this fringe system. It has, superimposed upon it, a strong photo- 
spheric interference spectrum of scattered light. On account of the great 
intrinsic width of the line, the fringes are extended in the direction of the etalon’s 
dispersion, while being narrowly confined perpendicularly to this direction on 
account of the very narrow slit. Hence they present (in negative prints) the 
appearance of black stitches. Plate 2(b) shows a still greater enlargement, in 
which the channelled spectrum and the plate grain appear. 

A careful inspection of the Ha fringes made it clear that one was observing 
slightly overlapping orders of a heavily self-reversed line. ‘To confirm this 
deduction, a further series of plates was obtained in the spring of this year, using 
a separator of o-5mm and Kodak IIN infra-red plates, in place of the Ilford 
R.40 plates used in the previous summer, since the sensitivity of the latter plate 
falls off rapidly at Hx. Plate 3(5) shows the self-reversed Ha line now com- 
pletely free from overlapping of orders. The reduction in the scattered light is 
due largely to the realuminising of the telescope mirrors during the intervening 
winter. The sensitivity of the IIN emulsion extends far into the infra-red, 
so that the spectrum of the solar centre shows the atmospheric A band near 7600 A. 
With the o-5 mm spacer and the increased wave-length, the range of the etalon is 
5°8 A, so that the channels are still well resolved and reveal the structure of the 
band in an unexpected way (Plate 3(c)). 

3- Photometric analysis of the interference fringe system.—We shall now proceed 
to consider how to interpret quantitatively the chromospheric fringe systems just 
described. It is evident that the intensity distribution along any fringe in the 
direction of etalon dispersion is closely related to the profile of the line that 
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(a) The superimposed system of cadmium fringes. 
(b) Ha with a separator of 0-5 mm showing self-absorption. 
(c) Bands of atmospheric oxygen (6870 and 7600A) in the interferometric photospheric spectrum. 
(Negative prints.) 
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produces it, and that the falling off in overall intensity from fringe to fringe is 
closely related to the chromospheric intensity gradient. However, since the 
Sun and the fringe system are imaged in the same plane on the slit, different 
points on a single fringe correspond not merely to different wave-lengths within 
the line, but also to different heights in the chromosphere. What is strictly 
true however is that if we denote the position of any point in the fringe system by 
x, where x is measured from the edge of the spectrum which is coincident with the 
Sun’s limb, then to each x will correspond an intensity /,, of which the plate 
blackening is a measure, which refers to the intensity of the chromosphere at a 
definite height /,, and at a definite wave-length A,._ The problem of interpreting 
the fringe system, therefore, reduces itself to two stages. First, we must deter- 
mine J, h, and A for a series of closely adjacent values of x covering the whole 
fringe system. ‘This part of the problem alone involves the experimental 
technique. ‘The second part of the interpretation is to use the purely solar 
data, J, h and A, so obtained, to construct the true profile at any given height, and 
to measure the intensity gradients. We shall consider these stages in turn. 

3.1. Determination of height.—Provided the origin of x measurement is 
truly the point of contact of the tangential slit with the limb of the Sun, a point at 
distance x along the slit corresponds to a height A above the limb, such that 


x*=2hR (4) 


where R is the radius of the solar image. A may therefore be determined at once. 
It must however be remembered that the height / is the height down to which 
the line of sight from the observer descends into the solar chromosphere, so that 
the light from x is produced by the integrated contributions of all heights 
above h. 

3.2. Determination of wave-length—The problem of determining wave- 
lengths from one side only of a Fabry-Perot ring system is treated in detail by 
Tolansky (1g). ‘The following simple treatment is quite adequate for the present 
case where only a wave-length scale is required. 

The dispersion formula (3) may be written in the form 


d 
5 =f?/Ar (5) 


where f is the focal length of the collimator and r the radius of the fringe. 


Integrating (5) we obtain 
A= kyr? + ky. (6) 


To determine the wave-length corresponding to a point x we must determine 
r for the point x, and evaluate k, and k,. Since the centre is not on the plate, 
r cannot be measured directly. The position of the centre may however be 
determined approximately. If, choosing the mth fringe as origin, we measure the 
coordinates R,,,,, R,,.. of the fringes n+1, n+2, it may easily be shown from 
the geometry of the fringe system that the radius of the nth fringe is given by 
~— 2Rin+1°— Rais® 4 (7) 
; Rn+2—2R,, +1 
In this way the centre may be determined in relation to the origin of x (preferably 
by measurements on a sharp fringe system such as that provided by the Cd lines), 
and the values of r computed for each value of x. Now consider two wave- 
lengths A, and A, +, where p is the wave-length range of the etalon. Denoting 
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the radius of a fringe of the former wave-length by r,, and the next fringe of the 
same wave-length by r., we have, by (6) and the definition of range, 
A, =k? +k, 
A, +pm=hyr? +h, 
whence 
ky =p /(r2? — 1?) =4/(r,, 4° -1,”). (8) 
Since the difference of square of successive radii is, for a given wave-length, a 
constant of the system, it is preferable to use a suitably determined mean value 
from several fringes to compute k,. ‘This is the mean value indicated in the 
denominator of the last term in equation(8). ‘To determine k, we need only 
identify and determine the radius 7, of a fringe of known wave-length A, (e.g. the 
centre of the D, line), yielding 
ky =Aq— Ry”. (9) 
Hence, finally, in any fringe 
(AA) = MP? = 167) ga? — Tn”) (10) 

3-3. Determination of the chromospheric intensity.—Consider the micro- 
photometer tracings obtained by running a very short and narrow microphoto- 
meter slit across the fringe system at each value of x for which we require a 
determination of J,, the chromospheric intensity at x. "The appearance of this 
tracing will resemble Fig. 2, which is reproduced from an actual record. ‘The 
tracing shows the density variations in the channelled spectrum of the scattered 
light, together with sections of the absorption line profiles, as in Shane’s method 
(17), of which this is in fact an application to an emission line. ‘The section of 
the D, fringe traversed in each run appears as a peak superimposed on the scatter 
spectrum. When this density curve is converted to an intensity curve, by means 
of the wedge calibration, the area under the section of the D, fringe, discounting 
the scattered light contribution, is a measure of the chromospheric intensity. 
It would be very laborious to have to reduce each trace in this way. Bishop has 
shown (2) however, from high magnification tracings from direct spectra of the 
D, line taken with the same spectroscope and measured with the same slit as that 
here employed, that the area required and hence the intensity of the chromosphere, 
bears a linear relation to the intensity difference between the peak and the level 
of scattered light. This result is not surprising, since it reflects the fact that 
the Moll microphotometer with a narrow slit (15) does not appreciably distort the 
profile of the D, peak, the form of which is constant, since it is determined almost 
wholely by spectroscopic slit width and plate grain. 

3-3-1. Corrections.—The simple method so described is sufficient when the 
scatter spectrum is free from absorption lines in the region of the chromospheric 
line profile. If these are present, the profile can be determined from regions 
of the scattered light spectrum where the fringe system of the chromospheric 
line is absent. After determining the wave-length of the chromospheric section 
traversed in any run, one would then know what percentage of the scattered 
light adjacent to the chromospheric fringe, and corresponding to the continuous 
background, should be subtracted from the peak intensity to give the chromo- 
spheric intensity. 

If there is any variation in instrumental transmission intensity along the 


spectroscope slit, this must be determined, and an appropriate correction to the 
intensity determined for each value of x. 





sc nheoe aecene aie Keanna 


No. 1, 1957 Chromospheric interferometry 29 


The intensity so determined is of course proportional to the flux received 
from the chromosphere over a short wave-length interval, determined here by 
the length of the microphotometer slit, and the interferometric dispersion. This 
latter quantity varies according to equation (5), and it would appear at first sight 
as though one should apply a correction for this, such as must be applied to 
allow for the variation of dispersion in the microphotometry of prismatic spectra. 
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Fic. 2.—Typical microphotometer trace across a chromospheric fringe. 





Further consideration of this analogy indicates that this is not the case. For it is 
a characteristic of the Fabry-Perot instrumental pattern, in contrast to the slit 
pattern of a spectroscope, that its width remains always a constant fraction of its 
range: consequently its half width in mm on the plate decreases with dispersion, 
from fringe to fringe, the monochromatic peak intensity remaining constant. 
The just analogy should therefore be a prism spectroscope whose slit width is 
made to vary directly with the dispersion, The writer has analysed this question 
quantitatively elsewhere (22). 
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the radius of a fringe of the former wave-length by r,, and the next fringe of the 
same wave-length by r,, we have, by (6) and the definition of range, 
A, =hyr? +k, 
A, +p= kyr.” + ky 
whence 
k= p/n —rP)=4/(r, 2° -1,")- (8) 
Since the difference of square of successive radii is, for a given wave-length, a 
constant of the system, it is preferable to use a suitably determined mean value 
from several fringes to compute k,. This is the mean value indicated in the 
denominator of the last term in equation(8). ‘To determine k, we need only 
identify and determine the radius r, of a fringe of known wave-length A, (e.g. the 
centre of the D, line), yielding 
ky =Ay — hy? (9) 
Hence, finally, in any fringe 
(A~Ag) = MP? = 16)? — Tn”) (10) 

3.3. Determination of the chromospheric intensity.—Consider the micro- 
photometer tracings obtained by running a very short and narrow microphoto- 
meter slit across the fringe system at each value of x for which we require a 
determination of J,, the chromospheric intensity at x. ‘The appearance of this 
tracing will resemble Fig. 2, which is reproduced from an actual record. ‘The 
tracing shows the density variations in the channelled spectrum of the scattered 
light, together with sections of the absorption line profiles, as in Shane’s method 
(17), of which this is in fact an application to an emission line. The section of 
the D, fringe traversed in each run appears as a peak superimposed on the scatter 
spectrum. When this density curve is converted to an intensity curve, by means 
of the wedge calibration, the area under the section of the D, fringe, discounting 
the scattered light contribution, is a measure of the chromospheric intensity. 
It would be very laborious to have to reduce each trace in this way. Bishop has 
shown (2) however, from high magnification tracings from direct spectra of the 
D, line taken with the same spectroscope and measured with the same slit as that 
here employed, that the area required and hence the intensity of the chromosphere, 
bears a linear relation to the intensity difference between the peak and the level 
of scattered light. This result is not surprising, since it reflects the fact that 
the Moll microphotometer with a narrow slit (15) does not appreciably distort the 
profile of the D, peak, the form of which is constant, since it is determined almost 
wholely by spectroscopic slit width and plate grain. 

3-3-1. Corrections.—The simple method so described is sufficient when the 
scatter spectrum is free from absorption lines in the region of the chromospheric 
line profile. If these are present, the profile can be determined from regions 
of the scattered light spectrum where the fringe system of the chromospheric 
line is absent. After determining the wave-length of the chromospheric section 
traversed in any run, one would then know what percentage of the scattered 
light adjacent to the chromospheric fringe, and corresponding to the continuous 
background, should be subtracted from the peak intensity to give the chromo- 
spheric intensity. 

If there is any variation in instrumental transmission intensity along the 


spectroscope slit, this must be determined, and an appropriate correction to the 
intensity determined for each value of x. 
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The intensity so determined is of course proportional to the flux received 
from the chromosphere over a short wave-length interval, determined here by 
the length of the microphotometer slit, and the interferometric dispersion. This 
latter quantity varies according to equation (5), and it would appear at first sight 
as though one should apply a correction for this, such as must be applied to 
allow for the variation of dispersion in the microphotometry of prismatic spectra. 
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Fic. 2.—Typical microphotometer trace across a chromospheric fringe. 


Further consideration of this analogy indicates that this is not the case. For it is 
a characteristic of the Fabry-Perot instrumental pattern, in contrast to the slit 
pattern of a spectroscope, that its width remains always a constant fraction of its 
range: consequently its half width in mm on the plate decreases with dispersion, 
from fringe to fringe, the monochromatic peak intensity remaining constant. 
The just analogy should therefore be a prism spectroscope whose slit width is 
made to vary directly with the dispersion, The writer has analysed this question 
quantitatively elsewhere (22). 
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3.4. Line profiles.—By the methods of Sections 3.1-3.3 we may extract all 
the solar information that the chromospheric fringe system will yield, in the 
form of tables of associated values of h, J and A for successive values of x. From 
these values, we may obtain true line profiles at required heights, by a procedure 
suggested by Plaskett, which we shall call the method of monochromatic in- 
tensity gradients. 
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Fic. 3.—Derivation of line profiles from chromospheric fringes. 
(A) Fringe profiles; (B) Monochromatic intensity curves; (C) Resultant profile. 


To carry out this method, we first plot from our tabular values the fringe 
profiles, with ordinates J, abscissa A. Fig. 3(A) presents such a plot obtained 
from one of our D, fringe systems. Now mark off from the line centre Ay a 
series of equidistant ordinates corresponding to A,, A, etc. Next read from the 
fringe profiles the intensity J at Ay on each fringe available, and obtain the corres- 
ponding heights h,,, A,, h, etc. by interpolation from the table. For each A in turn, 
we,can now plot height as a function of intensity. This stage is shown schem- 
atically for 4 wave-lengths in Fig. 3(B), and these curves give essentially mono- 
chromatic intensity variation with height in the chromosphere. From this array 
of curves, we can now read off the monochromatic intensity at any desired height, 
for the series of wave-lengths Ay, A,, A, and in this way construct the profile at 
any height represented on all the graphs. 

The profiles at different heights are therefore obtained independently from 
each plate. In combining them, profiles at some common height must be 
normalized, since the absolute intensity will vary from plate to plate owing to the 
varying conditions of exposure. 

3-5. Chromospheric intensity gradients.—Having obtained the profiles at a 
series of heights from individual plates, it is clearly easy to determine the integrated 
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chromospheric gradient for the line as a whole by plotting the area of successive 
profiles from the same plate as a function of height. ‘The mean gradient adopted 
must of course be the mean of the gradients so determined from individual plates. 


TaBLe I 
Computed and observed profiles of Ds 








A C, Cy C; C O O-C a 
5875°140 141 00°0 141 16°1 17°6 brs £35 
"190 20°9 00"! 21°0 23°% 26°8 +3°6 2°0 
"240 29°5 00°! 29°6 31°8 360°5 -4°7 1°8 
"290 40°0 00°2 40°2 42°0 42°4 +O"4 1°3 
"340 51°9 00°4 52°3 53°6 53°5 orl 2°3 
"390 64°4 00°7 65°1 65°3 65°0 —0°3 1°4 
"440 76°6 Or'! 79°79 76-9 75°! —1°8 1°4 
"490 87°3 o1'8 89°1 37 °1 84°5 2°6 3 
"540 95°2 02°6 97°8 94°5 2°32 ke 1°9 
‘590 99°7 03°7 103°4 98°3 97°8 -O°5 16 
“640 99°7 05°0 104°7 100°0 100°O fohze) 03 
‘690 952 06°5 IOI°7 97°7 07°4 -0°3 2°3 
"740 87°3 O81 95°4 g2°2 93°6 +34 a°7 
‘790 76°6 09°5 86:1 84°3 86°1 1°8 1°4 
"840 64°4 10°9 7§°3 74°5 76°2 the 1°6 
‘890 51°9 II'9 63°8 64°0 61°8 2°2 2°2 
"940 40°0 12°5 52°5 53°35 50°6 2°79 3°3 
"990 29°5 12°5 42°0 43°3 43°4 orl 2°5 
5876-040 20°9 11°9 32°8 34°4 34°4 oo 2°9 
*Ogo I4'l 10°9 25°0 26°5 26°2 O°3 2°0 
"140 0g"! 09°5 18°6 20°3 20°2 orl 3°0 
































Successive columns denote: (1) A, the wave-length ; (2) C,, the computed Doppler 
profile of the strong components ; (3) C2, the corresponding profile of the weak components ; 
(4) C;=C,-+-C,, the computed blend; (5) C, the computed blend, subjected to photometric 
distortion; (6) O, the observed blended profile: (7) O—C ; and (8) the standard deviation 
o of the observed values. 


4. The profiles of the Ds line 

4.1. Seven notably successful plates showing the fringe system of the D, 
line due to neutral helium at 5875-6A were selected for profile determination. 
Apart from its astrophysical interest, this line is particularly suitable as a test 
of the method, since it is not observed in absorption in the photosphere, which is 
disturbed only by water vapour lines in the wings (beyond the wave-length range 
measured) at 5876°162 and 5876-452 A, and, within the line, only by much weaker 
lines of water vapour the effect of which we may fairly safely neglect. 

The reduction was carried out by the methods of Section 3, and requires only 
a little further comment. ‘The choice of a tangential slit in making the exposures 
was taken only because a suitable device for accurate guiding, which would have 
permitted the use of a slit slightly inclined to the tangent position, and encroaching 
a little on the photosphere, was not available at the time. With the arrangement 
used, the limb does not appear on the plate, and we can do no more than to assume 
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that its mean position coincides with the base of the slit, which is what we en- 
deavoured to achieve by guiding. Consequently, our individual absolute height 
determinations are more uncertain than they need have been, but the relative 
heights and the mean heights from several plates can hardly be seriously wrong. 

The scheme for wave-length determination worked fairly well, except for 
the lowest of the heights measured, 1500km, where an apparent line shift of 
0-08 A appeared. It is believed that this effect is not real, but is due to the choice 
of the centre of the D, line to determine k,. Krat (11) has found that below 
1500 km the centre of the D, line is shifted by asymmetry towards the red. The 
use of fringes in this lower region as a standard, though not for profile measure- 
ment, accounts for our blue shift. Its effect on the profile determination is 
practically negligible. 

The scattered light was rather high, mainly due to the poor condition of the 
final Coudé mirror. This factor increases the random error, as can be seen from 
the very ragged profiles of Fig. 3(A), and also limited the range of A over which 
the microphotometer (a Moll thermocouple instrument) could usefully record. 
For the same reason, it was necessary to use a slit height corresponding to a 
wave-length interval of o-18A. A correction for this smearing is therefore 
required in interpreting the profile, but owing to the much greater width of the 
helium line, the effect is not too serious. The tracings were made over two, 
three or four fringes, at intervals of 0-126 mm (1 division on the head of the plate- 
raising screw). 

Profiles were determined for heights of 1500, 2250, 3000, 3750 and 4250km 
above the solar limb, and are displayed in Fig. 4. ‘Those for the extreme heights 
are means from only 2 plates; intermediate heights depend on 4 or 5 plates each. 
The peak intensities at the various heights have been brought into arbitrary 
agreement in the figure. Since the profiles show no systematic variation with 
height, they can be combined to derive a mean profile, corresponding to a mean 
height of about 3000km. ‘The mean points, together with their standard devia- 
tions, are shown in Fig. 5. In deriving this mean the extreme profiles have 
been given half weight. 

4.2. Interpretation of the profiles—The D, line of helium originates from 
the transition 2*P-3°D. The three components of the main triplet structure of 
the line have recently been completely resolved by Brochard, Chabbel, Chantral 
and Jacquinot (4) using a discharge cooled by liquid hydrogen. ‘The measured 
intensities agree very well with the theoretical values. The two stronger com- 
ponents are separated by only 0-04 A, and these together with three much fainter 
components (the line is really sextuplet) may be lumped together, to form with 
the third line of the triplet what is equivalent to a doublet with components at 
587561 and 5875-96 with intensity ratio 7°96:1. In our analysis we have used 
the now conventional value of 8:1. 

The fact that the profile of D, does not vary with height in the chromosphere 
within the range of heights measured makes it likely that self-absorption in the 
line may be neglected. From this it follows that the intensity ratio of 8: 1 may be 
assumed for the fine structure of the line under chromospheric conditions. It 
also follows that the broadening of the line is not a curve of growth effect, but 
represents faithfully the distribution of intensity with wave-length prevailing in 
the atomic population. The most obvious source of the broad distribution is a 
combination of thermal and turbulent velocities, and this explanation is rendered 
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Fic. 4.—Observed profiles of the Dg line at various heights. 
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Fic. 5.—Observed mean profile, compared with a computed Doppler profile. The symbols Cy, Cs, 
C and O have the meanings assigned to them in Table I. 
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that its mean position coincides with the base of the slit, which is what we en- 
deavoured to achieve by guiding. Consequently, our individual absolute height 
determinations are more uncertain than they need have been, but the relative 
heights and the mean heights from several plates can hardly be seriously wrong. 

The scheme for wave-length determination worked fairly well, except for 
the lowest of the heights measured, 1500km, where an apparent line shift of 
0:08 A appeared. It is believed that this effect is not real, but is due to the choice 
of the centre of the D, line to determine k,. Krat (11) has found that below 
1500km the centre of the D, line is shifted by asymmetry towards the red. The 
use of fringes in this lower region as a standard, though not for profile measure- 
ment, accounts for our blue shift. Its effect on the profile determination is 
practically negligible. 

The scattered light was rather high, mainly due to the poor condition of the 
final Coudé mirror. This factor increases the random error, as can be seen from 
the very ragged profiles of Fig. 3(A), and also limited the range of A over which 
the microphotometer (a Moll thermocouple instrument) could usefully record. 
For the same reason, it was necessary to use a slit height corresponding to a 
wave-length interval of o-18A. A correction for this smearing is therefore 
required in interpreting the profile, but owing to the much greater width of the 
helium line, the effect is not too serious. The tracings were made over two, 
three or four fringes, at intervals of 0-126 mm (1 division on the head of the plate- 
raising screw). 

Profiles were determined for heights of 1500, 2250, 3000, 3750 and 4250km 
above the solar limb, and are displayed in Fig. 4. ‘Those for the extreme heights 
are means from only 2 plates; intermediate heights depend on 4 or 5 plates each. 
The peak intensities at the various heights have been brought into arbitrary 
agreement in the figure. Since the profiles show no systematic variation with 
height, they can be combined to derive a mean profile, corresponding to a mean 
height of about 30cokm. The mean points, together with their standard devia- 
tions, are shown in Fig. 5. In deriving this mean the extreme profiles have 
been given half weight. 

4.2. Interpretation of the profiles—The D, line of helium originates from 
the transition 2*P-3°D. ‘The three components of the main triplet structure of 
the line have recently been completely resolved by Brochard, Chabbel, Chantral 
and Jacquinot (4) using a discharge cooled by liquid hydrogen. ‘The measured 
intensities agree very well with the theoretical values. The two stronger com- 
ponents are separated by only 0-04 A, and these together with three much fainter 
components (the line is really sextuplet) may be lumped together, to form with 
the third line of the triplet what is equivalent to a doublet with components at 
5875-61 and 5875-96 with intensity ratio 7-96:1. In our analysis we have used 
the now conventional value of 8:1. 

The fact that the profile of D, does not vary with height in the chromosphere 
within the range of heights measured makes it likely that self-absorption in the 
line may be neglected. From this it follows that the intensity ratio of 8: 1 may be 
assumed for the fine structure of the line under chromospheric conditions. It 
also follows that the broadening of the line is not a curve of growth effect, but 
represents faithfully the distribution of intensity with wave-length prevailing in 
the atomic population. The most obvious source of the broad distribution is a 
combination of thermal and turbulent velocities, and this explanation is rendered 
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| ae Yen cake Vee | | | | SD ick 
-O7 -06 -05 -04 -03 -O2 -0O! af tes o2 o3 o4 o's o6 07 


— 











Fic. 5.—Observed mean profile, compared with a computed Doppler profile. The symbols Cy, Ce, 


C and O have the meanings assigned to them in Table I. 


3 








34 P. }. Treanor Vol. 117 


immediately plausible by the approximately Gaussian form of the determined 
profile, even before correction for microphotometric distortion and for fine 
structure, since neither of these effects will greatly alter the form of the profile. 

Using these clues, we have sought for a Doppler profile of half-width slightly 
smaller than the half-width of the observed curve, which after smearing with the 
microphotometer slit pattern, and after the addition of a fine structure component 
of the same smeared form, and 1/8th of the intensity, distant 0-35 A, will reproduce 
the observed curve. For this purpose we have approximated the microphoto- 
meter pattern to a Gaussian form J = J, exp[—(6A/o"108)"]._ ‘The result is shown 
in Fig. 5, and Table I. A true profile of the form / = J, exp[ — (8A/o-340)*], after 
smearing and adding the fine structure component, reproduces the observed 
profile to within the limits of experimental error. Note that the asymmetry 
of the combined profile has moved the peak intensity towards the satellite by 
0-025 A. 

If the broadening of the line were entirely thermal, its profile for an assumed 
temperature of 4800deg. would be given by 


I =I, exp [—(8A/DA)*] (11) 
where ae 087 
amt — . (12) 


The distribution of turbulent velocities which in combination with this thermal 
distribution will reproduce the deduced profile is only slightly narrower than the 
profile itself, and corresponds to a turbulent velocity (at height 1/e in the distri- 
bution) of 16-7 km/sec. 

The profile that we have determined (DA=0-34) is slightly narrower than 
that of Peripelkin and Melnikov (14) who like us found Doppler profiles in- 
variant with height but with DA=0-377. ‘This is the only result with which 
strict comparison may be made, since Unsdld’s measures (24) yielding DA=0-316 
rested on rather slight and uncorrected observational material, and Krat’s recent 
profiles (10, 11) (DA=0-26), refer mainly to lower levels. 

4.3. A note on the D; gradient.—The absence from the plate of the solar 
limb, besides creating some uncertainty in the determination of heights, left no 
means for determining the effect of seeing on the measured gradient of the 
chromosphere. In view of both these limitations, the determination of the 
intensity gradient was not pursued seriously. A pilot investigation of our 
profiles, however, yielded an intensity distribution of the form J= J, exp(— Bh), 
with a mean value of B=0-42 x 10-*cm™!+0-06. Using ordinary spectra taken 
a few months previously with the same telescope and spectrograph, Bishop (2) 
found an identical mean value for the D, gradient, before correction for scintilla- 
tion. This correction, however, increased his value by 50 per cent. 

In conclusion, the author wishes to acknowledge the help and encouragement 
he has received from Professor Plaskett, most particularly in connection with 
the deduction of the profiles. Thanks are also due to Dr Kuhn and his colleagues 
at the Clarendon laboratory for the excellent silver films and the measurement of 
their optical constants; to Dr Bray for the calibration data of the step wedge, and 
to Mr Bishop for providing me with much useful information from his own 
researches, partly unpublished, on the D, line. Finally the author acknowledges 
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his indebtedness to Balliol College for the Research Fellowship during the tenure 
of which this research has been carried out. 
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OSCILLATOR STRENGTHS 
FROM ARC SPECTRA OF DILUTED COPPER ALLOYS 


C. W. Allen and A. S. Asaad 
(Received 1956 December 11) 


Summary 


The method of determining absolute oscillator strengths from arc spectra 
of diluted copper alloys has been described. ‘The scale has been averaged from 
various absolute measurements and calculations that are available. Absolute 
values are tabulated for certain lines in the spectra of Al, Si, Cr, Mn, Fe, Co, 
Ni, Cu, Ga, Ag, Sn, Pb, and Bi. A relation between oscillator strength and 
excitation potential has been found as in earlier work, but this is now shown to 
be due mainly to the weakness of those transitions for which the moving 
electron is equivalent in the lower level. Some empirical rules for estimating 
oscillator strengths in complex spectra are given. 





1. Introduction.—Oscillator strengths are essential to the interpretation of 
astronomical spectrum line intensities. More values are required and those 
available are strongly in need of comparative check. 

An independent measurement of the absolute oscillator strength (f-value) 
of a spectrum line involves the difficult technique of determining the number 
of atoms in the medium emitting or absorbing the spectrum. Relative measure- 
ments are much more straightforward. Inthe present work relative measurements 
only are made but these are standardized by any absolute measurements or 
calculations that are available. ‘The method allows the comparison of one atom 
with another to give absolute f-values for a large number of atoms. The essential 
requirement for this purpose is a source of light containing several atoms in 
known abundance ratios. Relative measurements then lead to (a) absolute 
f-values of atoms for which no values were previously available, and (6) checks 
on absolute values previously observed or calculated for different atoms. 

This procedure has been applied to the determination of absolute f-values 
of the neutral atoms Al, Si, Cr, Mn, Fe, Co, Ni, Cu, Ga, Ag, Sn, Pb and Bi. 
The source was a direct current arc between poles of copper alloys containing 
known and very small amounts (<o-o5 per cent) of the other elements. ‘The 
poles were obtained from Johnson, Matthey and Co., London. 

2. Experimental procedure.—The alloy arc spectra were obtained by means 
of a 10-foot Eagle-mounted concave grating spectrograph in the first order. 
The ruled grating area was 7°5 cm x 4:0 cm, the rulings 14 438 lines to the inch, 
and the first order dispersion 5:54 A/mm. A rhodium on quartz step filter was 
placed at the sagittal focus a little way in front of the slit to provide photometric 
calibration. A lens close to the step filter focused the arc onto the grating and 
the arc length was such that the arc image came just within the grating ruled 
area. The slit width was o-o6 mm—about three times the theoretical resolving 
distance. 
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‘The data for the absorptions of the filter steps were provided by the makers 
(Hilger), and checked as a function of wave-length by means of a photo-cell. 
The absorption of the six steps was rather accurately neutral. Any differential 
absorption by the quartz plate was eliminated by passing both arc and colour 
correction spectra through it. 

‘The spectrographic regions and the photographic emulsions employed were 
3200-4600 A on Ilford Zenith plates and 3g00-5400A on Ilford Pan. H.P.3 
plates. D76 developer was used at 65 deg. F. 

The arc was operated from a 240-volt d.c. generator, the arc current was 
3 amps, the pole separation was 3 mm, the poles were always similar and the 
negative electrode on top. ‘The copper alloy electrodes were 3} inches long 
(at start) and | inch diameter. ‘The alloy rods were available in two series, one 
labelled CA containing Sb, Bi, Cr, Ga, Pb, Ag, Sn, Al and the other labelled CB 
containing Co, Fe, Mn, Ni, Si. Each series was in nine concentration grades 
extending from grade o containing 0-05 per cent by weight of each alloy to 
grade 8 containing about o-ooo1 per cent of each. The compositions were given 
by the makers (Johnson, Matthey and Co.). Grade o was determined by chemical 
analysis and the other grades calculated from the amount of dilution with pure 
copper. 

‘Two colour standards were available. During the early part of the work 
use was made of a tungsten strip lamp operating at a brightness temperature 
that was measured by an optical pyrometer. This was satisfactory in the visible 
spectrum but for A<3800A the faintness of the source and the absorption of 
the glass bulb produced errors. Later, when a heavy electric generator became 
available, use was made of the positive pole of a carbon arc. For this purpose 
German ‘“ Spektra-Kohlen”’ rods of } inch diameter were inserted in an arc that 
held the positive rod horizontally and the negative at an angle of 120° to it. 
A high enough current (10 amps) was forced through the arc to bring the positive 
pole to the volatilization temperature of 3995 deg. K (1). The arc was focused 
onto a screen (at the spectrograph sagittal focus) which contained an aperture to 
select the positive pole. 

3. Measurement.—Photometric measurements were made by means of a 
Hilger Recording Microphotometer. With very many narrow lines to be 
measured it was necessary to run the photometer rather faster than would allow 
full photometric resolution, but this did not affect the photometry since the same 
speed on similar lines was used in the various filter steps for obtaining the 
calibrations. 

For each plate a calibration curve was prepared which enabled the micro- 
photometer deflection for the centre of any line to be converted into intensity. 
The continuous background was also converted into intensity and subtracted. 
The curves for various plates were on arbitrary scales but one plate was taken 
as standard and the others converted to this standard by comparing the faint 
copper lines. Comparison showed that the standardized individual lines had a 
probable error of 15 per cent. 

Alloy lines were all sharp enough for the observed width to be entirely 
instrumental and central intensities were used without correction. However 
the copper lines had a great variety of widths (because of autoionization) and 
for these the intensity was represented by the central intensity multiplied by the 
ratio of the observed to the instrumental width. 
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Curves giving a correction for sensitivity as a function of wave-length (colour 
correction) were prepared from the density measurements of the colour standard 
spectra. Independent calibration curves were prepared for this purpose using 
a continuous spectrum and step filter. ‘The densities at various wave-lengths 
were converted into intensities by the calibration curve for the plate and then 
compared with the calculated intensity. For the calculation of carbon arc positive 
pole intensity the black-body formula was used with a colour temperature of 
3950 deg. K (1). For tungsten strip lamp calculations the colour temperature 
was obtained from the optical pyrometer brightness temperature using the 
optical constants of tungsten (2). In the region 3900-5400 A the colour correction 
from the carbon arc and tungsten lamp agreed well and a mean was used. In 
the region 3200-4600A there was disagreement for A<3800A. However, 
throughout most of this region (A> 3400 A) it has been possible to compare Fe 
line intensities in the arc with values expected from King’s furnace absorption 
measurements (3). This comparison also leads to a colour correction which 
agrees well with the carbon arc but not with the tungsten lamp. ‘The colour 
correction for the 3200-4600A region has therefore been averaged from the 
carbon arc and from the Fe line intensities. Lines in the two regions are 
standardized to the same scale by means of the overlap 3900-4600 A. 

4. Conversion to wf-values.—The measurements lead to a set of line intensities 
I on a uniform but arbitrary scale. 

The emission intensity of a spectrum line from an upper level 2 to a lower 
level 1 is 

I=const. N4( —f2,)/A* = const. No(w,/a»)f,0/A*, 


where N, is the number of atoms in the upper level, w, and w, are the upper 
and lower statistical weights, f,, and f,, the emission and absorption oscillator 
strengths, and A the wave-length. We use wf (=a f,.= —@f2,) to describe 
-our results since this is symmetric in the initial and final levels. Then 


wf = const. (w,/N,)A*/. (1) 


Boltzmann’s equation relating N,, the population of atoms in the 2 level, with No, 
the population in the ground term o, may be written 


(w2/No)/(wo/ No) = exp (E,/kT) (2) 
in the usual notation. Hence 
wf = const. A3/(m,y/ No) exp (E,/kT). (3) 


This equation converts the measurements / into values of wf on a scale that 
is common for all atoms concerned. For its application we need relative values 
of Ny and also 7. 

The values of N, are in proportion to the concentration by mass of the element 
in the copper alloy divided by the atomic weight. If an appreciable fraction of 
the atoms is not in the ground term (at 4300 deg. K) we use the partition function 
in place of w. The partition functions have been obtained from Claas’ table (4) 
with a little extrapolation. 

The excitation temperature 7 was obtained by a comparison of our A*J with 
the wf values published by King and others (3, 5, 6, 7, 8, 9). For a particular 
atom (3) becomes 


log (A°7/af) = const. — E, . 5040/7. (4) 
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The slope of log (A*//wf) against E, gave the values of 7 in Table I. ‘The E, 
ranges for Mn and Ni were too small (<1 volt). ‘There were good ranges for Fe 
(2 volts), Co (3 volts) and Cr (3 volts) all of which appeared to give a satisfactory 
solution. However T from Cr is found to be anomalous; a similar anomaly 
has been found by Meggers et al. (10), by Sandage and Hill (11), and by Bretz 
(12). Cr, Mn and Ni were therefore given low weight and the temperature 
adopted was 4300deg. K. For comparison it may be noted that Meggers (10) 
found an excitation temperature of 4940 deg. K for a 10 amp copper arc. 





TABLE | 
| Atom Fe Co ce Mn Ni 
T in deg. K 4350 4250 3800 4500 4500 
Reference (3) (8) (5) (9) (6, 7) 








In the present work it has not been possible to use the method of temperature 
measurement devised by Schuttevaer et al. (13) since the spectrum lines used 
in their calculations (5153 and 5700 A) were too strong and liable to self-reversal. 
Attempts were made to determine 7 by comparison of measured arc line 
intensities with values calculated by various means. However, this led to some 
surprising anomalies which have already been described (14) and which invalidate 
any determination of 7 by this method. Attempts were also made to relate the 
excitation correction to E, by means of a curve instead of a straight line, but 
there appeared to be no evidence of curvature and the linear solution at 
T = 4300 deg. K was adopted throughout. If later work should show that the 
appropriate excitation temperature for our arc should be corrected to 7, the 
corresponding correction to our published values denoted log wf, , in Tables III 
and IV would be 

log wf — log mf, , = (5040/7, — 1-17)(E — 3-4). (5) 
In this expression 3-4 is taken as the mean E, in volts for the lines that provided 
the absolute standard. 

5. Derivation of the absolute scale—In ‘Table I] we compare log af,, 
representing our measures of wf on an arbitrary scale, with values of log wf 
from various absolute measurements and calculations. ‘The references are given 
for the measurements, and for the calculations (BD) indicates the use of the 
Bates-Damgaard tables (15). Lines for which our wf, values are not thought 
reliable (e.g. Al and Cu resonance lines) are omitted from the table although 
calculations and measures may be available. The differences log wf—log wf, 
represent the correction required to standardize our measurements. ‘These 
differences have been weighted and averaged to give a mean correction (— 1-11 
in the log) which is then applied to wf, to give wf, , representing our measurements 
on an absolute scale. 

The final results of (19) were received from Dr King after the final correction 
had been decided, but they do not demand any significant change in our constant. 
King also advises that the absolute values of (20) are being increased as a result 
of new vapour pressure determinations. 

6. Observational errors and tests.—The overall accuracy of the values wf,, 
can best be judged by the agreement of the various atoms as demonstrated 
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in Fig. 4 in Section 8. From this diagram one would judge the probable error 
of log wf,, to be about +0-2. However, some of the scatter is certainly due to 
the implicit assumption of LS-coupling, and the measurement error may be a 
little less. Although this final accuracy would thus appear to be satisfactory 
it is necessary to justify the reliability of assumptions and data that have been 
used. These are considered below. 








Taste II 
Atom | A log log log wf 
wha why. 
(BD) 
Si 3905 0°24 0°87 o'8o 
, (16) (17) (18) (19) 
Cr 4254 0°57 | —0°54 0°22 —o'50 —2°17 —0o'42 | 
4289 0°32 0°79 0°45 0°74 2°41 0°66 | 
3578 082 0°29 0°16 1°70 
3593 + 0°73 0°38 0°27 1‘80 | 
} 
(18) (19) 
Mn 4030 0°54 0°57 0°43 0°44 
4033 0°36 0°75 0°55 0°69 
4034 o'21 0°90 0'72 0°72 
(20) (21) (19) 
Fe 3719 060 O's! -1°O4 —O'4I —0'54 
3737 0°47 0°64 1°05 —0°54 
: (16) 
Ni 3461 + 0°54 0°57 1°19 
3524 0°94 —0O°!7 0-89 
3414 | +0°77 0°34 0-92 
(BD) 
Ga 4172 + 0°44 —0°67 0°64 
4032 +o°10 1‘O1 0°92 
(BD) 
Ag 3382 0°56 | —o'55 -0'26 



































(a) Atomic composition of arc plasma.—We have been obliged to assume that 
the atomic composition of the plasma is the same as that of the pole. Since the 
poles are systematically consumed the alloy atoms must be volatilized into the 
arc plasma at a rate proportional to the pole concentration. Furthermore, since 
the alloy atoms are in such small proportions one would expect their diffusion 
to be controlled mainly by the diffusion of the dominating copper atoms. If 
this were accurately true the composition of the plasma would be the same as 
the pole, and also the distribution of the alloy atoms throughout the arc would 
be the same as the copper atoms. One can make a rough check on this distribution 
by studying spectra in which the arc is focused onto the slit (or onto the sagittal 
focus). It has been found that the distribution of intensity throughout the arc 
is much the same for all lines concerned, although there is perhaps a slight 
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tendency for light atoms to be more concentrated near the negative pole than 
heavier atoms. In the circumstances the similarity of composition of pole and 
plasma seems to be a fair working assumption. We have yet to see whether 
the assumption can be justified for all atoms. 

The probable errors of the pole compositions are not stated but a study of 
inconsistencies in the tabulated values suggests that they are about 15 per cent 
for the grades that were employed. When line intensities 7 were plotted 
logarithmically against concentration C in the pole the slope dlog C/dlog J was 
found to be about 1-2 instead of the expected 1-0. No clear explanation of this 
effect was found but we are inclined to regard it as being due mainly to errors 
in the relative concentration from one grade to the next. Any such error will 
not affect our results provided the same grades have been used for all lines 
measured. For the faint lines that were measured in the more concentrated 
grades only, a correction for this effect was devised and applied. 

(b) Effective excitation temperature—The various parts of an electric arc 
are doubtless at rather different temperatures and not in thermal equilibrium. 
However if we consider Ny and N, in (2) to represent atomic populations 
throughout the whole arc we can define T for a whole arc. 7, thus defined, 
may vary with E, and give rise to a curved relation between log (A*//wf) and E, 
instead of (4). Our reason for choosing a single value of 7 was because no 
curvature could be detected in the data—not because any constancy of T or 
thermal equilibrium for the arc was assumed. Our results are not likely to be 
seriously influenced by errors in T unless E,>5 volts. For the Cu lines E, 
reaches as high as g volts and considerable errors associated with T are possible. 

Our adopted 7 depends almost entirely on the temperature standard of 
King’s furnace. As a check we have determined the excitation temperature 
of the Sun’s reversing layer by comparing Fraunhofer line intensities with the 
wf,, values. The solar temperature obtained is 4850+75deg.K. Solar 
excitation temperatures from other sources are summarized by Hunaerts (22) 
from whose tables we might take 4750deg.K as a representative value for 
molecules and 4400 deg. K for atoms. Our value is a little higher but not at 
all improbable, being the temperature for an optical depth of 0-06 (23). ‘The 
solar temperatures we obtained from individual elements were Cr 4754 deg., 
Mn 4996 deg., Fe 4917 deg. and Co 486gdeg. These are in good agreement 
and do not show much Cr anomaly. 

(c) Self-absorption.—It is rather difficult to prove conclusively that there 
has been no self-absorption but if it has only been slight it would produce no 
great differential effect from one atom to another. The chief evidence against 
self-absorption is that multiplets follow the LS intensity rules when expected 
to do so, and that our A*/ values are proportional to the absorption wf values 
measured by King and others (3, 5, 6, 8,9). Fig. 1 shows values of our A*/J 
plotted logarithmically against published values of wf for Cr, Fe, Co. The lines 
are in groups selected by E,. It is seen at once that there is no tendency for the 
slope to change with E, nor for the slope to be less than 1-0 as would be the case 
if self-absorption were effective. 

7. Results—The absolute values of wf derived in this work for the neutral 
atoms Al, Si, Cr, Mn, Fe, Co, Ni, Ga, Ag, Sn, Pb and Bi are set out in Table III. 
The first column contains the atom, the multiplet number from the Revised 
Multiplet ‘Table (24) and sometimes a letter a or B explained in Section 8. The 
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wave-lengths are also from the Multiplet Table. ‘The final absolute wf-values 
are denoted w/,, and expressed logarithmically. For use in our discussion we 
also quote an experimental value for o? defined (15) by 

cexp? = S|A(LYS(M) = tof, xN304A (LA (M) (6) 
where S is the line strength, o/(/) is the line strength relative to the multiplet 
strength (23, 25), 0%(M) which might be called the relative multiplet strength 


is available from tabulations (23, 26), and A is in A. The lower excitation 
potential E, is given in the last column. 








Log A*/ 
ie 
a 
i © 











Log wf 
Fic. 1.—Plots of log *I against published log wf for groups of Cr, Fe and Co lines. The numbers 
represent E, for each group. The straight lines drawn in are at unit slope. 


‘There is some doubt about the Al resonance lines. Al was detected in the 
alloys as an unexpected impurity and a value for its abundance was given in the 
maker’s report. However the low measured value of wf,, suggests that the 
measured concentration is not accurate. Sb was contained in some of the alloys 
but no lines were seen. 

The Cu results are quoted separately in Table IV. ‘The details for Cu are 
somewhat reduced because both the accuracy and astrophysical interest may be 
less than for other elements. ‘The three columns give (i) the lower term, (ii) the 
* whole number part of the wave-length, and (iii) logwf,,. The term data are 
from Shenstone (27). ‘The lines are not separated into multiplets although 
they are in the order used for other elements in the Multiplet Table (24). 

‘The copper resonance lines are far too strong for normal measurement but 
an attempt was made to correct for the severe self-reversal by evaluating the 
constants in a formula by Cowan and Dieke (28). The resulting measurement 
wf,, is probably too small. 

8. Discussion of results.—Our interests in the measurements are (i) to seek 
for useful regularities, and (ii) to investigate the excitation potential effect described 
earlier (14). 

All the wf, ,-values are plotted against E, (upper EP) in Fig. 2. They show 
the tendency to increase with E, that had been noticed before (14), but now that 
some high potential copper lines have been added it would appear that the increase 
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TaBLe IIT 

Atom dr log log E Atom \ log log E 
Muh. wf, Sexp" re Mult. 2 tf ,, Sexp* ™ 
Al(r) 3961°52 1°3? 0'8? oo Cr(51) 345560 0°07 2°5 
3944°01 1°5? ‘7? 0-0 3445 62 0°53 2°5 
(52) 3436°19 0°25 2°5 
51(3) 3905 °53 0°87 006 1°9 3441°44 o'51 2°5 
3447°43 016 2°5 
Cr(r)a 4254°35 0°54 0°34 00 3453°33. + 0°22 2°5 
289°72 "79 0°33 0'0 (67) 399284 0°09 2°7 
(2) 3732°03 2°44 ore) (69)A 3841°28 + 0°56 O°-71 2°7 
3730°81 2°59 °'o 3850°04 0°48 O'92 2°7 
(4)B 357869 0°29 1°16 oo 384898 +017) +0O'91 2°7 
3593°49 0°38 I'l5 oo 3854°22 O°17 +o8o0 2°7 
(10) 4496°86 1°12 °'9 3857°63 +O'13 -o-88 2-7 
454596 1°37 °'9 (70) 3819°56 0-22 2°7 
4580-06 1°71 °'9 (71) 3815 °43 0°09 a9 
{21)B 4600°75 1°21 0°48 1°'o (139) 3804-80 0°88 3°0 
4616°14 1°5 O's5I I°O 3797°72 + 0°46 3°0 
4591°39 1°77 ~0°54 Io 3793 °29 0°05 re 
4613°37 1°64 —O0°55§ 1°0 3793 °88 o'ol 3°0 
{22)B 4351°77 0°66 0°84 I°'o 3792°14 0°23 3°0 
4344°51 0°64 0°63 I°o 3791 °38 0°09 3°0 

4384°98 1°40 0°72 1°'0 
4371°28 1'2 —o'69 I°'0 Mn(2)B 4030°76 0°57 0°65 oo 
4359°63 1°22 —0°65 1°0 4033°07 0°75 o'71 oO 
4351°OS 1°38 0°67 1°00 4034°49 °'9g0 0°68 0-0 
(23)B 3919°16 0°89) = —0'95 _I°0 (5)A 4041 °36 068 0°89 2°1 
3908 °76 1°07 o'81 1°0 4055°54 0°28 +o'85 21 
3941°49 1°34 o'7I 1°'0 4070°28 o'59) «+ 0°88 2-2 
392864 1°25 0°75 1°0 4018°10 0704 +0°8Q 2°1 
3921°02 1°25 —0°69 I°0 4035°73 +028 0°99 2°I 
3916°24 92 og! I°0 4048-76 -Or13 0:87 2:2 
3886-79 1°32 e-69 I°0 4083 °63 0°00) «6—+-0°72 2°2 
3885-22 1°25 0°69 1°0 (6)A 3806°72 ~+0°'7g «-+-0°80 2°1 
3894°04 1°58 0°77 1°0 3823°51 ‘ 0°66 +o'82 2:1 
(24) 3849°53 1°71 1'o 3843 °98 O'21 o-99 2°3 
3852°22 1°98 1°o 3809°59 0°02 0°65 2°1 
(32) 4571 °68 0°54 2s 3823°89 + 0-01 0°65 2:2 
(33)A 4526°47 O'17 0°25 2°5 383386 0°08 o'82 2:2 
4535°72 0°03 + O'24 2°5 3839°78 0°07 0°83 2°2 
4544°62 0°25 +O°3I 2°5 (S)A 3577'°88 +-0°61 o'goI 2°1 
(38)a 3963 69 0°97. +0O°9OI 2°5 3607 °54 o°'15. +0°99 2°! 
3969°75 0°86 +O'QO 2°5 3608-49 0°42 Cree 23 
3976 °66 0°79 +087 2°5 (16)A 4823°52 O'rs 0°35 2°3 
3983 °91 0°69 0°92 2°5 4783°42 0°06 +0°36 2°3 
399112 0°63 «+ O95 2°5 4754°04 0°04 0°37 2°3 
3991 °67 0°30 + 0°87 2°5 (18)a 354780 o'81 O31 2°3 
(40)A 3816°17 0°13 0°03 2°5 (21)A 4762°38 0°54 0°73 2°9 
(438 3743°58 0°75 0°35 2°5 476643 0°45 o8r 2-9 
3749°00 0°34 0°05 2°5 476586 o'2. 079 2°9 
3757 °66 o'21 6+0°12 2°5 4727°48 O°14 089 2°9 
(47) 3639°80 o°41 2°s (22)A 4451°59 0°52 0°84 2°9 
3636°59 O'31 2°5 4464°68 0°28 o'89 2°9 
3649°00 —o-06 2°5 4470°14 0°07 088 2°9 
4472°79 0°28 =6+-0'87 2°9 








4436°35 
4453°00 
4502°22 
4498-90 
4490°08 


3859°91 
3886-28 
3899°71 
3824°44 
3856°37 
3878°58 
38905 66 
3922°91 
3930°30 
3927°92 
3920°26 
3719°94 
3737°13 
3748°26 
367992 
3705 °57 
3722°56 
3733°32 
3440°99 
3490°58 
3476°70 
5269°54 
5328-04 
5371°49 
5397°13 
5429°70 
3820°43 
382588 
3834°22 
3840°44 
3878-02 
3872°50 
3865 °53 
3734°87 
3749°49 
375824 
3763°79 
3767°19 
3709°25 
3727°62 
3743°36 
3798°51 
3799°55 
3795 00 
3787°88 
3812-96 
3790°10 


log log E Atom 
wf Fexp” ; Mult. 


—o'08 +0°93 2°99 ~=©Fe(23)a 
o'20) «6+0°96 2°9 

—0'13, +0°97 2°9 
0°05 0°97 2°9 


0°27 o'89 2°9 (24) 
0°63 —0'70 0°0 
1°05 0°78 ov! (41) 
1°42 o'71 Ov! 
IIS —O*52 0°0 (42)A 
1°28 -0'77 O'1 
1°36 0°79 OT! 
1°62 o81 ov1 
1°60 “0°95 O'l 
1°46 —o'95 Ov! (43)A 
1°49 0°92 O'1 
1°66 —o'85 oF 
o'51 —o'76 o-o 
0°64 -0°72 Ov! 
1°04 —0°70 O'F 
1°45 —o0'84 0-0 (45)A 
1°29 0°82 o'1 
1°22 —0°73 O'1 
1°48 —o°84 ov1 
—1°07 —O°7§ O°! (152)A 
1°24 —0°63 o'1 
1°58 —o'6o o1 
1°40 o'9 
1°44 °o'9 
1°53 1’o 
1°89 °'9 
1°74 1° 
0°33, +0°55 09 = C03) 
0°09 «+ 0°47 0°9 
+-Or12 +069 1°o 
0°29 0°53 ~«+1°0 (4)B 
—o62 O°3I 10 
0°58 +0°37 1:0 (5)B 
~0'68 0°43 10 
+0:38 +048 o9 
0°29 0°58 o-9 (6)B 
O°3I +0-76 1°0 
0°02 +060 1°o 
0°22 +0°57 1°0 (17) 
— 0°39 e'70 O90 (18) 
0°37 0°54 I°'O 
O'51 +o'58 1°o 
~O°7 0°28 09 (19) 
0°50 0°40 1'O 
—o'48 0°44 1°0 (20) 
-0°72 +0°38 1°0 
—0'68 I’o 
—F*s2 I'o 


3581°20 
3647 °84 
3585°71 
3586-98 
3570°10 
3565°38 
3558°52 
4383°55 
4404°75 
4271°76 
4307°91 
4325°76 
4202°03 
4250°79 
4045 82 
4071°74 
3969-26 
4005 °25 
4143 °87 
4132°06 
3815 "84 
3827°82 
3841°05 
3902°95 
4260°48 
4235°94 
4187°80 
4187-04 
4299°24 
4271°16 
4250°12 


3909°93 
3979°52 
4027°03 
3575 °36 
3594°87 
3513°48 
3529°03 
3533 °36 
3431°58 
3442°92 
3455°24 
3941°73 
3873°12 
3873°95 
3957°93 
3627°81 
356495 
3521°56 
3605 °37 
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Atom 
Mult. 


Co(21)a 


(22)A 


(23)A 


(25) 
(28) 
(29) 
(31) 
(32) 


(33) 
(34)A 


(35)A 


(36)A 


(58) 
(62) 


(64) 
(95) 
(98) 
(115) 
(145) 
Ni(3) 
(4) 
(5) 


_ (6) 


3502°28 
3506°31 
3523°42 
3574°97 
3560°89 
3453°51 
3529°82 
3509°84 
3495 68 
3367°11 
338522 
3405 °12 
3417°15 
3433 °04 
338816 
3474°02 
3462-80 
3412°34 
4121°32 
4118°77 
395292 
4092°39 
4110°53 
3995 °31 
4045°39 
393596 
3997°90 
3842°05 
3845 °47 
389407 
3569°37 
3587°19 
3704'06 
3489°40 
3518°34 
4086-30 
3732°39 
3730°48 
3543°26 
373593 
3708 °82 
3662-16 
3702'24 


3437°28 
3548°18 
3670°43 
3664-10 
3391 °05 
3571 °87 
3369°57 
3500°85 
3612°74 
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log 
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—0°02 
—o'16 
— 0°39 
—o-70 
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0°50 
—0°04 
—0°23 
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:°s3 
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—6°S7 
—0'28 
-0°88 
0°05 
—O°IS 
—O°I4 
0:26 
0°38 
1°17 
—0'95 
1°10 
—0°13 
—I'rl 
-0°83 
“Fo 
—o'62 
0°07 
o'17 
0°45 
+-0°27 
—0°79 
0°33 
—O'OS5 
—o-46 
— 0°34 
—o'si 
— 0°47 
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—0o'50 
—O'05 
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06 
°0°6 
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06 
06 
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O°5 
o'4 
06 
0°6 
06 
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Mult. 


Ni(7)B 
(8) 
(9) 

(17) 


(18)a 


(19)A 


(20)A 


(22 
(30) 
(31) 
(32) 
(33) 
(35)A 
(36)A 
(37)A 


Ga(1) 


Ag(1) 
Sn(2) 


(3) 
(5) 


Pb(1) 


3371°99 
3306°17 
3320°26 
3461-65 
3452°89 
352454 
349296 
3510°34 
3610°46 
3597°70 
3414°76 
3515 °05 
3458°47 
3433°56 
3361°56 
3392°99 
3446°26 
3423°71 
3472°54 
3548°18 
3315 °66 
3783°53 
3973°56 
3831 69 
3858-30 
3674°15 
3807°14 
3775°57 
3619°39 
3566°37 
3380°57 


4172°05 
4032°98 


338289 


3330°62 
380102 
3262°34 
4524°74 


3639°57 
4057°81 
3683-47 
3739°94 
3572°73 


3397 °21 
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Atom 
Mult. 


Mn(22)a 


Fe(4)B 


(5)B 


(6)B 


(15)B 


(20)A 


(21)a 





4436°35 
4453°00 
4502°22 
449890 
4490°08 


3859°91 
3886-28 
3899°71 
3824°44 
3856°37 
3878°58 
3895 66 
3922°9I1 
3930°30 
3927°92 
3920°26 
3719°94 
3737°13 
3748°26 
367992 
3705°57 
3722°56 
3733°32 
3440°99 
3490°58 
3476°70 
5269°54 
5328°04 
5371°49 
5397°13 
5429°70 
3820°43 
382588 
3834°22 
3840°44 
3878-02 
3872°50 
3865 °53 
373487 
3749°49 
3758°24 
3763°79 
3767°19 
3709°25 
3727°62 
3743 °36 
3798°51 
3799°55 
379500 
3787-88 
3812-96 


3790°10 


log 
Wha, 


—o'o8 
—-0°20 
—©@°53 
0°05 
0:27 


“46 
“49 
66 
“51 
64 
‘04 


89 


“09 
12 
29 


62 
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log 


exp” 


+093 
+ 0-96 
+ O°O7 
0°97 
+089 


-0°70 
-0°78 
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—0°52 
“OF? 
°°79 
-o'81 
0°95 
~0°95 
~O°92 
—o0'85 
—o-76 
-0°72 
—0o'70 
—o'84 
082 
—0°73 
-0°84 
—@°75 
—0°63 
—o6o 


+0'55 


E, 


2°9 
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Atom 
Mult. 


Fe(23)A 


(24) 


(41) 


(42)A 


(43)A 


(45)A 


(152)A 


Co(3) 


(4)B 


(5)B 


(6)B 


(17) 
(18) 


(19) 


(20) 


3581°20 
3647°84 
3585°71 
3586°98 
357010 
3565 °38 
3558°52 
4383°55 
4404°75 
4271°76 
4307°91 
4325 °76 
4202°03 
4250°79 
4045 82 
4071°74 
3969°26 
4005 °25 
4143 °87 
4132°06 
381534 
3827°82 
3841-05 
3902°95 
4260°48 
4235°94 
4187-80 
4187°04 
4299°24 
4271°16 
4250°12 


3909°93 
3979°52 
4027°03 
3575°36 
359487 
3513°48 
3529°03 
3533°36 
3431°58 
3442°92 
3455°24 
3941°73 
3873°12 
3873°95 
3957°93 
3627-81 
3564°95 
3521°56 
3605 °37 
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Atom 
Mult. 


Co(21)a 


(22)A 


(23)A 


(25) 
(28) 
(29) 
(31) 
(32) 


(33) 
(34)A 


(35)A 


(36)A 


(58) 
(62) 


(64) 
(95) 
(98) 
(115) 
(145) 
Ni(3) 
(4) 
(5) 


_ (6) 


3502°28 
3506°31 
3523°42 
357497 
3560°89 
3453°51 
3529°82 
3509°84 
349568 
336711 
3385 °22 
340512 
341715 
3433 °04 
3388-16 
347402 
3462-80 
3412°34 
4121°32 
4118°77 
395292 
4092°39 
4110°53 
3995°31 
4045 °39 
3935°96 
3997°90 
3842°05 
3845 °47 
3894°07 
3569°37 
3587°19 
3704°06 
3489°40 
3518°34 
4086-30 
3732°39 
3730°48 
3543°26 
373593 
3708 82 
3662-16 
3702'24 


3437°28 
3548:°18 
3670°43 
366410 
3391°05 
3571 °87 
3369°57 
3500°85 
3612°74 
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—0'02 
—o'16 
—0°39 
—o-70 
—o-72 
"50 
—0'04 
—0'23 
~0°40 
—1°93 

0°95 
0°31 
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‘ o'88 
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‘26 
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‘83 
‘73 
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log 
Sexp* 


+ 0°19 
+ O°21 
0°37 
0°35 
+044 

0°48 
- 0°06 

0°00 
—0°05 

o'o! 

0°27 
0°40 
—Oo'l4 
+O'25 
+ O25 
+1°07 
9°98 


+O°17 
+0°38 
" 0°64 
+058 
+0°83 
+063 
+04! 
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O'5 
06 
06 
06 
o'4 
O'5 
°0°6 
0°6 
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Atom 
Mult. 


Ni(7)B 
(8) 
(9) 

(17) 


(18)a 


(19)A 


(20)A 


(22) 
(30) 
(31) 
(32) 
(33) 
(35)A 
(36)a 
(37)A 


Ga(1) 


Ag(1) 
Sn(2) 
(3) 
(5) 
Pb(1) 


(2) 


Bi 


3371°99 
3366°17 
3320°26 
3461 °65 
3452°89 
3524°54 
3492°96 
3510°34 
3610°46 
3597°70 
3414°76 
3515 °O5 
3458°47 
3433°56 
3361°56 
3392°99 
3446°26 
3423°71 
3472°54 
3548°18 
3315 °66 
3783°53 
3973°56 
383169 
3858-30 
3674°15 
3807°14 
3775'57 
3619°39 
3566 °37 
3380°57 


4172°05 
4032°98 


3382 °89 


3330°62 
3801 °02 
3262 °34 
4524°74 


3639°57 
4057°81 
3683 °47 
3739°94 
3572°73 


3397°21 


log 
wf xs 


1°41 
— 1°45 
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— OS 
i 0°98 
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Tas_e IV 
Copper lines 


Low d log Low d log Low ) log Low A log 


term wf a. term Tf as term WH as term WH aa 

4’°S 3247. —1°2? 2*F° 4354 —0°33 2*D® 3488 +0739 8©602°F° 3800 +017 
3273. —1°5? 3365 + 0°95 3472 + 0°61 3527 + 0°86 

m’D 5105 —3°5? 3396 —0-'46 3613 +0°91 3797. —0°28 
3594 —4°42 3384 + 0°05 3632 +0330 «-3*P® 4525 + O1L 
3727. —4°09 3354 +061 3656 0°32 4003 +0°06 
3457 —3°56 3335 10°50 3734 10°42 3979 0°67 
3337. —2°68 3381 + 0°90 3650 —o'1? 3964 —o'18 
3530 2°65 3349 + 1°06 3602 + 1°15 3811 +018 

4p*P* 4480 0°94 3329 + 0°49 3614 + 0°53 3659 052 
4015 1°77 3375 +or°rr 3721 —o'16 3837 0°42 
3825 1°05 3319 +089 3632 0+ 0°33 3817. —0°36 
3687 +0°03 z‘D® 5292 +0°92 3546 --0°08 3671 0°48 
3566 = —1°43 5391 0°42 3610 + 0°53 3645 ~ 0°63 
3481 o°75 5250 ~~ 0°40 3498 9+ 0°35 3635 + 086 
3463 —1°71 5111 0°24 3459 = O'S 3420 91 O25 
3433 2°36 5144 +0°47 3474 —0°3! 2*D®° 4866 0°83 

aT" 4a%s o’7I 5016 +0°37. 2*F® 5200 +0°34 4767 =} 0°51 
4378 +0°67 5432 +018 5034 +0°28 4050 0°08 
4415 +0°36 5354 {0°05 5076 -++-0°64 4121 0°03 
4104 —0°44 4230 —0°25 4642 +041 4080 +} 0°40 
4177. +0°06 4242 —0'09 4776 —o°40 4010 o’4I 
4248 + 0°26 4336 —o-41 4123 + 0°04 4075 O52 
4°73 —0°40 3799 +0°02 4507 + -0°98 3771 + 0°63 
411I —0°49 3780 0°39 3860 +0-91 3699 —o'09 
3975 o's! 3507 0°63 3633 10 3946 o°25 

z*F® 4651 +1°12 3655 +-0°59 3695 0°46 3888 o731 
4586 0-92 3759 + 0°06 3665 +0°09 4052 0°25 
4539 +064 3624 0°40 3652 + 0°94 3860 = O'91 
45099 + 0°50 3745 +0O°4I 3820 —0°32 3925 +028 
4704 +0°48 3627 +089 3648 +o-12 3712 +0°40 
4674 0°55 3764 —0-'24 3517 +016 3921 +O°17 
4697 ~+-0°28 3741 +0°62 3813 —o'1o 3487 +0°52 
4842 —0°57 3621 + +0°96 3544 10°44 3436 +0°32 
4328 —o'6o - 3620 +0°37 3641 +0°23 


in log wf is not linear (as our previous analysis seemed to imply) but is bounded 
by a curve with a limiting log wf of about 1-0. This is in better agreement with 
f-sum rule expectations. 

In our previous work it was seen that the scatter could be reduced by comparing 
observed wf with the values calculated on the basis of the Bates-Damgaard tables 
(15). ‘The corresponding relation using the wf,,-values is shown in Fig. 3. 
We find the slope, the absolute values, and the interagreement of different atoms 
much the same as before. However the high level copper lines do not fit onto 
the linear curve. 

A more useful relation is revealed if we convert wf,, by means of (6) into 
exp’. ‘This can only be done for certain multiplets permitted in LS-coupling. 
We plot the resulting values against lower excitation potential E, in Fig. 4. It is 
now found that the values segregate naturally into two groups A and B. A study 
of the spectroscopic transitions shows that a-lines are for transitions in which 
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Fie. 2.—Comparison of tf,, with upper excitation potential. Each point represents a line. The 
curve represents the suggested maximum wf, and the straight line is the relation given earlier (14). 
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Fie. 3.—Comparison of calculated and experimental values of wf. Each point represents a 


multiplet. 


The straight line shows the relation found previously (14). 


the moving electron is not equivalent in either the upper or the lower level, 
while for the B-lines the moving electron is equivalent in the lower level. 

From Fig. 4 we might derive the following approximate empirical rules. 
For the a-lines o? increases with E, in agreement with o,2=1-0+2-5E,. For the 
B-lines there is not much evidence of the change of o? with E, but we might 
write provisionally, log o;,?= — 0-95 +0°33E,. For complex atoms these relations 
could be used as a general method of predicting oscillator strengths although 
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it is not at present known how widely the rules are applicable. ‘The relations 
could not have been expected from any studies of the f-sum rule nor the 
Bates-Damgaard procedure. It might be noted that o,° is in the range 1 to 20 
whereas values of 1 to 5 are usually found when the Bates-Damgaard method is 
applied to strong lines. There is therefore order-of-magnitude agreement 
between A-group observations and calculations. 
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Fic. 4.—Relation betzveen experimental o* and lower excitation potential. Each point 
represents a line. 
A and B empirical curves for A and B line groups. 
----— division between a and B line groups. 


The a-group are mostly d"s—d"p transitions allowed in LS-coupling. In 
one case of inter-parentage, Cr(69), o? is found to be normal instead of very 
weak. There are some cases of d"sp — d"ss and d"sp — d"sd transitions ; Mn (16) (18) 
and Fe(152) tend to be weaker than the a-group normal, and the Cu multiplets 
exhibit a variety of intensities (due in part to the failure of ZS-coupling). There 
are some other multiplets that are weaker than expected, Cr (33) (40), but much 
of this can be ascribed to departure from LS-coupling. 

The B-group are mostly d"s?—d"sp transitions but there are also some cases 
of d"s—d"—'sp, Cr(4)(43) and Fe(15). 

The discussion of Fig. 4 shows clearly that the excitation potential effect 
of (14) is due mainly to the existence of low d”s* configurations throughout the 
iron group of atoms. 

The results for the non-iron atoms do not appear to be very regular. af, 
values are compared with Bates-Damgaard calculations in Fig. 3. Si and Ga 
agree. For Al wf,, is low, but this is thought to be due to error in the composition. 
The Ag value of wf,, is also a little low. For Sn and Pb the wf,, values are 
much higher than expected. Fig. 3 suggests that this could possibly be 
another manifestation of the excitation potential effect. On the other hand the 
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applicability of the method to the measurement of non-iron atoms cannot be 
regarded as established. 

‘The authors wish to acknowledge that the idea of measuring absolute oscillator 
strengths by this method arose through a conversation with Dr W. F. Meggers 
in Rome, 1952. One of us (A. S. A.) has been in receipt of a research grant 
from the Egyptian Government. : 


University of London Observatory, 
Mill Hill Park, 
London, N.W.7: 
1956 December 11. 
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Summary 


In Part I we discuss the velocity distribution of a cloud of non-interacting 
particles in the gravitational field of a point mass. It is assumed that at 
great distances the velocity-distribution is Maxwellian, and that the cloud and 
the point mass are in uniform relative motion. Values are obtained for the 
density, for the mean velocity, and for the components of velocity-dispersion 
at points on the axis of symmetry near the centre of attraction. 

These results are applied in Part II to the problem of accretion by stars 
of interstellar matter. It is found that, in the physical conditions usually 
assumed for interstellar matter, collisions between particles are relatively 
unimportant, and the rate of accretion is so extremely slow as to be negligible. 





Part I1.—'THE MOTION OF A CLOUD OF PARTICLES IN A GRAVITATIONAL FIELD 

1. Introduction.—'Vhe problem of the motion of a cloud of particles in the 
gravitational field of a star arises in connection with the theory of stellar accretion. 
It was considered many years ago by Eddington (8); later it was treated very 
differently by Hoyle and Lyttleton (9, 10), and their theory of accretion has been 
followed up more recently by Bondi (2, 3), Dodd (5, 6) and McCrea (12). As 
first presented, both theories of accretion were purely dynamical in form, but 
the assumptions underlying the mechanism which leads to accretion were so 
different that the conclusions reached were diametrically opposed. 

Both theories were concerned with the behaviour of a cloud of particles in 
the gravitational field of a single star. Both suppose that the small particles in 
interstellar space whose behaviour they consider all have the same vector velocity 
at great distances from the star. ‘his is clearly an over-simplification, and the 
main purpose of this paper is to examine the two theories when we assume that 
the particles have a Maxwellian distribution of velocities at great distances. 
Eddington assumed that the particles were effectively non-interacting, and he 
calculated the number of particles, per unit of time, which, moving in the 
gravitational field of the star, are deflected so as to strike its surface. For densities 
sufficiently low it is reasonable to ignore collisions, and Eddington’s result is 
readily generalized (in Section 8) to cover the case in which the particles have 
a velocity distribution at great distances. 

In the Hoyle—Lyttleton theory, it is supposed that particles will collide with 
each other on the leeward side of the accretion axis—the line through the centre 
of attraction parallel to the initial velocity vector of the particles. The effect of 
these collisions, which are inelastic, is to annihilate the velocity components 
perpendicular to the axis, and the consequent loss of kinetic energy will prevent 
many of the particles from escaping from the star’s gravitational field. It is 
clear that if initially there is no dispersion of velocities then all particles will 
intersect the accretion axis and the density on the axis will be infinite. However, 
if there is a velocity dispersion, then the particles will not in general intersect 
the accretion axis, so that the density near the axis will be greatly reduced, and 
the likelihood of collisions is correspondingly lowered. : 
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For densities sufficiently high it is certain that collisions are the controlling 
mechanism; in that case the problem of accretion is one of gas dynamics, and 
the methods of Bondi and Dodd are applicable. But the circumstances which 
lead to so high a density have still to be established, and this is our immediate 
aim. We first determine the velocity distribution and density of the interstellar 
particles near the centre of attraction in terms of the velocity distribution and 
density at great distances, on the assumption that the effects of collisions between 
particles may be ignored. A special case of this problem, when the cloud has 
no mean motion relative to the attracting mass, has already been solved by Eakin 
and McCrea (7). 

2. The dynamics of a system of particles in the presence of a point source of 
aitraction.—We first examine the behaviour of a system of particles under the 
attraction of a point mass. It is supposed that the particles do not affect each 
other, either by collision or by gravitational forces. ‘This is clearly an over- 
simplification, but it is equivalent to assuming that the mean free path between 
such encounters or collisions is comparable with the radius of the sphere within 
which the attraction of the point mass seriously affects the behaviour of the 
particles. What this radius will be is not known, of course, until the problem is 
solved. It is also assumed that the motion of the point mass continues to be 
uniform, so that we are ignoring the gravitational attraction of the particles on 
the point mass, and also the transfer of momentum which would arise if particles 
were to strike the attracting mass. In these circumstances we can suppose that 
the centre of attraction is at rest, and at great distances the particles have a 
Maxwellian distribution of velocities, with the same mean vector velocity at 
each point. ‘The system is also assumed to be steady, as if a cloud of infinite 
extent were streaming past the point mass. 





Q P 





FIG. 1. 


The density of the particles at any point can be calculated in two distinct 
ways. In the more direct method we consider first a stream of particles, all 
with the same vector velocity at infinity, and moving under the attraction of the 
point mass. At infinity all particles have velocity V in the direction parallel 
to OO, where O is the centre of attraction. If v is the number of particles per 
unit volume at infinity, then vV will cross unit area of the plane at infinity in 
unit time. All particles cross the cone of vertex O, axis OO and semi-angle 
at some point. We first examine those which enter this cone at a distance between 
a and a+6a from the vertex. They all come from the cylindrical shell of radius R 
and thickness 6R at infinity. If the gravitational force is u/r? per unit mass, then 
the equation of a typical orbit in polar coordinates is 

pe 
77 RV 





(1 —cos @) + psind. 
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Consequently the relation between a and R is 


a 
ai p(t + cos yp) + Risin’. (1) 
a 
Further, ™ 
; = SR. 
= = at + cosy) + aa siny. (2) 
2? 


Now the total number of particles entering the element of the surface of the 

cone in unit time is 27R3RvV, and their normal velocity on crossing is RV/a. 

Thus the density at this part of the cone due to particles entering is 
27RdRvV/{(RV/a)(2z7a sin ys da)}, 


or 


a 


V 
—— cosec yb. 
0a 


Using equations (1), (2) we obtain for the density of entering particles the 
expression, 


2p -3¥\*? m . pe intl aan bn sin?® ne 
y >; + sin? — cosec — 4——— +sin?— +sin— — +sin?+ . 
- — 3) 2 |aV* 2 2\aV? 2) | 


‘There are also particles which leave the cone through the element of surface 
27a sin da. By the same method we find that the number of particles per unit 
volume which are leaving the cone through the selected element is 

ae. Cg f\- 12 * , 2 . 9 P\ 1?) 
ly { +e? en, 5 +ent® ~ sat <n +eint=) +. 
* \al? 2 2 |aV* 2 2\aV* 2) | 
Thus we obtain the particle density at (a, ) as 


2 : 12 oad 
»(— +sin?? ouut Sie +sint® . 
av? 2 2 





aV? 2 





This result, which refers to a stream of particles all with the same vector 
velocity at infinity, can be extended to cover a much more general case. Let us 
suppose that the stream of particles which moves in the gravitational field of 
the point mass has a Maxwellian distribution of velocities with the same mean 
vector velocity at infinity. This system clearly will have axial symmetry about 
the line through the centre of attraction parallel to the mean motion. We take 
cartesian coordinates, x, y, 2, with the origin, O, at the centre of attraction, and 
with the axis Oz directed towards the oncoming stream. The corresponding 


components of velocity are X, Y, Z: the velocity distribution function at infinity 
is 


(27)-pyh? exp {— }h?[X? + ¥2+4(Z+c)%]}. 


In this, py is the density at infinity, (—c) is the mean velocity, directed down 
the z-axis, and h-? is the velocity-dispersion. 
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We shall confine our attention to the determination of the density of particles 
on the axis of symmetry. ‘This is partly because we are reduced to the evaluation 
of a multiple integral, and on the axis of symmetry the evaluation is much 
simplified. It is also because on the axis the density will presumably take its 
maximum values. 

We shall therefore evaluate the density at the point A, distant a from the 
centre of force, on the axis of symmetry, and on the leeward side. If in the 


velocity distribution function we put X= —Vsin@cos¢d, Y= —Vsin@sind, 
Z= —Vcos@, we obtain the density at infinity of a stream of particles all moving 


with speed V along paths inclined at an angle @ to OA. In fact the number of 
particles per unit volume at infinity with speeds between lV’ and V+dV, and 
moving in the solid angle sin @d@d¢ is 
(27) 3*p,h? exp | — Sh2( V2 +c? —2Vccos@)}V* sin 6 dé dd dV. 
‘Thus from (3) we see that the contribution which these particles make to the 
density at A is 
(277) *?ph3 exp | — Lh?(V?2 + ce? —2Vccos8)} 


9 


@\-12 @ a See ; 
(7 + sin? ;) cosee S(t + sin? 4 V*sin@dédpdV. 


aV® aV? 2 
‘The density at 4 is therefore the integral of this expression, and as the integration 
over ¢ is trivial, we have 


P@) _ anys (f(a gine? ao sin?) 
ny =(27)'*h | 2 | a \aV sin s\ar sin : 


a) 


< exp | — Sh(V?+ c? — 2Vccos@)} sin 6 cosec 


- 


V? dédV. 


tol @ 


It is convenient to put AV =u, and to use two dimensionless parameters, defined 
by p=hce, q?=2uh*/a. In terms of these it is found that 


to 2 = | 2 _ . @\-12 
(a) =(27)'? | | (f +2 sin? ale + sin? >) 
w \@ 2 u- 2 


Po J ua0l @ 2 2 
x exp |— }(u® + p? — 2up cos @)} cos Y d0 du. 
. 2 

Finally, if we put v=sin }@, we obtain 

p(a) 

Po 
It is evident that we can obtain an expression for the density at distance a on 
the ‘‘windward”’ side by similar methods--or we can simply reverse the 
direction of the mean motion. It is sufficient to change the sign of c, that is 
of p. On the windward side the density is given by 


g? + 2u*x* 


2px fl 
= (2) | | exp{—4(u—p)*} exp { — 2upx?} (g? + wate dxudu. (4) 


7 ~ Mee a 


p(a) a ee eee ee a C+ 2urx* c 
= = (2) | P ! Pe i a(ut p)*; exp (2upx*; (2s uta dxudu. (5) 
The behaviour of these two integrals is examined in some detail in Section 4. 

3. An alternative treatment.—If we wish to find other properties of the 
stream of particles it is convenient to determine the velocity-distribution function 
at all points. Let us take cylindrical polar coordinates, w, ¢, 2, with the origin 
at the centre of attraction, and the s-axis directed towards the streaming cloud 
of particles. ‘The corresponding velocity components are I], ®, Z. From the 
rotational symmetry of the problem, the velocity-distribution will be independent 
of the angle ¢. The velocity-distribution function f(a, z, Il, ®, Z), must satisfy 
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the generalized Boltzmann equation of continuity, which in these coordinates 
takes the form, 


F) , Of 2 S) aH i® of ps of e 
ow 





a ow OO OZ 
In this 7? =m? + 2”, and y/r? is the force of attraction per unit mass. 

It was pointed out by Jeans that f is simply a function of the first integrals 
of the set of equations, 


r3 


— T 
+é 





agen’ a aa 
TO 2 Be-pof Does /r** 


The appeal to Jeans’ theorem has not been very productive in stellar dynamics, 
because it is rarely possible to find more than two of the integrals of the auxiliary 
equations. In the present case there are four integrals, all of which can be 
written down. They are: 

T, = 1? + @? + 77-27, 

1,=o®, 

I,=(Ilz— Za)? + r@?, | 

1,= Zam — 2( (1? + @?) + pa). 
The first is related to the energy, the second to the angular momentum about 
the axis of symmetry, and the third to the modulus of the angular momentum 
vector. 

The next problem is to find the function f(J,, Js, Jy, 1,) which reduces at 

great distances to 


(6) 


(277)? poh’ exp | — SA*[ 11? + ©? + (7 + c)?]}. 

‘This can be done by using equations (6) to express Il, ® and Z in terms of 
T,, Is, Is, 14, and r. The reduction is tedious, though quite straightforward. 
It introduces the angle 6, where rcos@=z, rsin@?=a, and we solve a quadratic 
equation for cos@. This introduces a square root, the sign of which is ambiguous, 
and we therefore give this root a coefficient 8, which will be either +1 or —1. 
Apart from this ambiguity, we can express 

1? + ©? + (Z +c)? =c? + 2cF(dj, 1, 15, 14,7) +1, + 2p/r. 
Finally, if F(J,, 12, 13, 14,7) tends to G(1,, 1, I, 1) as r tends to infinity, we have 
that (c?+2cG+J,) is the combination of the four integrals which tends to 
[1]? +@?+(Z+c)?] at great distances. We now replace the symbols J in this 
expression by their equivalent combination of space- and velocity-components. 
In this way, after much reduction, we determine the velocity-distribution function 
at any point to be 
{277)-*?p,h* exp [ — th? { e+, 


\ 


+ 





= 1,Z — BI,\* ur cos 6 + B1,"?Z? cos 6+ B1,12ZTI sin @ | (-) 
I, + pr + BI,1?Z cos 6+ BI,'711 sin @ 2 fies 

where /, has been written for (I1?+@? + Z?—2y/r), and 62=1. With c=o, this 
gives the expression obtained by Eakin and McCrea (7). It is obvious on 
inspection that this function gives the correct distribution of velocities at great 
distances. From the way it has been obtained, it must satisfy the equation of 
continuity—and by direct substitution it will be found to do so. 

Clearly we can obtain the density at any point, or the mean velocity, or the 
dispersions in velocities, by making use of this velocity-distribution function, 





iin 
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but it is also plain that integrals involving the expression (7) are unlikely to be 
very tractable. We therefore confine attention to the axis of symmetry, on 
which sin#@=0, and cos@ is +1 on the windward side and —1 on the leeward. 
It should be noticed that, with sin@ zero, the coefficient 8 is always multiplied 
by cos @, so that if we change the sign of 8 we are changing from one side of the 
axis of symmetry to the other. Unfortunately it is not clear which side of the 
axis corresponds to B= +1. For the present, therefore, we shall put cos@=1, 
and retain 8; so we have that on the axis of symmetry the velocity-distribution 
function is 





; 7 rk oe [,Z — BI,'*(ur-! — 2?) 
f= (27) **poh* exp | - Sh? 5 +1, +2¢ T, . - T+ BL ®Z } 


The density is given by 
p= || | fdll d@dZ. 
‘To perform the integrations it is convenient to put 
Il = V sinAsin v, ® = V sinAcos p, Z=V cosa. 
The integration over v is then trivial and we obtain the equation 
Pp =(27)-"2h3 | | exp | - bh? {e+ y2— = 
Po J. ( r 
IV cosA— BI,"*%(pr—! — V2 cos*A))\] ,.. .,,. 
+2¢-1 Bi Ne = \) V2dV sin) da. 
I, + pr + BI,"?V cosa 
‘The angle A takes all values from o to 7, and the speed V is positive and takes 
all values greater than (2/r)'®. After some changes of the variables, and the 


introduction of the dimensionless parameters p and g used in the previous 
section, we obtain the expression for the density 





—— 3 8.1 
Sas (2) ‘ [ exp | — 3(u—p)?} | exp ( — 28pux*) foe dx udu. (8) 


he a g? + uex2)e 
On comparison with the formulae (4) and (5), it is seen that the sign B=1 
corresponds to the leeward side, B=—1 to the windward side. This 


determines the sign of 8 in the formula (7), the expression for the velocity 
distribution. 

4. The determination of the density.—Vhe integral for the density can be 
handled simply in the case when p=o, so that there is no mean motion. This 
is the case treated by Eakin and McCrea (7). We have 


2 -@ “ 2 2 2 
me (2) | exp (— }u*) | Baboon dx udu, 


Po v / u=0 pmo (G2 + u2x® 2 
= (2) [exp (tet) tage + wey 
id | u=0 


(2) “| exp(— but) (gt +02)! udu, 


(2)'* 7 expt-ve ante 


For large values of q, the integral is expanded in asymptotic series, and we obtain 
the formula, 


p 2 1/2 . a 
P - (2) lgt+q°—q*+3¢°—-3-5¢" +...) 
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For small values of g, the integral can be expressed in ascending powers of q, 


in the form 
p 2 1/2 G)- 18 3 5 K 
mi | - exp (39°) — 9°/3 — , ef 
Pa (2) > p(3q)— 9/3 -9°/3-5 


When p is not zero, so that there is a streaming motion past the centre of 
attraction, we are primarily interested in the behaviour of p for large values of , 
and small values of p. (It will be remembered that q? is inversely proportional 
to the distance from the centre of force.) Now the exponential terms in (8) are 
negligibly small except near u=p, x=0, in the case B= +1, or nearu=0, x=27'? 
in the case B= —1. Thus, in the ranges of u and x where the integrand is not 
negligible, u2x? is small compared with g*?. We shall therefore expand the 
integrand in descending powers of g, obtaining 


p 2\12 “2 ne 
a (?) q| exp[—i(u—fp)*]| —_ exp(—28pux*) 





Po 7 oF ace | 
P . YY _, (2n—2)!(an+ - i i 
x sis +2 = 22"-In'(n—1)! ge” dx udu. 
We write 
p 2 12 a ' (2n om 2)! ! 
—_—_ = =» A oth N\ — J n+] Be 2 ASE fn A an 
Po (;) . [ ot ei ee tn'(n—1)! "4 (9) 
where 


-x +] 


A,(p)=(2n+1)| — exp[— 3(u—Bp)?] | exp (— 2Bpux*) x°" dx u?"*! du. 


~ u=O »- r=0 


Of the several ways in which the functions 4,,(p) can be found, the simplest 
seems to be to make use of the differential equation satisfied by 4,(p), namely, 
7A, 








4p" (an 1+ pty Fe +(4n?—1)A, 


dp* 
=(4n? —1) | exp[—}(u+Bp)*] u?""'du. — (10) 


The initial conditions are 


Pe) a (=> 
. "n! \an+3/’ 


when p=o. This equation has been solved numerically for the special case n =o. 
For such other values of n as may seem necessary, say n= 1, 2, 3, 4 and 5, it is 
best to make use of the three relations, 


Any 144A, +{ . an+1\ , _ anti 
x: 





B,,, a on B,, +29 








an+3 4p dp " \4° 8p? | A, 8p 4p 
1 dA,,,_ _—‘1 dA, antl , + 4 
(an+3)° dp 3p? dp rop* 7" * ~gpr BBs, 





2n+iI I I 
— pes Ae 
> (11) 











1 @A,,, -{4s re (2n + 1)(2n+ 5) 
my ele ta et 
_ P_ {(6n+ sant) , (2n+1Xan+ 3X(an+ 5)} 
2n +3 | 8p 16p* Be 


_ (2n+1)(2an+5) ant+3, 1 )\ 











2n+3 | 3ap* = 8p? Boy +1 
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In these equations the quantities B,, which occur are functions of p, defined by 
the integral 


m 


*o¢ 


: B.,(P) = | exp [ “x s(u T Bp)?] a” du. 


It will be seen that 


B,(p)= Me —Bp (" exp (— bx?) dx, 


- 0 

es: m7 
B\(p)=p | exp(—1x?)dx + exp(— |p?) pf. 

- oO & 
so that the values of these functions for different values of p can be readily 
obtained from tables of the error function. For the subsequent functions we 
make use of the recurrence relations 

B,, =(m—1)B,,.— BPB,,,,, m>1. 

‘Thus the procedure is to find values of the first few B’s for ditferent values of p. 
From these, and from the values of 4,(p) and dA,/dp obtained by numerical 
integration, we can use relations (11) to get values of A, and its first two 
derivatives. For each value of p these results can be checked from equation (10), 
the right-hand side of which is (4n?—1)B,,,.,.. These values of A, and dA,/dp 
are now used to find A, and its first two derivatives. ‘The procedure can be 
extended without difficulty, and it has been carried out as far as A;. 

It will be realized that the process has in effect to be carried out separately 
for the two values of 8. The pairs of functions A,(8= +1) and 4,(B= —1) 
behave quite differently. For 8= + 1, that is, on the leeward side, the functions 
B, tend to zero rapidly as p increases. Indeed for large p, 

B,,~m!p-""exp(— $p?). 

It is also found that as p increases, Ay tends to a constant (which we shall 
show subsequently is 7/2) and so we have from (11) that 4,,,,/A,,>(2"+ 3)/4 
as p increases. Indeed it is found that, whereas A,(0)=2"n!, the limit of 4,, 
as p tends to infinity is (2v+1)!7/n!2*"*!, which is smaller forn>1. Hence 
the series (g), which converges asymptotically quite rapidly for g> 3, and p zero, 
will converge even more rapidly for g>3, and p large. We therefore obtain the 
functions 4,(p) for large p by solving (10) in descending powers of p. 

For B= —1, the problem cannot be tackled in this way. The integrals B, 
increase rapidly with p, and we can show that 


Lt {B,,/p"}=(27)*?. 
p= 
The function A,(p) decreases as p increases—in fact for large p, 
Ay~3 a ', 


but the other A, increase with p. It is found in fact that 


Lt {A,/p*3}=1 fens 1). 


p > 
On comparing this result with equation (g), we may surmise that for large values 
of p and q, and for 8 = — 1, the approximate value of p/p, is }(q* + 2p*)p-(q* + p?)-**. 
This will be confirmed by other methods. 
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5. The determination of the density when the velocity of streaming ts high.—'Vhe 
method described for determining the density for small values of p would be 
difficult to justify if the value of p were large. An alternative method has been 
found, but it appears necessary to treat the windward and leeward sides of the 
axis separately. On the leeward side we have 





a i | ; exp[— }(u—p)*] udu | exp( ~ 2pux*) EAL ale dx. : 
Po \a] Jo “Pls | gy (g2 + u2x2)t2 | 
2\ U2 ; ~ ¥ q? +23 2 . | 

Ks (2) | . exp [—3(u—~p)] | De — 2py"/u) q?ayye® du. 


In the inner integral we can replace the upper limit by infinity* thus adding 
to the right-hand side a quantity not exceeding (27) '*gp *exp(—}p*). This 
is clearly negligible for g small and p large, and even for q and p large, but of 
the same order. ‘The inner integral now becomes 

9 24 208 

2 q° + 2y 
ex — 2pVv"u l —___—. dy, 
| p(— 2p: get yay . 


~ y=0 


I 


ig? | exp (— 2pq?u"' sinh? }t) cosh t dt, 
(=0 


*¢ 


sqexp(pq?u-') | exp(—pq?u-' cosh t) cosh ¢ dt, 


= 2g exp (pqru")K,(pqru). 
Thus, for large values of p, the density is given by 
V2 2 
P= (2)'" |" expl—1(u—p)] Lat exp (pata!) Ky pau) de 
0 “0 
If we now put u=p(1 +3), the expression can be written as 


5, (2n) pg? | exp(—pis*)exp[g*/(1 + s))Kitg?i(1 + 3)} dz. 


This integral can be expanded in inverse powers of p, using Watson’s Lemma. 
The leading term is independent of p, so that for large values of p, 

eo 8) (4! 

= TAS). 

0 
and if g is also large, we obtain 


Ae 5 ag ie: Bae 
é (=) | 4+ 34 = t+... 


‘This result confirms that the method described in the previous section can be 
applied to find the density on the leeward side even when p is large. We can 
therefore find (and use) the asymptotic form of A,(p) for large p, determined 
by equation (10). 

With B= —1, the value corresponding to the windward axis, the integral 
must be treated rather differently. We write g=scosa, p=ssina, and for 
©<a< 7/2, and for large values of s, we examine the integral 

p 2\l2 5% rl 

—= (=) | exp[—}(u+ssinx)?] | exp (2sux? sin «) 
Po 7 / u=0 / 2=0 
s* cos® a + 2u*x* 


* (costa + u2x?)l2 en 








* The authors are indebted to Professor M. J. Lighthill, F.R.S., for suggesting the method 
which follows. 
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After some changes of variable, this can be written in the form 


p a\is _ he : “*'2 ; 
—= (=) ssin?%| exp(— }s*3* sin® x) (1 + s) | exp [| — 2s? sin® «(1 + 3) sin? 6] 
B pene “0 


> 





cos? «+ 2sin?a(1+2)®cos?@ 
—s ~y 3 aps Sin 6 db dz. 
{cos? x + sin? «(1 + 2)? cos? 6} 
‘This also can be evaluated for large values of s, again using Watson’s Lemma, 
and, after considerable algebra, we obtain the result 


p i+sin?« costa(1+6sin?«)  costa(7+gsin*«+210sin® a) 


= 4. 





Po 2 sin % 8s? sin? x 32s? sin’ « 

It is now convenient to tabulate p/p) as a function of p and g. It should be 
remembered that p= /c and is thus the ratio of the mean speed to the root-mean- 
square speed of the particles at a great distance. ‘The quantity q is \/(2uh?/r), 
where y/r® is the attraction per unit mass at a distance r. The first column of 
the table gives values of r in astronomical units. ‘They have been calculated 
for the case in which the gravitational field is assumed to be the same as that of 
the Sun, and the root-mean-square velocity at infinity (A~') is assumed to be 
1 km/sec. (This corresponds to a temperature of about 120deg.) With these 
values of A and y the value of g which corresponds to r=1 astronomical unit is 


TaBLe [ 
The variation of density (p/ po) as a function of distance from the centre of attraction and of the mear 
speed at infinity 























IN p 
. oo 0°5 1'O 2°0 30 4°0 7‘° 10°O ao 
astr. units) q i 
100 4°21 3°54 2°87 2°26 1°47 - _ — oe 1°00 
50 5°96 4°88 3°91 3°02 1°87 - — _ 1°00 
25 8-43 6°82 5°42 4°14 2°46 1°72 1°40 re 1°O4 1°00 
‘= 10 13°3 10°7 8:46 642 3°70 248 1:92. «1°31 I‘l4 1°00 
= 5 18'S 15‘I 11°9 g*02 5°14 3°38 2°56 1°62 1°30 1°00 
= 29'°8 23'8 18:8 14°2 8-04 5°32 3°90 2°31 93 1°00 
2 I 21 33°6 26°5 20°! 113 7°34 5°44 3°15 2°26 1°00 
EZ Os 59°6 47°5) -37°5 283 16-0 10°3 7% £37 ..3°33 1°0o 
z o-2 94°2 7§°2 59°3 44°38 25°3 16°3 12°0 6°32 4°8c 1°00 
o'l 133°2 106°3 83:9 63°3 35°7 23°0 17°0 9°61 6°73 1°0°0 
0°05 188-4 150°3 118°6 89-6 50°5 32°6 24°0 13°6 9°48 1‘0o0 
o'o! 421°3 490° 2063 200-3 1129 72°38 53°7 30°3 2e*1 1°00 
O°or | 421°3 330°2 401°7) 452°5)  504°7) 519°3)  523°5 §26°7 527°4 = 5258°0 
0°05 188°4 rse°3 179% 208°4 225°7 .233°3 236°% -296°5 235°5 290° 
orl 133°2 106°3 12770 «14371 «#3915976 164°2 165-5 166°5 166°7  166°9 
2 o'2 94°2 75°2 SoS rors rIz8 1rtG71 rm7°O 1979 3179 410"! 
4 O0'5 59°6 47°5 568 64°0 71°4 734 74° 745 746 74°7 
; 21 336 402 45°3 50°95 519 524 §2°7 $277 528 
. 3 29°58 2338 44 oe UT UP. UO. UCSD. 3S 
» ¢ 18:8 15°1 18:0 20°3 226 2372 23% 236 236 236 
= 10 13°3 10°7 12°7 14°3 16°0 16°5 16°5 16°7 16°7 16°7 
25 | 8°43 6:82 S11 Q°'I2 10°2 10°74 10°5 10°6 10°6 10°6 
50 1 ot! 23 ote Op Fee OS TS: FE. 2a eS 
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42:1. In the third and subsequent columns are given the values of p/pp at different 
points on the axis for various values of p. ‘The upper half of the table corresponds 
to the windward side of the axis, and the lower half to that part which lies beyond 
the centre of attraction. 

There are several features of this table which call for comment. It is seen 
that on the windward side the density increases extremely slowly as we approach 
the centre of attraction unless the mean speed of the stream is very small. On 
the leeward side the increase of density is fairly rapid and becomes rather more 
rapid as the speed increases. However the density distribution is altered very 
little for values of p greater than 2. For a given q, p/p) decreases monotonically 
as p changes from — « to + ©. 

It is not clear from the table whether the condensation which is particularly 
noticeable on the leeward side extends far beyond the range of distances included. 
The methods used do not lend themselves to the determination of density for 
small values of g except in the three cases of p=0, and p-> oo on either axis. 
Some results in these cases are given below. It will be seen from this table that 
the density at distances beyond 1000 astronomical units is very little in excess otf 
the density at infinity. 








Tas_e II 
The variation of density at large distances 
r (astr. units) 500 1000 2000 
q 1°884 1°332 0'942 
p 00 (windward axis) 1°00 1°00 1°00 
p=o 1°86 I*4I 1°29 
p <0 (leeward axis) 2°59 1°98 1°55 




















‘The assumptions which underlie these calculations may be repeated. It is 
assumed that a uniform stream of particles has taken up a steady state, solely 
under the attraction of a point mass which is at rest. In these hypothetical 
conditions we ignore for the moment such things as radiation pressure and the 
possible effects of collisions. The distances along the axis to which the tables 
refer are determined by the values of g, and also by the values assumed for the 
mass of the attracting particle and the velocity dispersion of the stream at infinity. 
A change in either of these quantities would lead to a different relation between 
r and q. It is for this reason that it was thought necessary to determine the 
density for a very large range of values of g. 

6. The velocity dispersion.—Since we shall need to consider the possibility 
of close encounters between particles, it is necessary to get a measure of the 
velocity dispersion, and how it increases near the centre of force. Again we 
confine attention to the axis of the system, in order to be able to evaluate the 
integrals without too much difficulty. ‘To obtain pZ, pIl? or pZ?, we multiply 
the velocity-distribution function (7) by Z, I? or Z?, and integrate over all values 
of the velocity components. The integrals are evaluated by the methods used 
in Section 4 to determine the density. Attention is confined to small values of 
the parameter p, and large values of g, that is, to the case when the mean motion 
is relatively small, and to the region near the centre of force. 
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Taste III 


The variation of the mean speed (hZ) along the axis as a function of distance from the centre of 
attraction and of the mean speed at infinity 

















a p 
r 
P 2 I ° 
(astr. units) q Pie 
100 4°21 2‘980 1°268 ° 
10 13°3 -3°167 1°296 ° 
I 42°! 3°221 1*302 ° 
I 42°! 1°930 —0'974 ° 
10 ke 1°930 —0'974 ° 
100 | 4°21 1°935 0°977 fo) 











Some values of the mean velocity and of the velocity dispersions parallel to 
the axis (o,) and perpendicular to it (o,,) are given in Tables III and IV. The 
qualitative behaviour of these quantities could have been predicted. On the 
windward axis there is a streaming motion towards the centre of attraction; the 
mean velocity is greater than at infinity. On the leeward side the streaming 
is perpendicular to the axis, and the mean velocity is reduced. ‘The tables also 
suggest that for values of r exceeding, say, 200 astronomical units, the velocity 
distribution will be effectively Maxwellian. 


TaBLe IV 


The variation of the velocity dispersions at points on the axis 
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r (astr. units) o=2 p=! p=o 
| 

hor lioz hex hoz hog hoz 

100 1°62 3°23 2°11 =t7 2°57 

10 4°13 52°97 6°11 10°18 7°74 

I 12°73 38°02 19°17 34°27 24°34 

I 28°81 10°82 27°38 16°68 24°34 

10 9°15 3°55 x 8-69 5°35 7°74 

100 3°00 1°46 2°86 Ig! 2°57 

| J 





Part IJ.— THE ACCRETION OF INTERSTELLAR MATTER 


7. The application to particles in interstellar space near an attracting star.—We 
now apply the results of the preceding sections to different types of particle present 
in interstellar space. Of the two main points to consider, the first is whether the 
assumption required in the earlier sections—of non-interaction between particles 
—can be justified. If it can, we are able to estimate the density and velocity 
dispersion near the star. The remaining question to discuss is whether conditions 
near the star lead to so high a collision rate that accretion is necessarily a 
phenomenon in gas dynamics. In fact, it appears that, in the physical conditions 
normally assumed for interstellar matter, our hypothesis that the particles do 
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not interact has provided us with a very good approximation for the density 
distribution of matter in which collisions occur. Further, even in the region 
near the star where the approximation is not so convincing, it seems unlikely 
that collisions are so frequent as to necessitate the use of gas dynamics. 

Data about the interstellar particles have been conveniently summarized by 
Spitzer (13). We shall need to distinguish between the H1 regions, in which 
the hydrogen atoms are almost all non-ionized, and the Ht! regions, where they 
are ionized. We shall suppose, at any rate in the first instance, that the attracting 
star has the same mass as the Sun. In the H1 regions the density of hydrogen 
atoms is estimated at 10 per cm*; this corresponds to a mean free path of 20 ast. 
units. The temperature, originally given by Spitzer as 50deg.K, is now 
considered to be about 125 deg. K (15). We assume a value of 120 deg. K, which 
gives the value of h-' as 1 km/sec. With a gravitational field due to one solar 
mass, and with this temperature, the parameter q of the previous sections is 
related to the distance from the attracting star as in the first two columns of 
Table I. From the earlier calculations we conclude that at distances exceeding 
200 ast. units the velocity-distribution of the hydrogen atoms is effectively 
Maxwellian. Although individual atoms will be involved in collisions, these 
will serve to maintain the Maxwellian distribution which is predicted by the 
present treatment, in which collisions have been ignored. Within 200 ast. units, 
the theoretical distribution departs further from the Maxwellian form, whereas 
in practice collisions become rather more frequent, for the mean free path within 
200 ast. units is rather less than 10 ast. units. This implies that the method .we 
have used is not strictly applicable at this short range. ‘This can be verified 
much more directly. If we assume that the mean motion is zero at infinity, we 
can calculate (see the Appendix) the average number of collisions per unit volume 
per unit time, assuming that the velocity distribution is given by (7), with the 
constant ¢ put equal to zero. We can also find the average number of collisions 
per unit volume per unit time in which one of the two particles loses sufficient 
energy to be temporarily “captured”, i.e. so that its speed W becomes less than 
(2pu/r)'?. Itis only within about 200 ast. units that a large proportion of collisions 
lead to captures and therefore set up a velocity distribution seriously different 
from (7). The proportion of collisions which lead to captures is independent of 
the density, but the average number of collisions per unit volume per unit time 
is proportional to the square of the density. Consequently the expression (7) 
is reasonably accurate for r>200 ast. units even at much higher densities than 
we have considered. On the other hand for lower densities the number of 
collisions per unit volume per unit time falls off rapidly, so that although a high 
proportion of collisions will lead to temporary captures, the actual velocity 
distribution will not be affected seriously even at distances less than 200 ast. units. 

Although the method we have used is not strictly applicable for un-ionized 
hydrogen atoms within 200 ast. units of the attracting star, we can make some 
useful deductions about the hydrogen atoms within this distance. In the first 
place it is clear that collisions are almost perfectly elastic. This contradicts the 
Hoyle—Lyttleton theory, in which collisions on the accretion axis are supposed 
to cause an actual loss of kinetic energy. In fact, if the collisions were completely 
inelastic, then the kinetic energy destroyed in this way—and presumably 
reappearing as radiation—would be sufficient to ensure that four cubic parsecs 
of un-ionized interstellar hydrogen would radiate more energy per second than 
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does the Sun. ‘his seems to be far in excess of the energy output in the 21 cm 
band, or of any other interstellar radiation. 

It is also apparent that although within 200 ast. units most collisions lead to 
captures these will usually be of a temporary nature. Many of the atoms which 
are so captured have kinetic encrgy little below that required to escape (see 
Appendix). ‘They are likely to be involved in further collisions, most likely 
with atoms of much higher energy, and may well therefore recover much of the 
energy they originally lost. There is the further point that the average number 
of collisions experienced by any one particle is relatively small. If we confine 
our attention to the case of zero mean motion we can make a rough estimate in 
the following way. The mean free path is inversely proportional to the density 
which, at small distances from the centre of attraction, varies as r}*. ‘Thus 
A(r), the mean free path at a distance r, is A(a)(r/a)'*. Let us consider the average 
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number of mean free paths per chord of a sphere of radius a. The number of 
mean free paths along the chord AB is given by 








~(a*@—y?)'? : V2 -(a*—y*))* 
2 ll aa | ‘io 
Jeo Ar) A(a) J 2-0 
The average number of free paths along any chord is then found to be 
2q'2 "a 2v “(a*— y*)'* , ) 8 a 
2y | r dx} dy, =-~-. 
A(a) J y-9 2 \J 2-0 ; 5 A(a) 








Hence the average number of mean free paths per chord of a sphere of radius 
200 ast. units is about 30. ‘This overlooks the fact that the actual orbits of 
particles between collisions are curved towards the centre of attraction, and that 
the direction of motion of a particle is changed at each collision. Of course, if 
the stream of particles has a mean motion, then this figure would be considerably 
reduced. ‘The frequency of collisions depends on the velocity dispersion, but 
not on the mean motion. If the latter is large, the average distance travelled 
between two collisions would be very much larger than in the case of zero mean 
motion. For these higher speeds it would be better to estimate the number of 
collisions by dividing the average time taken to cross the region of relatively high 
density by the average time between two collisions. ‘This gives a much smaller 
figure for the average number of collisions in a sphere of radius 200 ast. units. 

The grains of interstellar dust in the H1 regions present no serious difficulty. 
Were gravitation the only force acting on them, they would certainly need to be 
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taken into consideration, not only on their own account, but more particularly 
because of the possible effects of collisions between hydrogen atoms and dust 
particles. However, it is clear that the force of radiation from the star on these 
grains will in general greatly exceed the gravitational attraction, except in regions 
where the density of grains is so high that a certain amount of screening occurs. 

So far, we have only considered un-ionized hydrogen atoms and dust grains. 
However it has been shown by Strémgren (14) that the brighter stars at least 
are surrounded by a large sphere in which all the hydrogen is ionized. For 
O and B type stars the radius of this sphere is of the order of parsecs, rather 
than ast. units. It now seems clear that for stars with an ionized envelope 
exceeding about 200 ast. units in radius, the effect of collisions in the un-ionized 
material outside can be neglected. Near the boundary of the ionized sphere 
there will be interaction between the neutral hydrogen from outside and the 
protons and radiation inside. The neutral hydrogen cannot penetrate far betore 
it is ionized; the protons cannot escape far before recombination occurs. ‘The 
boundary between the outer Hi and inner HI regions is quite sharp. ‘The dust 
grains do not come into consideration here, for in a region where the radiation 
is sufficiently energetic to ionize all the hydrogen, the radiation pressure will 
certainly repel the dust particles. 

Within the Hu region the density is about 1 proton per cm*. As a result 
of the strong radiation the temperature is of the order of 10‘ deg. K, so that the 
root-mean-square velocity is about g km/sec. Unless the Hil region is very 
small indeed, the ionization of incoming atoms will re-establish a Maxwellian 
velocity-distribution. Outside the region there is a mean streaming motion of 
the interstellar matter; this may persist to some extent inside, but it is likely 
to be checked by radiation pressure, and it seems more probable that the mean 
velocity within the sphere will be reduced. The sphere itself may be somewhat 
deformed by the external streaming. Even if streaming persists within the 
H 11 region, the parameter p, by which it is measured, will be small, for p= sc, 
and even if c, the mean velocity, is comparable with its value in the H1 regions, 
the constant / is only one-ninth of its value outside. ‘This reduction in the value 
of h enables us to deduce the velocity distribution in the HI regions at any rate 
approximately. Whereas in an H1 region the parameter g takes the value 1-88 
at 500 ast. units, in an H11 region the value at this distance is only 0-21. If 
therefore the H 11 region has a radius of 500 ast. units—and this is a very modest 
figure—the Maxwellian distribution which is set up there differs insignificantly 
from our theoretical distribution, (7), which is only truly Maxwellian at infinity. 
We can therefore deduce that the more serious departures from the Maxwellian 
distribution which occur at g=3 correspond to a distance of 2:5 ast. units— 
which is less than the radius of many O and B type stars. 

In the H11 regions, actual collisions between protons are extremely rare, 
but there will be close encounters between them, in which the electrostatic 
repulsions will produce deflections. The mean free path for the protons can 
be calculated by the method Chandrasekhar (4) has used in finding the mean 
free path for star systems, in which a star is deflected by the gravitational attraction 


of a near neighbour. The method leads to a formula for A, the mean free path, 
namely 


A= Ko* m*/pe! log {4mDo*|ne2}, 


where e and m are the charge and mass of the proton, v? is the mean square 
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velocity, p is the number of protons per cm’, D is the mean distance between 
protons, and K is constant. In the outer part of the H11 regions, the value of A 
is about 0-5 ast. units. However, it should be noted that, whereas for uncharged 
particles the mean free path between collisions varies inversely as the density, 
and is independent of the speed, for protons the mean free path between encounters 
varies inversely as the density, and more or less directly as the fourth power of 
the root-mean-square velocity. In these circumstances, the mean free path 
actually increases as the protons approach the star. At 2-5 ast. units (q¢= 3), the 
mean free path has already risen to 2°5 ast. units. For smaller values of r, we 
have that por”, and v=ar!. If we use these approximations to find the 
average number of mean free paths per chord of a sphere of radius 2-5 ast. units, 
the value comes out to be o-g. At r=1 ast. unit, the average number of mean 
free paths per chord is 0:12, so that only one proton in eight which enter this 
sphere will experience an encounter before leaving. It therefore seems to be 
established that encounters in the H11 regions will not be frequent, and a theory 
which ignores encounters completely will not be seriously in error. 

The calculations of this section have been based on the assumption that the 
attracting star has the same mass as the Sun. Now the parameter gq is proportional 
to (m/r)'?, so that doubling the mass of the attracting star will produce the same 
density and velocity-distribution at twice the distance. For the bright stars 
which are surrounded by H1I regions, an increase in mass would lead to an 
increase in the mean free path at a given distance. We may therefore conclude 
that collisions scarcely affect the motion at all in the Ht regions, whatever the 
mass. 

For those stars whose radiation is insufficient to ionize the hydrogen round 
them, collisions may be more effective for the more massive stars. As we have 
seen, for hydrogen atoms, the average number of mean free paths per chord of a 
sphere of radius 200 ast. units round a star with the Sun’s mass is about 30, and 
the number varies as the square root of the mass. ‘This, however, rather under- 
estimates the effect of increasing mass, for the radius of the sphere in which a 
collision is likely to lead to a temporary capture also varies as the square root 
of the attracting mass. ‘Thus the average number of mean free paths per chord 
of the sphere within which collisions are likely to lead to temporary captures 
varies directly as the mass of the attracting star. 

It therefore seems possible that some hydrogen atoms will lose sufficient 
energy by collision to be captured by the gravitational field of the star, and that 
the rate of capture will be proportional to the mass of the star. The method we 
have used in this investigation does not reveal the rate of this accretion, but 
there are two reasons for supposing that it will be small. In the first place, 
since the kinetic energy is only exchanged at a collision, the average kinetic 
energy in any region will certainly exceed the energy of escape. ‘Thus in the 
region within 200 ast. units, where, as we have seen, there will be some atoms 
which have been temporarily captured, further collisions between these and 
the fast-moving atoms will often raise the kinetic energy of the captured atoms 
above the escape level again. It is extremely important to distinguish between 
capture and accretion. A particle only undergoes accretion when it strikes the 
attracting star. Its kinetic energy may or may not exceed the escape energy. 
A particle is captured when, as a result of collisions, its kinetic energy falls below 
the escape energy. It is not automatically doomed to fall into the star, but at 
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least until involved in another collision it cannot escape to great distances. Since 
the paths of captured particles are entirely within regions of rather higher density, 
the frequency of subsequent collisions will be higher, but such collisions may 
raise the kinetic energy again. It is, of course, true that loss of energy is likely 
to make a particle pass closer to a star, where it is more likely to collide with other 
low-energy particles. And each such inward move increases the probability that 
a further collision will set the particle off on an orbit which meets the star. 
(In this connection it seems reasonable that the rate of accretion should be a 
function of the star’s radius, and should increase with radius.) If we knew more 
precisely how the velocity-distribution is modified by collisions, we could estimate 
the rate at which atoms become more or less permanently captured. Presumably 
very near the star the atoms approach a Maxwellian velocity-distribution, provided 
the original density at infinity is sufficiently large, but the region which seems to 
be vital to the problem is the one in which collisions affect the velocity-distribution, 
but are not sufficiently frequent to produce a Maxwellian distribution. 

The second feature which we have left aside, but which may well prevent 
any rapid accretion is the existence of the dust grains. ‘These, we have explained, 
will actually be repelled by radiation pressure. They will scarcely be affected by 
collisions with the hydrogen atoms which have a much smaller mass, but the 
atoms may well be affected by such collis'ons. The grains in fact serve to check 
the inflow of hydrogen, and in clouds where the density of dust grains is high, 
the grains may prevent most of the hydrogen from approaching the star. 
Circumstances which favour this situation are either weak radiation or large 
mass, so that a barrier of dust grains forms relatively close to the star. 

We can therefore conclude that in Hi regions generally, and in the H1 
regions surrounding stars of relatively small mass, collisions do not appreciably 
affect the motion of the interstellar particles. In the H1 regions round massive 
stars, collisions are likely to affect the motion, but they will not necessarily lead 
to a large scale accretion of hydrogen by the star. 

8. The modification of Eddington’s treatment to allow for a velocity-distribution 
in the cloud.—Since we have shown that collisions will not in general seriously 
affect the motion of particles in an interstellar cloud, we can calculate the rate 
of accretion by generalizing Eddington’s result, so as to allow for a Maxwellian 
distribution of velocities at infinity. (The rate of accretion for non-interacting 
particles has in fact already been given by Eakin and McCrea (7), but the method 
which follows is more direct.) Let N(V,0,¢) be the number of particles per 
unit volume in the cloud approaching the star from infinity with velocity V in 
the direction (8, ¢). The rate of accretion for these particles is, by Eddington’s 
formula, 


mr*[V + 2u/r,.V]N(V, 8, 4), 


where 1, is the radius of the star, and y/r? is the gravitational intensity which it 
produces. If we assume a velocity distribution at infinity of the form 


(277)-*poh® exp [ — 4h>(V? + c?— 2Vccos)], 


where c is the mean velocity of the star in the cloud, then 


N(V, 6, b) =(27)-*poh? exp [ — 44°(V2 + 2 — 2Vccos6)]V? sin 8, 
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and the total rate of accretion is 
mrHany ph | | | g(t 2H/r0¥) 
x exp [— $h7(V?+c?—2Vccos6@)|V?sin6dV dé dd, 
12 70 R 
= () ro"poh® | | (V? + 2/19) exp [— $h°(V2 +c? —2Vccos6)|V sin 6 dV dé, 
en “0 “0 


5 (h?Vc eel 


n= 


exp [— $h?(V?+ c? —-2Vccos6)] =exp[—}h(V?2+e 
P ( 


so the integral becomes 


1/2 
(=) rotol? exp (—sa%e) 


« 2np-n -o cr 
= | | exp(—3h8V2)(V""8 + 24 V"*1/r,) cos" 0 sin dV dB. 


n=0 n! “0 





This may be integrated term by term, and reduces to 
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As before, we put p=he, and qo? = 2uh?/r,, and then the rate of accretion becomes 


(277)"?r5? poh [2f,(P) + qo7fo()].- 


For small values of p, f,; and f, may be expressed as 
I p ; : — {)"p2n 
SilP)= 5 5 {exp —}p*)+(p+p") [" exp (— 1y*)ayl = a Biss ol 


n\(4n*—1)’ 
— 7\"p2n 
fp)= = exp(— 49") dy= een ame 


2"n!(2n+1) 
if p is large, then we have 


f(D) =(4n)"" (r+ =) i - aoe 


fp) = (47)? 57 pexP(— 8 »p*)+.. 


For normal values of M, ry he h, qo >1 (see Table I), so the rate of accretion 
may be taken as 


(277)"*19°poh™"qo'fa(P), 


(877)"*MGropohf,(p). 
Graphs of h.f,(p) are given for various values of h and c in Fig. 4, and it will be 
seen that the introduction of the velocity distribution to Eddington’s model 
(which is equivalent to h~'=o) reduces the rate of accretion, and that this rate 
is a maximum for p=o, or c=o, when it is given by (87)"*MGr p,h. This is 
the rate of increase of the mass of the star due to accretion, so we may take 


u/>- (87)2(Grohpy) 


Values of this function for standard spectral types are given in Table V. The 
quantities r, for the different spectral types are those provided by Allen (1) 
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Fic. 4.—The rate of accretion, measured by the function hf, (jw), as a function of the mean speed, 
c, for different values of the velocity-dispersion h-*. 
TABLE V 
; ' . dM . 
Inverse rates of accretion for different spectral types. The quantity tabulated is M1 | ——, in 
€ 
years, and the separate columns refer to supergiants, giants and dwarfs 
Sp S-—g g d 
Os ew 6-7 x 101 ) 
Bo 6-7 x 101% 1°5 X 10!? 
Bs 4:3 10" 3°5 X t0"" Stars in 
Ao 3°4 x 1018 6:0 « 1017 . H 11 regions 
5 2°7 x 10° 7°7 x10" 
Fo 2°1 X 1016 g*1 x 1017 J 
Fo 1°3 X10! 5°5 x 1015 7 
Fs5 10 X 1014 6°7 x 1015 
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These results show that if collisions are neglected, then Eddington’s conclusion 
that in general accretion is quite insignificant is fully confirmed. They also 
suggest that, since the case where there is spherical symmetry (i.e. c=0) gives 
the maximum rate of accretion, a more exact treatment, including collision 
effects, would be useful for this case. 


Appendix.—Collision frequencies in a rarefied gas at rest in a gravitational field. 
—We have seen that a gas in which collisions can be effectively neglected and 
which has a velocity distribution function at infinity of the form 

poh®(27) >? exp { — $h°(X? + V2 + Z?)}, 
will have a distribution function of the form 
poh®(27)-3” exp { — $h2(X2 + V2 + Z?2—2y/r)}, (X?+ Y?+Z?-2y/r)>0, (12) 
at a distance r from the centre of attraction. This result has been obtained on 
the assumption that collisions never occur. However, if the molecules are of 
finite size, there will be collisions, though, for low densities, they may be very 
infrequent. We shall therefore examine the collision frequencies when the 
velocity-distribution is given by (12). The methods used are standard ones (11), 
but of course the velocity-distribution is not Maxwellian. 

We seek both the total number of collisions per unit volume per unit time, 
and also the total number per unit volume per unit time of those collisions in 
which the speed of one particle is reduced to a value less than (2Ap/r)"*, A<r. 

Let us specify two particles, A and B, of speeds W and W’, which collide 
in such a way that their line of centres lies in a certain small solid angle ’w. The 
line of centres is directed in such a way that the direction of motion of A is inclined 
to it at an angle «, o<a<z7/2. The direction of motion of B is inclined at B 
to the line of centres, whereo<B<z. Then if ois the diameter of the particles, 
the number of particles such as B which collide with A in the chosen solid angle 
in time df is o7Sw|W cos x— W’ cos B|dt.f(W’), where 

f(W’) = poh'(27)" exp { — )h2(W"? — 2p/r)}W2 dW’ sin B dB, 
and W">2y/r. This expression is integrated over dw, B and W to give the 
probable number of collisions in which A is involved in a time interval dt. The 
integrations are not difficult, but cases such as W’ = Wcos« have to be treated 
separately. We find that the probable number of collisions is 
2) 
if 





, ‘ 2 
(827)'*07Stpyh { W? cos? « + ~ + 
for W cos « <(2u/r)"*, and 


12 
(827)! 207Stp,h {, W cos x (*#) + 2h exp [— }h?(W? cos? « — 2u/r)] 
; -W cos x 
+2Wcos« | exp [— $h>(W’? — 2/r)] aw’\, 
. (2 rye J 
tor Wcosa>(2/r)'*. ‘To obtain the total number of collisions per unit volume 
per unit time, we multiply these expressions by the number of particles per unit 
volume such as A (namely poh*(27)-'* exp { — $h2(W? — 2u/r)}W? dWsina da, 
W?>2y/r) and integrate over « and W. In so doing, of course, every collision 
is counted twice. Again the integrations are tedious, but we finally obtain for the 


total number of collisions per unit volume per unit time the quantity 


socorh Leggs [” exp(g—s)de} (13) 


qd 
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where g?=2ph?/r. As r— oo, q->0, and the expression tends to the Maxwellian 
value of 27'/20%p,2h-!. 

In very similar fashion we can calculate the number of collisions per unit 
volume per unit time in which the speed of one particle is reduced to a value 
less than (2Au/r)!?, A<1. Ina perfectly elastic collision between A and B the 
components of velocity perpendicular to the line of centres are unaltered, and 
the components along the line of centres are interchanged. It follows that if 4 
is to lose kinetic energy of amount $U? per unit mass, then 


W? cos? « — W”* cos? B = U?. (14) 


Thus the probability that A will be involved in a collision in unit time in which 
it loses energy of amount between }U? and }U?+ U8U per unit mass is 


o75w | poh*(277)-"? exp { — $h?( W’? — 2u/r)}|W cos «— W’ cos B|W’? dW’ sin 8 5p, 


where f is given by (14), and 
W” sin Bcos B58 = U SU. 
This leads to probabilities 
a7 8wU 8 Upp h*(27)-*? 
W’ cos « — (W? cos? a — U?)!? 


(W? cos? « — U?)"2 W' dW’, 





x | exp { — $h>(W’? — 2y/r)} 


for particles which A overtakes, and 
o* 8wU 5 Up h*(27)-"? 
; Wcosa+(W?cos*a—U?)? 
* [ exp {- MW? — 2y/r)} Wate ott W'dW’, 
for particles which A meets. Thus the total probability is 
a W cosa ; ges 
207 $wU 5Up,h(27) '” CW costa— Ta | exp {— $4 W’?— 2y/r)}W’ dW”. 
In both integrals we require 
W"?>W?cos?a—U*, and W’?>2/r. 
Thus the probability is 


20% 8wU 8Upyh(20)-?? 








W cosa 
(W? cos? « — U?)!2 


’ 





if W*cos*a<U*+2y/r, or 


W cos a 
—1'2 
20° dw U 5 Up,h(27) (Wcosta— Uy exp | 
if W*cos*a>U*+2yu/r. The integration over all elements of solid angle d5w is 
trivial. Next we multiply by 


poh®(27)-" exp { — $h?(W? — 2u/r)}W?8Wsin « da, 





— sh>(W? cos? « — U? — 2yu/r}, 


which is the number of particles per unit volume with speed W inclined to the 
selected direction at an angle «. Integrating over «, we obtain the total number 
of collisions per unit volume per unit time in which particles of speed W are 
reduced in speed to (W?— U)"?._ The answer will depend on whether W?> or 
<U?*+2p/r; however, we are going to apply these results to the case in which 
a capture is effected, so that W?< U? + 2y/r. It follows that for all a ino <a <7/2 


W? cos? « < U? + 2y/r. 





saws 
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The number of collisions of the kind described is therefore 
“cos ~! (U/W) W cosa 


2=0) (W? cos? « — U2)? 





407U 5Up,*h' exp { — $h?(W?— 2 /r)}W2dW | sin « da, 
or 
40° U 8Up,7h' exp | — 1h2(W? — 2u/r)} WSW [ —(W? cos? a — U?)¥2] ces 9 O/) 
or 
40°U 8Up,7h' exp | — $h?/(W? — 2u/r)}W(W? — U*)!28W. 

‘Yo tind the number of collisions per unit volume and per unit time in which 
particles of speed W are reduced in speed to a value less than (2Au/r)"?, we 
integrate with regard to U, for values of U lying between (W?— 2Au/r)"? and W. 
This gives 

<0%p,2h* exp | — }h?(W? — 2y/r)}W(2Au/r)>? dW. 

Finally we integrate over all values of W, to get the total number of collisions 
per unit volume and per unit time in which the energy of one of the particles is 
reduced to Au/r per unit mass, or less. The number is found to be $07p97h*(2Au/r)>?, 
or, in terms of the parameter gq, 

50° py *h 82g. (15) 

A comparison between (13) and (15) shows that the first term in (13) gives 
the number of collisions which lead to captures (A= 1), and the remaining terms 
give the number of collisions which do not lead to captures. Only for g=2 and 
larger values is the proportion of collisions leading to captures significant. Even 
then, as is clear from (15), the distribution with regard to kinetic energy amongst 
those particles which are captured favours energies just below the energy of 
escape rather than small values. 


Department of Mathematics, 
The University, 
Manchester, 13: 
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Summary 


In a series of papers which are now being issued by the Astronomical 
Institute of the Ankara University (18) the autocorrelation and subsequent 
power analysis is applied to large numbers of RR Lyrae and cepheid 
variables in the Milky Way cloud, globular clusters, etc. ‘The present paper 
contains a brief outline of the method which was used in these papers and also a 
brief discussion of the significance of the intensities of the overtones obtained in 
that way. For statistical purposes only the relative intensities A(1) of the first 
overtone are important. The values A(1) are plotted against the logarithms 
of the period. With periods P<10' there is a general decrease of A(1) 
with increasing period, while with periods P~10' A(1) increases with the 
period. There is a tendency of the points A(1) to concentrate along certain 
lines, but it is difficult to say whether these lines indicate certain levels or 
correspond to the limits of fairly broad belts in which there is a concentration 
of the points A(1). With the shorter periods there is a discontinuity near 
log P 0°30, where a sudden increase of the values A(1) occurs. This 
discontinuity could be related with a transition from Bailey’s subtype ¢ to a—b 
orvice versa. The points A(1) corresponding to the RW Draconis variables are 
confined to a narrow vertical belt around log P 0°30. Some systematic 
differences seem to occur between the distribution of the points A(1) in the 
different galactic star clouds and the Magellanic Clouds. 





1. Outline of method.—'The discussions in this article are based on the results 
of an autocorrelation analysis of the light curves of different groups of cluster 
variables and of cepheids. The autocorrelation method for analysing time 
series has been developed by Kendall (16). For a time series with equidistant 





terms u), Us, -.., u, for which Xu,=o, the autocorrelation coefficient r, is 
defined by 
7. Luu; +4 (1 1) 
a ys Ly 9 2 . 
(Su. Su; ..2)" 


where all summations are from o to N—k. Here we use the notation introduced 
by Ashbrook, Duncombe and van Woerkom (2) who give the convenient 
approximation 


Vy, = (Luu; ,4)/ (Lu? — ku,?} (1.2) 
which is adequate if N>k. The disadvantage of the correlogram is to mask 
a secondary periodicity. ‘To overcome this disability the correlogram can be 
transformed into the corresponding power spectrum. By Kintchine’s theorem 


(6) the power spectrum is the Fourier transform of the correlation function. 
It is computed from the relation 


am f)= = Er, cos 2nfk, (1.3) 


where f is the reciprocal of the trial period and x(f) indicates the mean squared 
amplitude of each frequency. 
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As a rule the requirements of extent, continuity and reasonable homogeneity 
of the light curve to be analysed are difficult to meet. With the RR Lyrae stars 
and the cepheids the light curve remains almost unchanged, period after period. 
So it was supposed that the light curve continually repeats itself and as a 
consequence the time series can be considered to be infinite. The period of the 
correlogram is the period of the variable. Obviously the introduction of a constant 
period and light curve imposes some restrictions. The power series is converted 
in something very much like a Fourier series. A consequence of our assumption 
is that for any value of f which is not commensurate with the fundamental period 
we have 7(f)=0. From the relation (1.3) we obtain the squares of the intensities 
of the frequencies v, 2v, 3v, ..., that is, of the fundamental and its overtones. 
These correspond to x(0), 7(1), 7(2), etc. In this paper we use the relative 
intensities of the overtones, e.g. the values 

A(fy= {2D 
(7(0) ) 

With all light curves which were analysed the series A(o), A(1), A(2), ..., is 
represented by a sequence of numbers with rapidly decreasing values. ‘The only 
term which is of interest for a general statistical study is the term A(1). The 
terms of higher order were found to be of the same order of magnitude as the 
probable error in their determination. By this method the amplitudes of all 
light curves are reduced to a uniform scale, while the shape of the light curve is 
described in terms of the numerical values of a few parameters. Starting from 
the assumption that the atmosphere of a star is affected by strictly periodic 
disturbances, while each disturbance sets up a damped oscillation (23), it is not 
difficult to construct artificial light curves which closely resemble that of a cepheid. 
With such an artificial light curve only one period is present but in the analysis 
several terms A(f) of non-zero value will turn up. The description of a light 
curve in terms of the numerical values of A(/) therefore is a purely formal one. 
It is not necessary that the overtones have any physical importance. Finally the 
following remarks should be made: 

(i) On the one hand it is an advantage that in this method all light curves 
are reduced to a uniform scale. On the other hand the amplitude of a light 
curve is a valuable part of the observational evidence and this part is sup- 
pressed. 

(ii) The light curve was supposed indefinitely to repeat itself so that small 
variations in the period and in the shape of the light curve are suppressed. 
Especially with the variables of the RW Draconis type this is a disadvantage. 
Several authors consider these stars to be affected by two vibrations of almost 
equal periods resulting in the occurrence of a beat period. Within each beat 
period the shape of the light curve varies considerably. When possible the 
light curves corresponding to different phases of the beat period were analysed 
and the range of the values A(f) was determined. The periods of the RW Draconis 
variables are around } day and only a few have been observed sufficiently 
to allow of a complete analysis. For some additional variables light curves for two 
different epochs are available, so that part of the range in A(1)is known. For 
the majority of the RW Draconis variables only certain mean values A(f) could 
be obtained. 

(iii) While the light curve completely determines the correlogram, the 
correlogram does not uniquely determine the light curve. 
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2. A survey of data.—Using the method described in the preceding section, 
large numbers of variables were analysed. Here only a short survey is given 
of the different groups of variables which were analysed with a reference to the 
sources from which the light curves were obtained. 

(a) Variables in globular clusters ——A large number of variables in w Cen, 
for which the light curves were taken from the work of Martin (20(a)), was 
investigated. First, a number of variables was selected at random but in such 
a way that all periods are represented. In addition, a number of variables was 
analysed which either have: 

(i) a very short period ; 

(ii) a comparatively long period; or 

(iii) a period near the transition from the Baily subtype ¢ to a—6. 

All stars with P>1" have been included. In the same way a number of 
variables in M5 were analysed for which the light curves were determined by 
Oosterhoff (22(a)). In M3 forty variables were analysed for which the light 
curves given by Hett (14) or by Greenstein (11) were used, except for those 
variables for which Martin (20(b)) has given improved light curves. Next, 
all in M53 were analysed. As a rule the light curves given by Grosse (12) were 
used, but for some the improved ones as determined by Oosterhoff (22(c)). 
Of the 86 variables which Wright (33) has studied in NGC 3201, thirty-eight 
were used and also all variables which this same author has studied in NGC 4833. 


For the clusters Mg, M22, NGC 6218, NGC 6254, M13, M80, NGC 5053," 


NGC 5466 and M2, the light curves given by Sawyer (26) were analysed. ‘The 
curves which Sawyer indicates as uncertain were excluded. 

(6) Variables in the Magellanic Clouds._-For the Small Magellanic Cloud 
Shapley, Carlston and Nail (28) have given a sequence of light curves for different 
periods which they consider as being representative for this Cloud. For the 
Large Cloud a similar sequence is given by Shapley and Nail (29). All variables 
contained in these two sequences were analysed except those for which the 
authors indicate that there was a faint companion in contact with the variable 
or a doubling which distorted the image and decreased the range. A few variables 
were included of which for special feasons Shapley (27) has separately published 
the light curves. 

(c) Variables in galactic star clouds.—In the Sagittarius cloud the values A(/) 
were derived for 70 RRLyrae stars and cepheids. The light curves were 
taken from the publications by Uitterdijk (32), Martin (20), Swope (31) 
and Ferwerda (7). From the light curve given by van Gent (9) for the variables 
in or near the constellation Corona Australis fifty were analysed. For a large 
number of variables in or near the constellation Carina, the light curves were 
determined by Plaut (24), Gaposchkin (8) and Payne-Gaposchkin (24). For 
fifty-five of these variables the values A(f) were determined. 

(d) Galactic cepheids—The analysis was applied to the variables in lists 
collected by Hertzsprung (13) and Joy (15). All thirty-eight variables contained 
in the list of Hertzsprung were used. From the list of Joy fifty-five systems 
were selected in such a way that all periods are represented. In addition to these 
the light curves based on “‘ mean coordinates’’ were considered. These ‘‘ mean 
coordinates” given by Payne-Gaposchkin (24) for different periods represent 


the mean of light curves which previously had all been reduced to a uniform 
amplitude. 
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(e) RW Draconis variables.—An analysis of the light curve at all phases % 
of the beat period was possible only for AR Her, which variable has been 
extensively observed by Balasz and Detre (3). For the RW Draconis variables 
in w Cen Martin (20) gives the light curves at two or three epochs. All of these 
were analysed. Consequently for these variables part of the range over which 
A(f) varies is known. Moreover all stars were analysed which in the lists 
mentioned under (a), (4), (c) and (d) are indicated to be or are suspected to be 
of the RW Draconis type. 

In some cases one and the same variable occurs in several of the lists 
mentioned in (a), (4), (c), etc. So we have several systems for which the light 
curve has been determined by at least two different authors. ‘This makes it 
possible to evaluate the probable error in the determination of A(f). Here only 
the probable error in the numerical value of A(1) is given, viz. p= + 0-040. 
With recent light curves the value of p may be smaller, but as our lists contain 
many older light curves here the value p= + 0-040 must be used. ‘This value 
includes influences which arise from the differences between visual and photo- 
graphic curves, from instrumental differences and from the way in which the 
mean light curve was drawn through the observed points. With the majority 
of the variables the numerical values of A(z), A(3) were found to be of the same 
order of magnitude as p. Therefore in the following only the values A(1) are 
used. In the Communications of the Astronomical Institute of Ankara the values 
A(1) have been tabulated for several hundreds of variables (18). ‘There is no 
reason to repeat this tabulation, and for the numerical values of A(1) the reader 
is referred to the original lists. 

3. The relation between A(1) and log P.—Fig. 1 gives a plot of the values A(1) 
against log P. All variables contained in the lists mentioned in Section 2 were 
used with the exception of: 

(a) the irregular and RW Draconis stars ; 

(5) two stars with period around 100 days which fall outside the limits of the 
figure (one in the Magellanic Clouds, one in Ms); 

(c) the values A(1) based on ‘‘ mean coordinates ”’. 

In Fig. 2 the open circles correspond to the values A(1) derived from the 
mean coordinates given by Payne-Gaposchkin (24). Black dots are cepheids 
with a distance over 250 parsecs from the galactic plane which according to 
Payne-Gaposchkin (/oc. cit.) must be considered as Population 11 stars. Triangles 
correspond to cepheids in globular clusters. If the shape of the light curve is 
well determined, black triangles were used. If the shape of the light curve was 
uncertain an open triangle was used. From Figs. 1 and 2 it is apparent that the 
distribution of the points A(1) in the correlation plane is not an haphazard one. 
The points A(1) seem to be concentrated along certain levels. ‘These different 
levels are poorly defined and it is dubious whether they must be considered as 
being levels in a strict sense or merely as the limits of fairly wide belts of preference. 
It would not be difficult to draw another set of lines a little higher or lower in the 
figure, or witha slightly different inclination, which are equally probable. However, 
provided that their relative value is kept in mind, these levels are convenient for 
discussing the main features of Figs. 1 and 2 and later on will enable us to compare 
the distribution of the points A(1) within the different groups of stars. From 
Fig. 1 it appears that, for periods P<104, A(1) decreases along all levels with 
increasing period. For evident reason there is an upper limit for the values A(1) 
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which in our diagram is around +0-60. Near P=10% small values of A(1) 
prevail. In other words here the light curves resemble a pure sine curve. 


At) — 





40°- 


SOUrwe Ieee eS | 


biti piy 


8 








. 
-10 log P 
Fic. 1.—In this diagram all available values A(1) are plotted against their corresponding 
value log P. The variables seem to be arranged along a number of levels which do not all have 
physical significance. The period P=10 indicates a reversal. Before this A(1) decreases with 
increasing P. Beyond P=10 the value A(1) increases with increasing P. Except the levels c, all 
levels to the left of P= 10 start in the upper left part of the diagram and can be traced up to P= 10. 
There are no discontinuities except near log P= — 0-3. 
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bic. 2.—In this diagram the large open circles indicate the values A(1) derived from the mean 
coordinates of light curves. Black triangles correspond to cluster cepheids with well determined 
light curves, open triangles to those with poorly determined curves. Black dots are cepheids with 
distance over 250 parsecs from galactic plane. The shaded area corresponds to the region occupied 
by the RR Lyrae stars. 





With the periods P>10% there is a general tendency of the values A(1) to 
increase with period. With the periods P<104 a discontinuity occurs around 
log P= —0-30. ‘The majority of the stars for which log P< — 0-30 are arranged 
along the levels c, while those for which log P>0-30 are along the levels a. The 
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stars on levels c belong to Bailey’s subtype c, those on a, and a, to the Bailey 
subtypes a and 5. 

These irregularities in the distribution of the points A(1) correspond to the 
period intervals detected by Ludendorff (19) and Hertzsprung (13). As appears 
from the set of curves collected by Campbell and Jacchia (5), each interval sets 
in with asymmetrical light curves (high value of A(1)) and ends with symmetrical 
light curves. From our diagram it would appear that this trend is reversed with 
the longer periods. In Fig. 1 the final interval sets in around log P= 1-0 with 
low values of A(1), while the asymmetry increases with increasing period. 

Fig. 2 leads to identical conclusions. In this figure the variables of Bailey’s 
subtype c have not been included. The cross-hatched area in the upper left part 
of the diagram indicates the area which is densely crowded by variables of the 
subtypes a and 6. From Fig. 2 it appears that the levels a, and a, (the belt a,—a,) 
are typical for variables of the Population 11 type, but some P11 variables also 
occur along 6, and 6,. Conversely the majority of the variables on the levels 5, 
and 54, (belt d,) is of Population 1 type. 

For periods P>10" the situation is not clear. ‘The levels d,, d, and d, are 
very indistinct. In Fig. 1 the majority of the stars are along d, and in Fig. 2 the 
points corresponding to the “‘ mean coordinates” are also along ds. 

Froma Fourier analysis of the light curves of ’ Cepheid stars Payne-Gaposchkin 
(24(d)) finds that up to periods 164 there is a difference between the variables of 
Population 1 and those of Population 1. From Fig. 2 it would appear that in 
the mean with P>10¢ values A(1) for Population 11 variables are smaller than 
those of Population 1 variables, but the evidence is not conclusive. All variables 
along d, are of type 11. But along d, and d, the distribution of the values A(1) 
seems to be the same for the two population types. ‘This may be due partly to 
the fact that we have included some dubious values in our material (the open 
triangles). 

4. The RW Draconis stars and the irregular variables.—'The values A(1) 
corresponding to the irregular variables and RW Draconis stars are plotted 
in Fig. 3. ‘The levels c3, a, and a, are borrowed from Fig. 1. From Fig. 3 it 
is evident that the RW Draconis stars are confined to a vertical belt which is 
roughly limited by the values log P= —o-40 and —o-20. This vertical belt 
therefore also roughly coincides with the discontinuity in the distribution of 
the values A(1) around log P= —o-30. It might therefore tentatively be suggested 
that the RW Draconis stars are variables in a state of transition between the 
mode of vibration corresponding to the levels c and the one corresponding 
to the levels a. Temporarily both vibrations are present. From the work of 
Schwarzschild (30) and Kluyver (17) we know that the simultaneous occurrence 
of two vibrations is possible. Also, in their study of the light curve of AR Her, 
Balasz and Detre (3) have shown that its complicated variations can be explained 
if it is accepted that simultaneously two vibrations are present. ‘The RW Draconis 
stars and the irregular variables are not evenly distributed over the whole interval 
between the levels c and a. ‘This may be a question of selection, the amplitude 
of the variables on or near the c levels being much smaller than of those near 
the levels a so that the latter ones are more easily discovered. 

Further evidence that with the RW Draconis stars and the irregular stars 
we may have to deal with variables in a state of transition is given by Fig. 4. 
In this latter figure we have indicated the intervals in log P which with the various 
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Fic. 3.—The irregular variables of the RW Draconis type. The arrows indicate the range 


of the values A(1) through the secondary cycle. The irregular variables are confined to a narrow 
interval of log P and may represent a group which is transitional between c and a. 





groups of stars are occupied by the variables on the levels c and those on the 
levels a. For the first the intervals are indicated by dotted rectangles and for 
the second by black rectangles. The cross-hatched vertical belt indicates the 
interval in which the RW Draconis and irregular variables occur. From the 
figure it is evident that with all groups of stars the end of the levels ¢ and the 
starting point of the levels a is in or at least very near the interval occupied by 
the irregular variables. With the Carina cloud we have not been able to find 
any variables which could be identified with the levels c. 
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Fic. 4.—The black rectangles indicate the range in period of the levels a, the spotted rectangles 
that of the levels c. The cross-hatched strip indicates the period range of the RW Draconis stars. 
For Carina the results are uncertain. It is noticed that in the Milky Way clouds the a levels extend 


to far shorter periods than in the clusters. Ali levels a and c start or end near the interval 
occupied by the irregular variables. 























fr wie ASA Hv oe ES. 





TK RTS 


No. I, 1957 A statistical study of pulsating stars 79 


5. The distribution of the points A(1) in Carina, Sagittarius and the Magellani: 
Clouds.—When plotted separately some typical differences appear between the 
distributions of the values A(1) in the Carina cloud, the Sagittarius cloud and the 
Magellanic Clouds. 

For the Sagittarius cloud the distribution of the values A(1) is given in Fig. 5 
and that for the Carina cloud and the Magellanic Clouds in Figs. 6 and 7 
respectively. In this case only the variables with P>14 are used. From the 
figures it would appear that in the Sagittarius cloud the relative frequency of 
the cepheids with P> 104 is slightly larger than in the Carina and the Magellanic 
Clouds. 

In Sagittarius the levels 6, and 6, (the belt d,) end abruptly towards the left, 
but in Carina they extend to much shorter periods. ‘This is even more so with 
the cepheids in the Magellanic Clouds. The large extension of the cepheids 
towards short periods, especially in the Small Cloud, has been discussed by 
Shapley (28). This abrupt end of the levels 5,, 5, in Sagittarius may explain 
why the frequency of cepheids with P< 104 is relatively small in this cloud. 
Between the distributions prevailing in Sagittarius and those in Carina and the 
Magellanic Clouds there is a second difference. In Sagittarius several variables 
are on the level d,. On the other hand, in the Carina cloud and in the Magellanic 
Clouds no variables occur which can be identified with the level d,._ It is evident 
that the differences mentioned before are in the nature of small shifts. But it is 
interesting to notice that the stars in Sagittarius are considered to represent 
Population 11 and those in Carina and the Magellanic Clouds Population 1 
stars. 

From Fig. 2 it appeared that the variables corresponding to Population 
11 have a preference for the levels a and d, (see Section 3). This might 
explain why the level d, occurs in Sagittarius but is absent both in the Carina 
and in the Magellanic Clouds. Payne-Gaposchkin (24) has pointed out that 
there is a much closer resemblance between the cepheids in Carina and those 
in the Magellanic Clouds than between the cepheids in Carina and those around 
Sagittarius. 

This conclusion is confirmed by the distribution of the points A(1) in 
Figs. 5,6 and 7. There are some differences between the distributions in Carina 
and the Magellanic Clouds. When we consider the cepheids with P<10%, it 
is evident that with the Magellanic Cloud variables the values A(1) in the mean 
are slightly larger than with those in Carina. 

6. The different levels —When considering Fig. 1 it was remarked that the 
different levels are badly determined and that it is doubtful whether we have 
to deal with real levels or simply with certain regions or belts of preference. 
In the subsequent figures the levels were convenient for discussing the general 
distribution of the values A(1) but in these additional diagrams also the levels 
are poorly determined. We must therefore consider the question whether or 
not any real value can be attached to them. 

We first consider the levels c. These are the only levels along which the 
variables of Bailey’s subtype c occur. From Fig. 1 it might appear that there 
is a faint preference for three different sublevels c,, c, and c,, but c, is especially 
badly determined. 

In view of the large probable error in the determination of the values A(1) 
it is dangerous to attach much importance to these three levels. But, according 
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Fic. 5.—The variables in the Sagittarius cloud and surroundings are plotted. The leveis 
3, and 8, are separated but poorly developed. The presence of dg should be noticed. The galactic 
longitude is from l= 310°-350 . 
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Fic. 6.—The variables in and near Carina are plotted separately. The levels 6, and 8, are 
more pronounced but are not separated. The level dg is missing and perhaps also d2. 
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to Martin, the c-type variables in w Cen fall apart into two different classes 
which he indicates as the classes c and c’. It appears that all variables to which 
he assigns the class c’ are on the level c,. Of the variables which he indicates 
as being of type c only a few are on the level c., the majority being along c, and c;. 

The variables along the levels a are mainly those belonging to the Bailey 
subclasses a and 6 and the cepheids occurring in globular clusters. ‘lhe levels a 
are Clearly distinct from the levels c. Only in some globular clusters do these 
variables seem to fall apart into two different groups. Of this no trace remains 
in the composite diagram (Fig. 2). 

The variables of the Bailey subtype a are found along both levels a, and a, 
and this is also the case with the variables of the subtype 6. So it is improbable 
that the levels a, and a, are real. From a study of the period—amplitude relation 
in the outer regions of the globular clusters w Cen, M3 and Ms, Belserene (4) 
concludes that the period-amplitude diagram suggests that there are two groups 
of stars, each represented in the different clusters but with a relative population 
that differs markedly from cluster to cluster. 

From our study of the distribution of the values A(1) in the different clusters 
a similar result was obtained. ‘The paper from Belserene does not indicate the 
individual variables which belong to the one group or the other and a direct 
comparison is not possible. 

Consequently we will consider the variables along a, and a, as members of 
one and the same group, but will admit the possibility that further research 
may reveal the existence of two different levels. From Figs. 1 and 2 it is evident 
that this group can be traced almost up to the period of reversal at P= 104, 
Apparently this group is typical for stars of Population type u. However, a 
few stars also of the Carina cloud and the Magellanic Clouds seem to belong 
to that group. This suggests that close relation may exist between these stars 
and the cluster variables. The levels 5, and 6, (the belt d,) are well developed 
with the stars in the Carina cloud and the Magellanic Clouds (type | stars). 
In Sagittarius the levels 5, and 6, are short, but it is only in this cloud that there 
is a separation between 5, and 6,. In the other clouds 8, and 4, overlap each 
other. The stars along 5, and 6, are to be taken together as a group which is 
indicated by the belt d,. ‘The belt d, seems to be typical for Population | stars. 
For the periods P>104 the level d, is typical for the variables of Population 11. 
Only cepheids in globular clusters and some cepheids in Sagittarius fall along d, 
and not one variable from the other groups which we have considered. So it 
seems possible that d, represents a separate level. 

Along d, and d, the scatter of the individual values A(1) is so large and the 
separation of d, and d, in the different diagrams so unsatisfactory, that these 
variables must be considered as one group. ‘The group (d,, ds) contains variables 
from both population types and there are no indications that within this group 
the distributions of the two types are different. 

It has been suggested that whole groups of variables may vibrate in an overtone 
rather than in the fundamental. Of this Schwarzschild (30) thinks the variables 
in M3 to be a notable example. 

Also Arp, Baum and Sandage (1) conclude that there is no unique 
period—luminosity relationship for the type 1 cepheids as exists for the 
classical cepheids. The points indicating the absolute magnitudes of type II 
cepheids fall along a number of lines which are more or less parallel to that 
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indicating the period—luminosity relation for classical cepheids. ‘hey can 
be made to coincide with the curve for classical cepheids by a horizontal shift 
over a distance log2, 2log2 or 3log2, suggesting that the variables vibrate 
in the first, second or third overtone. We have considered the possibility that 
the apparent arrangement of the points A(1) along different levels is due to a 
similar effect. 

If this were the case the period—-luminosity relations corresponding to the 
different levels would not be identical. This possibility is tested in Fig. 8 where, 
for all cepheids in the Magellanic Clouds for which values A(1) are available, the 
apparent magnitude is plotted against the logarithm of the period. 
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Different symbols were used for indicating the variables on the different 
levels A(1). ‘Those on the levels a are indicated by crosses, those on or near 
5, and 6, by open circles and black dots respectively. The same symbols have 
been used for indicating the variables along d, and d,._ It was supposed that both 
clouds have the same distance modulus. The results of this test are negative. 
The diagram merely indicates the mean period—luminosity relation, which has 
been given repeatedly by Shapley, and for which he has discussed the large 
scatter of the individual points around the mean curve. 

Between the distributions of the different groups of variables in Fig. 8 no 
systematic differences are apparent. This confirms that the variables along 3, 
and 6, must be considered as one group and also the variables along d, and d3. 

The variables which we have used in Fig. 8 are Population 1 cepheids. 
Consequently our negative results do not invalidate the conclusions which 
the aforementioned authors have obtained for the variables in globular clusters. 
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7. Final remarks and summary.—The material upon which Figs. 1 and 2 
are based is far from homogeneous. Therefore an attempt to give an explanation 
of the general distribution of the points A(1) in the correlation plane would be 
premature. We therefore only give a brief summary of those results which 
seem rather positive. 

(a) With the variables with P< 104 the asymmetry of the light curve decreases 
with increasing period; with those for which P>104 the asymmetry increases 
with increasing period. Apparently around P= 104 a point of reversal occurs. 

(6) Around log P= —o0-30 a discontinuity occurs. Here there is a sudden 
increase of the asymmetry. ‘This discontinuity may be related to a transition 
from the variables of type c to the types a and b, while the RW Draconis variables 
may be stars in a transitional stage. 

(c) With periods log P< —0-30 the variables are arranged along the levels c 
and there is a possibility that we have to deal with at least two different levels. 

(d) With the periods 1 >log P> —o0-30 the variables are arranged along the 
belts a and d,. A division in sublevels has no physical reality. The levels a 
are mainly occupied by variables taken from groups of stars which are considered 
to represent the Population type 1. ‘The levels 5 are mainly occupied by variables 
from groups of stars representing the Population type 1. However, a few variables 
of type I also are on the levels 5 while conversely a few type I variables are on the 
levels a. 

(e) With periods log P>1 the level d, is exclusively occupied by variables 
of Population type 11. ‘Thus d, may represent a separate level. Along d, and d, 
variables occur from both population types. ‘The subdivision of these variables 
into the separate levels d, and d, does not seem to be real. 


Astronomical Institute, 
Ankara: 
1956 December. 
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RADIUS VARIATION AND POPULATION ‘TYPE OF CEPHEID 
VARIABLES 


A. W. Rodgers 


(Communicated by the Commonwealth Astronomer) 
(Received 1956 December 4)* 


Summary 


The method of Wesselink for the determination of the radii of pulsating 
variables is applied to the southern cepheids« Pavonis, f Doradus and / Carinae. 
‘The radius of 8 Doradus is found to be abnormally large for its period. The 
results are combined with the radii of other cepheids found by previous workers 
using the same method and a plot of relative radius amplitude against period 
is made. __I[t is apparent that, for classical cepheids, the light variations are 
mainly due to changes in effective temperature, and the relative radius 
amplitude increases slowly with period. 

However, for the type 11 cepheids, the relative radius amplitude increases 
rapidly with period until for stars of the longest period W Virginis type the 
variation in radius is of the order of the mean radius of the star itself. 





1. Introduction.—A method, due in its final form to Wesselink (1), for the 
determination of absolute diameters of cepheids from a combination of photo- 
metric and radial velocity data has been applied in recent years to 5 Cephei, 
7 Aquilae and RR Lyrae by Stebbins (2), and to Y Ophiuchi (3) and W Virginis (4) 
by Abt. 

Recent radial velocity determinations by Stibbs (5) and photometric work by 
Eggen and Gascoigne (6) have made it possible to extend this work to the three 
southern cepheids « Pavonis, 8 Doradus and / Carinae. « Pavonis (g-1 days) was 
selected because there are a number of indications that it belongs to Population 
type 1 (in which case incidentally it is easily the brightest of the type 11 cepheids) ; 
8 Doradus because of its unusual light curve, high galactic latitude and period of 
g°8 days, and / Carinae because it has a much longer period (35-5 days) than any 
other cepheid for which a radius has been determined. 

The validity of Wesselink’s method rests on the basic assumptions (a) that at 
two points on the light curve where the colours are equal the effective temperatures 
are also equal; and (6) that actual radial displacements may be obtained from the 
observations of radial velocity. More recently Cavanaggia and Pecker (7) have 
presented an exhaustive discussion of the photometric and radial velocity data as a 
test of the pulsation hypothesis. Using the results of their computations on 
supergiant model atmospheres they found disagreement with the radius variations 
found from Wesselink’s method. Whitney (8, 9) presented a similar analysis 
of 5 Cephei and » Aquilae and found that, when corrected for the effects of line 
absorption, the use of model atmospheres gave results not inconsistent with those 
obtained by Wesselink’s method. 


* Received in original form 1956 May 8. 
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2. Radii Determinations.—Following Stebbins (2), we may write the luminosity 
of the cepheid at any phase as 


Lp=47(Rmin+ SR} TA (1) 
‘Thus the ratio of light at any two phases of equal colour is given by: 
L, L.=(Rmin + 6R,)?/(Rmin +6R,)* (2) 


and if m,, m, are the observed magnitudes corresponding to L,, L, then the 
minimum radius of the star may be written 


Rmin=(85R,— A 5R,)/(A —5) (3) 


where A =antilog (m,—m,)/5. 

Observations of the three Cepheids, « Pav, 8 Dor and / Car were reduced 
according to equation (3). 

‘The photoelectric observations of the light and colour curves on the P, | 
system were made by Eggen and Gascoigne (6). ‘The radial velocity variations 
have been observed by Stibbs (5) at epochs overlapping the photoelectric observa- 
tions. Earlier, the radial velocity curve of 8 Dor was determined by Applegate 
(10), and Jacobsen (11) published velocity curves of « Pavonis and / Carinae. For 
all three stars the measures of Stibbs were consistent with the curves drawn through 
the earlier observations and confirmed that the form of the curves and systemic 
velocities were unchanged over the interval between the observations. ‘lhe 
combined observations from both epochs were used to define the form of the 
velocity curves. ‘The relative phases of the photometric and velocity data were 
found using Stibbs’ epochs of velocity minimum and the photometric ephemerides 
of Eggen and Gascoigne. An accuracy exceeding 0-05 P was thus obtainable in 
relating the two sets of curves. 

The radial velocities were integrated over the period to obtain the systemic 
velocity, and then cumulative integration at steps of 0-05 P was performed to obtain 
the radial displacements. A correction to the displacements by a factor of 24; i7 
(12) was made to allow for intermediate limb darkening and the projection effect. 
A summary of the data taken from the photometric and displacement curves is 
given in ‘l'able I and plotted in Fig. 1 against phase. 

For each star twenty pairs of phases of equal colour were taken and the radius 
obtained in the manner described above. A mean radius was formed, weighted 
according to the probable errors of the radial velocity and light curves, and the 
probable error of the weighted mean computed. 

For adjacent pairs of phases of equal colour having low weight, i.e. small light 
or radius differences in equation (3), systematic deviations from the final mean were 
noted. This effect was probably due to the fact that the light and radial displace- 
ments were read off smoothed curves drawn through the observed points. ‘lhe 
radii and other relevant data on these stars are given in Table 2. 

3- Results.—In this section we shall determine the absolute magnitudes of 
B Dor, « Pav and / Car and compare these with the values obtained for other 
galactic cepheids and cepheids in the Magellanic Clouds and M31. ‘The 
absolute magnitude is determined from equation (1), taking the solar radius, 
effective temperature and absolute magnitude as standard. ‘The radii we have 
obtained in the previous section; the effective temperatures we obtain from the 
observed spectral type for each star. 
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«x Pav 6 Dor 
Phase V P-V Radius V P-V Radius V 
< 10° km < 10° km 
©"00 3°97 + 0°39 = 15°03 3°44 O55 72°61 3°40 
0°05 4°05 0°44 «15°68 3°53 0°60 «73°41 3°28 
o'10 4°15 0°50 «61622 3°56 0°63 «74°20 3°30 
O'1s 4°27 0°56 «61666 3°63 0°67 3=74°83 3°38 
0°20 4°36 0°66 =17°00 3°69 o'7r 76°92 3°44 
0°25 4°42 O71 17°22 3°77, 076-7551 3°52 
0°30 4°52 0°75 17°31 3°84 080 75°35 3°59 
0°35 4°62 0°78 = 17°27 3°92 085 74°85 3°65 
0°40 4°69 080 «17°12 3°98 0°87) «74°15 3°72 
0°45 4°73 080 = 16°83 4°04 0-89 =673°23 3°79 
o'50 4°74 0°76 =616°40 4°04 0°88 = =72°12 3°84 
0°55 4°74 o-71 15°85 3°98 0°83 70°93 2°91 
0°60 4°70 0°64 15°16 3°92 O75 70°13 3°95 
0°65 4°63 0°57 14°38 3°85 0°69 69°80 3°98 
0°70 4°54 O°51 13°59 3°76 0°66 69°80 4°00 
0°75 4°42 0°45 12°88 3°68 0°65 70°25 4°01 
o'8o 4°22 0°37. «12°50 3°60 0°65 70°67 4°01 
0°85 4°04 0°32 «12°84 3°57 060 7JI'II 3°98 
0-90 3°98 O32 13°57 3°53 O49 «(71°54 3°87 
0°95 3°94 0°35 14°29 3°39 0°46 72°03 3°64 
0°00 3°97. +0°39 15°03 3°44 +0°55 72°61 3°40 
Taste Il 

Dimensions, luminosity and galactic data for « Pav, f Dor 
Star «x Pav 8B Dor 
R.A. (1900) 18h 40'9-6 5h 32-7 
Dec. (1900) 67° 22° 62° 33’ 
Galactic longitude 296 239 
Galactic latitude 26°3 32°7 
Period 4.9730 94842603 
Systemic radial velocity 38 km/sec 8 km/sec 
Spectral range F4-G5 F6o-G2 
Light range |’ 3°94-4°74 3°39-4°05 
Colour range, P-V’ -0°31-0'80 -~0°45-0°89 
Median radius 14°5 x 10°km 73x 10%km 
Median absolute photographic magnitude — 1°63 4°52 
Distance 116 pe 295-468 pc 


Population type 


TABLE | 


population type of cepheid variables 





Light, colour and radius variations in x Pav, B Dor and | Car 


/ Car 
P-V Radius 
10° km 
+ 0°95 84°00 
0°87 85°73 
oor =: O51 
1°00)=—- 93°67 
1°08 =—-97°65 
I*I5 10104 
1°22 103°77 
1°29 ©6105 89 
1°32 107°41 
1°35 108+23 
1°36 §=6108-36 
1°37 107°84 
1°37 106°67 
1°36) «= 104°81 
1°35 102°34 
1°33 99°39 
1°31 96-01 
1°27 92°37 
116 = 88-64 
1°O5 85°17 


+ 0°95 34°00 


and | Car 
/ Car 
gh 420-5 
62° 03’ 
251 
72 
359-5323 
3 km/sec 
F8-Go 
3°28-4°00 
089-1 °38 
96 « 10° km 
4°03 
179 pe 


I 


(a) B Doradus.—The mean minimum radius of 8 Dor was found to be 


(70+ 5)x 10° km. 


The radius at maximum light is 72-1 x 10*km, at which 


phase the star has a visual magnitude of 3-39 and a colour of +0-46 on the P, V 


system. 


Miss Applegate gave the spectral type, at light maximum, as F5 and Gratton 


(13) gives F6-7. 


Recent classification of spectra by Mme. A. de Vaucouleurs 


from plates obtained by Dr S. Gaposchkin at Mount Stromlo give the spectral 
type variation to be from F6la at light maximum to G2Ia,b at minimum on the 


revised Yerkes system. 


Of note is the very loose correlation found between 
the observed colour and spectral type variation throughout the period. 
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For the spectral type at light maximum, then, we adopt F6la, giving from 
Keenan and Morgan’s table (14) an effective temperature of 5800deg. K, 
assuming the la and Ib temperatures to be the same. ‘Taking the standard values 
of absolute bolometric magnitude, effective temperature.and radius from Allen’s 
book (15), we obtain a maximum absolute bolometric magnitude for 8 Dor of 
—5:47. Applying a bolometric correction of 0-17, again from Allen, we 
finally find the median absolute photographic magnitude to be — 4°52. 

(6) « Pavonis.—The reductions gave the mean minimum radius of « Pav as 
(12-5 +0°8)x 1o*km. Its radius at light maximum is then 14°5 x 1o°km. ‘The 
magnitude at light maximum is V = 3-95, and colour P— V = +0°35. Jacobsen’s 
estimate of the spectral type at light maximum is Fs, agreeing with the relation 
between spectral type and period found by Code and given in Aller’s book (16). 
We adopt therefore F5Ib as the spectral type, which gives an effective temperature 
of 6200deg. K from Keenan and Morgan’s table. The resulting maximum 
absolute bolometric magnitude of « Pav is —2-28. Applying a bolometric 
correction of o™-08 from Allen’s tables, we find a median absolute photographic 
magnitude of — 1-63. 

(c) 1 Carinae.—The mean minimum radius of / Carinae was found to be 
(84+6)x 10%km. Its radius at light maximum is then 87 x 10°km at which 
phase the magnitude is 3:26 and the colour +0-87 on the P, V system. The 
spectral type at maximum is given as F'8 by Jacobsen. ‘The spectral type at 
maximum corresponding to a period of 35:5 days is given by Code’s relation 
also as F8. ‘The effective temperature for a spectral type of F8Ib is 5300 deg. F 
and the median absolute photographic magnitude of / Car is then — 4:03. 

In order to compare the absolute magnitudes derived above with those of 
other cepheids, we must allot the three stars to either of Baade’s two population 
types. The criteria of population type membership can be divided into two main 
kinds, (a) spatial, i.e. whether the star is situated in or apart from the outer disk 
of the galaxy, whether or not it is associated with systems such as globular clusters, 
the galactic nucleus, spiral arms or galactic clusters and (6) kinematic, dependent 
on the velocity of the star with respect to the local standard of rest in the solar 
neighbourhood. ‘These criteria are correlated in some instances with the intrinsic 
properties of the stars themselves, e.g. the blue supergiants of the young Popula- 
tion | or the metallic line deficiences of some high velocity Population 1 stars. 
The cepheids form two parallel groups belonging to each population type and 
the differences between them has been well summarized by Aller (16). However 
as we shall see, these purely spatial and kinematic criteria may be ambiguous. 
B Doradus has a small systemic velocity relative to the Sun; but it is at a 
distance of about 200 parsecs from the central plane of the galaxy, a distance 
large for the general class of type 1 cepheids. Its period of 9:8 days places it 
at a minimum of the period distribution function for classical cepheids and this 
may be noted when consideration is given to the star’s peculiarities. ‘The 
period is remarkably constant, this constancy being a feature of the classical 
cepheids as opposed to the type 11 cepheids, which exhibit much more marked 
irregular fluctuations in the length of the period of oscillation. In view of the 
above we will regard 8 Dor as an anomalous type 1 cepheid. Eggen (17) has 
asserted that « Pavonis is a cepheid of Population type 11 because of the remarkable 
similarity of its light and colour curves to those of AL Virginis, a star definitely 
belonging to the high velocity group. « Pav has a systemic velocity of + 38 km/sec 











go A. W. Rodgers, Radius variation and Vol. 117 


and a distance from the plane of the galaxy of approximately 50 parsecs. ‘The 
spatial and kinematic criteria are again ambiguous. However, on the basis of 
Eggen’s photometry and because « Pav shows the variations in period more 
commonly observed in the Population type 11 variables, we will take it to be a 
type 11 cepheid. / Carinae is at a small distance from the galactic plane, has a 
small systemic velocity with respect to the Sun and can be regarded as a classical 
cepheid. 

Baade (18) has recently discussed the history of the derivation of the form 
of the period—luminosity relation. He had shown the need for, and others, 
including Blaauw and Morgan (19), have since given the value of, the revision of 
the original Shapley period- luminosity relation derived from the Magellanic 
cepheids for slope and the cluster variables for zero-point. ‘The period luminosity 
relation for classical cepheids with the Blaauw- Morgan revision of zero point ts 
M,,,= —1"-4(+0°3)—1™-74log P where M,, is the median absolute photo- 
graphic magnitude and P is the period in days. Again, if we take the apparent 
magnitude period relation for the Magellanic cepheids from Shapley’s work as 
m,,, = 17°04 — 1°74 log P, then from the modulus of the Magellanic clouds, given 
by de Vaucouleurs (20) from a summary of all determinations of 

m— M = 18-75 +0°10, 
we find the period luminosity law to be M_,,,= —1-70—1-74log P. ‘The agree- 
ment of this determination with the Blaauw- Morgan zero-point revision within 
the errors of the individual determinations affords some evidence for the intrinsic 
similarity of galactic type 1 and Magellanic cepheids. It should, however, be 
kept in mind that there are notable differences both in the amplitudes of the 
light and colour curves and in the frequency distribution of their periods. 

The zero-point of the period-luminosity law for the type ul population is 
fixed by the cluster type variables. After Baade’s suggestion of the need for 
revision of the zero-point of the type 1 period-luminosity law it was thought 
that the type 11 cepheids formed a continuous and single valued sequence of 
luminosity with period with a slope parallel to that found for the classical cepheids. 
Ina study of variables in several globular clusters, Arp (21) proposed the existence 
of a multivalued period-luminosity relation for the type 11 cepheids, each sequence 
being displaced from the other by a factor of two in the period. ‘The slope of 
each sequence is the same as for the type | relation. Recently Reddish (22) 
has suggested that, from consideration of the correlations between period, lumino- 
sity and radii, no such multivalued relationship exists. His form of the period 
luminosity relationship has then a significantly greater slope than that found 
for the Magellanic cepheids. Further photometric observations of type 
cepheids are required to settle this question. In Fig. 2 the absolute magnitudes 
of the cepheids for which the radii have been determined by Wesselink’s method 
are plotted against the logarithm of the period. ‘The two straight lines drawn 
are represented by the relations 


M,,=0"0-—1™74logP and M,,=—1™-6—1™74 log P. 
The stars 8 Dor and Y Oph lie 1-4 magnitudes from the brighter sequence on 
which fall the type 1 stars. 
In order to determine the distances of the stars x Pav, 8 Dor and / Car we 
must obtain the amount of interstellar extinction affecting the apparent magnitude 
‘The extinction is determined from the colour excess, the ratio of 


of each star. 
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absorption to reddening being assumed known. Knowledge of the colour 
excesses presupposes a knowledge of the intrinsic colours of these variables. 
‘T'wo recent determinations of the intrinsic colours of cepheids have been made. 
Stibbs (23), from considerations of a model of the galactic absorbing medium 
based on the colour excesses of B stars, gives the relation 

(P-V)max = +0™-17 +0™-18 log P. 
From observations of the colours of Magellanic cepheids Gascoigne (private 
communication) has suggested that the colours at maximum of the type I cepheids 
may be represented by the relation (P—V)n, max = +0™-10+0™13log P. ‘This 
equation represents the line which is the blue envelope of the observed colours 
on the colour-period diagram and thus neglects interstellar reddening within the 
galaxy in the direction of the clouds, but takes into account the possible reddening 
of the Magellanic cepheids by material in the Clouds themselves. In this 
paper we will take Gascoigne’s relation to represent the intrinsic colours of 


type 1 cepheids. 
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Pic. 2.—-Median absolute photographic magnitude plotted against log pertod. 


In determining the colour excess of 8 Dor we find another abnormality of 
this star. From Gascoigne’s relation, at maximum light, the intrinsic colour of 
B Doradus is +0™-23 on the (P—V), system. ‘lhe colour excess is then o™-23 
and the photographic absorption is one magnitude. However, de Vaucouleurs 
(24) has shown that £ Dor is seen projected against the edge of a spiral arm and 
thus would appear to be in a relatively transparent region. Also it would be 
expected that any galactic reddening suffered by 8 Dor would equally attect the 
Magellanic cepheids ; thus it would appear that the intrinsic colour is anomalous. 


o™-23 the distance is 295 parsecs. 
Assuming type I and type 1 colours to be similar, « Pav has an intrinsic colour 
at maximum of +o0"-35 and the colour excess is then o™-14 leading to a photo- 
graphic absorption of o-6 magnitudes. The resulting distance is found “ 
116 parsecs. y 
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/ Carinae has from Gascoigne’s relation an intrinsic colour of +0™-30. The 
star is reddened by 0°47 magnitudes and the distance is found to be 179 parsecs. 

4. Radius variations and population type-—Using methods similar to those 
described above, several workers have determined radii of other cepheids. As 
the method demands observations pushed to the limit of their accuracy it is 
essential to have available good velocity curves and first-class photoelectric 
observations at or close to the same epoch. In Table III are listed the stars 
whose radii have been determined, the maximum and minimum values thereof 
and references. 


Tasie III 
Cepheid variables with known maximum and minimum radu 
Star Period Population Rnin Rusax Reference 
Type < 10° km 10° km 
d 
RR Lyr 0°57 i 4°77 5°23 (2) 
3 Cep 5 I 3571 38°9 (1), (2), (8) 
(9) 
7 Aqu 7°18 I 44°6 49°4 (2), (8), (9) 
Kk Pav 9°07 II 12°5 s9°3 
B Dor 9°84 I 69°8 75°5 
Y Oph 17‘I I 47°! 52°5 (3) 
W Vir 7°3 Il 18 40 (4) 
/ Car 35°5 I 84 108 


In Stebbins’ paper (2) the mean radii only were given. ‘lhe radii extremes 
were taken from his Fig. 1. Stebbins’ radii gave confirmation of the new zero 
point of the type I period-luminosity curve. 

Abt has investigated the cepheids Y Ophiuchi and W Virginis. Y Oph he 
found to be an anomalous type I cepheid which showed the opposite effect to 8B Dor 
of being fainter than predicted by the period—luminosity relation. His compre- 
hensive discussion of W Virginis led him to the view that this star “expands 
periodically from about 18 to 40x 10%km. At each expansion it throws off a 
shell which continues to rise for 35 x 10°km...” W Virginis is definitely a 
type 11 star and its behaviour has been regarded as typical of the type 11 variables 
of similar period. 

From the data of ‘Table II] we can form the quantity 

5R/R=2 x (Rmax — Rmin)/(Rmax + Rmin) 
which we shall call the relative radius amplitude. 


Fig. 3 shows the plot of relative radius amplitude against the logarithm of 
the period for the stars listed in Table II. 


‘The stars regarded as classical cepheids show quite a slow variation of 5R/R 
with period. The type 11 variables however show a rapidly increasing 5R/R 
with period. An implication of Fig. 3 is that the mechanism of pulsation, and 
consequently the structure of the star, is quite different in the two population 
types. Of note is the convergence of the two curves for the cluster type variables. 
Population type criteria as applied to a single star have been in the nature of 
spatial or kinematic correlations with the colour luminosity arrays of a given 
stellar system. We see however that there are, at least for intrinsic variables 
in the period range considered, structural differences for an individual star 
dependent on the population type but largely independent of its present spatial 
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or kinematic situation. his is clearly portrayed by those cepheids having 
periods between 2 and 10 days, such as « Pav, whose light curves do not show the 
characteristic hump of the W Virginis stars and which have no distinguishing 
spatial or kinematic property. 
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Fic. 4.—Light difference at phases of mean colour plotted against log period. 
Determinations of absolute radius by Wesselink’s method are limited to very 
few stars due to the stringent observational criteria mentioned above. However, 


a measure of 5R/R can be obtained purely from the photometric data. If for 
each star listed in the paper of Eggen and Gascoigne (6) we find the difference 
in brightness from the light curve of the two phases where the star’s colour 
equals its median colour, we are essentially determining the part played by the 
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variation in radius in producing the observed light curve. This difference we 
denote by (65M)z and plot against the logarithm of the period in Fig. 4. 

Fig. 4 shows, with more scatter, the same features as Fig. 3. ‘The relative 
radius amplitude is of course a function of the form of the velocity curve as well 
as of (5M)z and the increased scatter of Fig. 4 is to be expected. 

AL Vir falls near « Pav on the type 11 curve and the star RT Aurigae, which 
according to Eggen and Gascoigne has.a large light amplitude for its period of 
3°7 days, also falls near the type I! curve. 

Arp (21), has obtained the photographic and photovisual light curves of 
cepheids in globular clusters with periods greater than one day. Using the data 
given in his paper, (8M)z was computed for his stars of period greater than five 
days. ‘These stars are plotted along with the other type 1 variables as open 
circles in Fig. 4. 

From Arp’s stars plotted in Fig. 4 the following is evident. At the point 
corresponding to the period of the longest period stars showing the hump-backed 
light curve of the W Virginis type, (6M)g ceases to increase. At this period of 
about twenty-five days, 5R/R is approximately unity, i.e. Rmax/Rmin=3, so the 
possibility arises that the amplitude of oscillation becomes so large that stars of 
longer period having homologous structure are unstable. 
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VARIATIONS IN THE LIGHT OF SIGMA SCORPII 


A. R. Hogg 
(Received 1956 October 19) 


Summary 


Photoelectric light curves of the 8 Canis Majoris type variable, Sigma 
Scorpii, were obtained nearly simultaneously with a series of radial velocity 
measures made by Struve and collaborators. The photoelectric results 
show the existence of a beat period of 84-028 superposed on the known 
short period variation of 04-246846. 

The star is shown to be brightest and bluest when at or near minimum 
size. 





Introduction.—Sigma Scorpii forms ADS 10009 with an 8™-7 common proper 
motion companion 20” distant (1). The bright star, Br III (2), is both a single 
lined spectroscopic binary with a period near 34 days (3) and a 8 Canis Majoris 
type variable with a period near 04-25. 

Dr O. Struve suggested that it would be of interest to re-investigate the light 
variations of this star (which had previously been observed at Mt Stromlo (4)) 
simultaneously with radial velocity measures to be carried out by him and 
Dr McNamara at Mt Wilson and Mt Hamilton in 1954. (The star is badly 
placed for precise photometry in California but passes near the Mt Stromlo 
zenith.) ‘The measures could not be made exactly simultaneously because of 
the longitude difference at the two observatories. Photometric observations were 
made in 1955 as well as 1954. 

Observational.—The measures were made in three colours at the Cassegrain 
focus of the 20-inch reflector at Mount Stromlo using an unrefrigerated EMI 
Type 5659 11-stage electron multiplier connected to a balanced d.c. amplifier 
and Brown recorder. Three filters were employed, viz. : 

U, Corning 9863 (standard thickness) ; 
B, Chance OB: (1 mm)+ Schott GGr5 (2 mm); 
V, Chance OY4 (2 mm). 

The comparison star, 7 Sco, had been used in the earlier Stromlo work and 
was found to be satisfactory. Although of nearly the same spectral type as o Sco 
and near to it in the sky, it was distinctly bluer than the variable. Because of the 
relatively narrow range of each of the filters employed the effect of the difference 
in star colours on absorption was neglected. Small absorption corrections, 
generally less than o™-o1, were applied, using mean values of the absorption 
coefficient found for other work with the filters. Measures were made through 
a diaphragm subtending a diameter of 27” which was small enough to exclude 
the light from the 1oth magnitude companion of the variable. Observations of 
the comparison and variable were alternated and magnitude differences between 
the reading for the variable and the mean of the immediately preceding and 
succeeding readings of the companion star were computed. 
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Six substantially complete light curves (i.e. curves covering the whole of 
one cycle with readings of the variable at average intervals of 12 minutes) were 
obtained during 1954 and a similar complete curve was obtained in 1955. 
Fragmentary curves designed to give colour at maximum were also obtained. The 
“‘complete”’ curves are shown in Figs. 1 a, 1 6 and 1c where the points are current 
means of three observations. 
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Fic. 1 6.— Short period light variation of « Sco-B filter. 
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Fic. 1 c.—Short period light variation of a re V filter. 


Periodicities —The system is known to exhibit three periodicities, viz. : 

(a) a periodicity of 34 days in the value of the radial velocity of the centre 
of mass ; 

(5) a short period variation of o“-25 in the radial velocity and also in the 
light output; and 

(c) a period of 8 days in the amplitude of the radial velocity curve recently 
discovered by Struve et al. and probably the result of a beat between two 
interfering periods close to 0-25. 

The current observations are not sufficiently numerous to detect any 34-day 
periodicity in the light. The o4-25 period is however clearly shown ‘in both 
light and radial velocity curves. Observed heliocentric times of blue minima 
are shown in Table I together with their differences (O—C) from the values 
derived from the following expression, viz. : 

Minimum light = J.D. (Hel.) 2 432 306-795 + 0"-24 684E. 


TaBLe I 
Observed and computed minima for a Scorpii (Blue Light) 


Minimum J.D. (Hel.) 2 432 306°795 + 04-246 846E 


: Observed ‘ Observed 
).D. (Hel.) sian Fhese® }.D. (Hel) O-C _Phase* 
i d d 
: 2432308 °038 0006 68 2434903121 0003 °'9 
é 2313°956 0°000 4°7 4904°106 0:006 1°8 
: 2350°007 ~O'Oll o°6 4906°062 0°024 38 
2365,040 0'O14 76 4911°026 0003 o'7 
pi 3396°107 0°022 31 4912°979 0019 2°7 
3397°106 O'oll 4'1 5265-982 0006 2°5 


*Computed from T’=J.D. 2 434 878-2 + 84-028E. 
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The standard deviation of a single O-C value in Table I is +04o11 but, as 
will be discussed later, systematic variations occur in the O-C values. 
A. D. Williams has privately communicated a result J.D. =2 434 878-92 which 
agrees exactly with a value computed from the above expression. The period 
of 04-246846 agrees very well with the figure of 04-246844 derived from the 
spectroscopic results. 

The results shown in ‘Table I refer to blue light and were obtained with three 
different photometers (4). Although the time of occurrence of minimum does 
depend slightly on the wave-length of the light in which it is observed, the 
different effective wave-lengths of the three blue photometers are sufficiently 
close to allow the effect to be neglected. 

The earlier results (4) gave a period of 04-246841, 1.e. 04000005 shorter 
than the present figure. Inspection of the O-C values in Table I does not 
suggest any certain systematic deceleration in the system, as the mean O-C for 
the first six entries is — 0-005 + 04-005 s.d. compared with — 04-009 + 04-004 s.d. 
for the second six. ‘The rather large standard deviation could be in part due to 
the small amplitude, rendering difficult the exact determination of the precise 
time of minimum. Approaching minimum the rate of change of blue light is 
about o™-o05 per o%-o1 so that a determination of minimum even to 04-005 
requires a precision of nearly o™-o02 in the photometry. This of course is at 
or near the present practical limit. 

The variation of the time of minimum of the light curves with the effective 
wave-length of the photometric system is shown by the results of a harmonic 
analysis for the curve obtained on 1954 June 17, viz.: 


U B V 
Minimum light J.D. (Hel.) 2 434 000 + gi 17016 QII°O17 gi 1-026 
Amplitude + 0™-053 +o0™-o41 +0™-026 


The results refer to the first harmonic and show that as the wave-length increases 
the amplitude decreases and the phase becomes later. 

The 8-day periodicity found by Struve et al. (3) in the values of 2K is also 
exhibited by other quantities connected with the star. The list of such periodicities 
is thus as follows: 

(a) a periodicity in the magnitude reached at maximum light in successive 
04-25 cycles (Ammax) ; 

(6) a periodicity in the difference ‘‘observed minus computed” time of 
minimum light (O-C), in the 04-25 cycle; 

(c) a periodicity in the difference ‘‘observed minus computed” time of 
maximum radial velocity (O-C)y of the 0-25 cycle; 

(d) a periodicity in 2K of the 04-25 cycle. 

The periodicity in Ammax is presented in the data collected in Table II 
where column 3 shows the phase computed from the expression 


T’ =J.D. (Hel.) 2 434 878-2 + 8%-028E. (2) 


‘The tabulated figures refer to magnitude differences (in the natural scale of the 
photometer) between o Sco and 7 Sco and are in the sense that numerically 
larger figures mean that o has become fainter. The results in the individual 
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Taste II 
J.D. (Hel.) 4m at max. light (o-r) Phase* 
2430000 +- U B V 
4903008 0°642 0°346 0°058 0°76 
04°244 0°625 0°333 0°050 1°99 
06°218 0°640 0°330 (0°057)t 3°97 
10°177 0°690 0°379 O13 7°93 
II‘1I50 0615 0°333 0°055 o’9g! 
12°138 0°635 0°332 0°065 1°87 
5266-120 0°634 0°344 0-069 2°58 
67°115 0°655 0°355 0°076 3°61 
69°085 0°672 0°368 0°077 5°53 
86-121 0°656 0°342 / 0°072 6°53 
g0°076 0°647 0°336 0°057 2°47 
95018 0°677 0°364 0080 7°41 


* Computed from 7” =J.D. 2 434 878-2 + 84-028E. 
+ Estimated from U, B and mean colour. 
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z Fic. 2.—Diagrams showing the variation of the maximum light attained in the individual 04-2468 
4 cycles of « Sco, Ammaz; and similarly the variation of the observed minus computed time of 
; maximum light, (O—C),; the observed minus computed time of radial velocity, (O—C)y; and 
g the amplitude 2K of the radialvelocity. These quantities found in individual 01-2468 cycles have 
been plotted against phases derived from the expression J.D. 2 434 878-2 +-8-028. 

In the 2K diagram the dots represent observations by Struve et al. and the crosses results 
obtained at Mount Stromlo (see text). 
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colours are rather erratic but the periodicity is shown better if a running triple 
mean for the three colours is plotted as in Fig. 2. 

‘The (O-C), results for blue minima are presented in ‘Table I, where column 3 
gives phases computed by expression (2) above. ‘These results, after smoothing 
by taking running means of three observations, are shown in Fig. 2. 

The (O-C)y results were obtained by scaling off the diagrams in the paper 
by Struve et al. (3). The computed times were derived from the expression 


T” =J.D. (Hel.) 2 434 861-980 + 04-246 846E. 


The periodicity in 2K is illustrated by what is mainly a reproduction of Fig. 3 
of Struve et al. but plotted with a period of 8*-028 in place of the figure of 84-o 
used by these authors. Three additional results obtained by courtesy of 
Drs Gollnow, Hagemann and Przybylski in the winter of 1955 have been 
included in the diagram and are indicated by crosses. ‘These new results are: 


J.D. 2K 
243528902 102 Km sec 

94°93 84 55 

97°16 97 » » 


The times refer to the occurrence of the maximum velocity and the ranges 
to the difference between it and an adjacent minimum. Actually the crosses 
would fall somewhat more closely amongst the dots if a period of 840g was 
adopted but with the use of the 84-028 period the three 1955 measures fit 
reasonably well amongst the more numerous 1954 observations and could well 
be part of the same series of measures. 

No one of the four curves shown above, with the possible exception of the 
(O-C)y result, presents by itself very convincing evidence for the 8-day 
periodicity, but together they indicate the reality of this variation. It is not 
improbable that more extended and precise work would show that the 8-day 
cycle itself does not repeat exactly. 

Change in dimensions.—Various workers, by combining light and velocity 
curves, have determined dimensions of pulsating systems. Because of the 
obvious changes from cycle to cycle in o Sco it would be most desirable to compare 
light and velocity curves taken simultaneously. With the present data it was 
not possible to do this because of the longitude difference between the observatories 
at which the bulk of the spectroscopic results and the photoelectric observations 
were obtained. Results in two adjoining cycles however were examined, 
viz. J.D. 2434910°657 to 2434910910 (spectroscopic) and 2434910°923 to 
2 434.911°161 (photoelectric). Comparison of the above radial velocity curve 
with that during J.D. 2434911740 to 2434911-g00 shows no appreciable 
difference, so it seems fair to assume that no great changes took place in the 
intervening period and that the above spectroscopic and photometric results 
refer to essentially the same condition of the system. 

The results are summarized in Fig. 3 which gives light, colour and radial 
velocity curves as well as a curve showing the displacement obtained by numerical 
integration of the velocity curve, assuming that limb darkening and projection 
effects (K = 0-6) are accounted for by the equation V = 24/17V, (where V = radial 
velocity to be used to obtain the displacement curve and V,=the observed radial 
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velocity change (§)). ‘The graph is arranged so that luminosity, temperature 
(represented by colour) and radius (indicated by the displacement curve) all 
increase upwards. For the particular cycles investigated the curves show that: 
(1) at or near minimum size the star is brightest and hottest; and 
(2) at or near maximum size the star is faintest and coolest. 
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biG. 3.—Showing changes in light (U, B and V), short wave and long wave colours (U-B) and 
(B-V), radial velocity RV, and displacement during the o'-246 846 day cycle of o Sco. 


Both of these statements require elaboration. In the first place only the 
B maximum coincided with minimum size. Maximum U light occurred about 
o-1 period (38 min) before this and maximum V light about 0-05 period after 
minimum size. The observed U, B and V minima were not definitely separated 
in time although the V minimum may have slightly preceded the U and B minima, 
and all minima may have been about 0-02 period in advance of the maximum 
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size. Statements about the temperature are based on the U-B and B-V 
colour differences. 

To the above should be added the statements relating to changes within the 
8-day period and displayed in Fig. 2, viz.: 

(3) the brighter maxima occur somewhat before the predicted time and the 
fainter maxima occur somewhat later than the predicted time ; 

(4) the brighter maxima are associated with a smaller range of variation (2A) 
of velocity in the radial velocity curve, i.e. are associated with less displacement ; 

(5) the bigger ranges of velocity tend to occur a little later than predicted. 

The data given in Fig. 3 should lead to an estimate for the radius of the star 
using well-known methods as exemplified in the work of Stebbins (6) on cepheid 
variables. If pairs of phases of P and P! be selected so that the colours at these 
points are respectively equal (i.e. (U—B)=(U-B)' and (B-V)=(B-V)') 
then it may be assumed that the surface temperature of the star is the same at 
the two phases and that any change in light is due solely to a change in dimensions. 
The fractional change in light is known (from the light curve) and the absolute 
change in radius can be found by integration of the radial velocity curve. ‘These 
two quantities then lead to a figure for the radius. In the cycle discussed only 
one pair of phases of the required nature could be located. ‘hese occur at 
04-082 and 04-229 in Fig. 3. Inspection of the figure at these two points shows 
that an increase of radius of 140000 km at a nominally constant temperature 
is accompanied by an apparent 1 per cent increase of luminosity. Whether 
this result is to be attributed to observational inaccuracy, the use of photometric 
and spectroscopic results covering different cycles, or to deficiencies in theory, 
is an open question, but it is of interest to note that Walker (7) in a similar 
investigation of 16 Lac also found luminosity and radius to vary in a manner 
opposite to that to be expected. 

Colour and magnitude.—Except for Fig. 3 the foregoing results have been 
given in the “natural system” of the photometer. This corresponds closely 
enough in the V magnitude to the V system of Johnson and Morgan but V and 
B natural systems differ somewhat from the standard U and B systems. 

Comparisons of o Sco and 7 Sco against nine Johnson and Morgan standard 
stars on three nights gave the following results reduced to the Johnson and 
Morgan system : 


V B-V UB 
ao Sco (mean) 2°92 -or'16 0°64 
7 Sco 2°84 —0'27 —o'98 


Using the relations given by Johnson and Morgan it is found that o Sco is 

distinctly reddened, the intrinsic colours being 
(B-V),;=—0-26 and (U-B),= —0-94. 

Conclusions.—'There is nothing in the present observations to cause o Sco 
to be regarded as other than a typical 8 Canis Majoris star. The intrinsic colour 
fits reasonably well into the colour—period relation obtained by McNamara and 
Williams (8). The results show the same phase relations between radial velocity 
and light curves that have come to be associated with BC Ma stars and thus 
confirm the work of Williams and Struve (9). Although it has not been possible 
to obtain a quantitative confirmation of the pulsating nature of the system, the 
light, colour and displacement curves are in agreement with the idea of pu!sation 
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in that the star is bluest and brightest at minimum size. Also beat periods in 
light curves have been observed in the case of 16 Lac, 12 Lac and v Eri (de Jager 
(10) and Walker (7)) and may exist in 8C Ma and BW Vulpecula which show 
changes in their light variation. 
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from both Mrs A. de Vaucouleurs and Mrs C. Wehner in the reduction of the 
observations and also to thank Drs Gollnow, Hagemann and Przybylski for 
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Summary 


Jeans’ formula gives the condition that a gravitating mass of gas shall be 
stable to small fluctuations in the density. It was originally obtained for a 
static mass of uniform density, and the proof is questionable. In this paper 
the formula is derived by studying small perturbations of the density and 
velocity in an expanding Newtonian world-model. It turns out that, given a 
suitable interpretation of the quantities involved, the formula is valid in all 
cosmologies except the steady-state theory. 

Although instability occurs in the models as suggested by jeans’ formula, 
the condensation process is too slow to account for the formation of the 
nebulae—a result which has now been reached in several ways, and which may 
be regarded as well-established. 

In the Newtonian version of the steady-state theory small perturbations of 
the average density die out, and so Jeans’ formula does not apply. The 
relevance of this to the work of Sciama, who assumes the contrary, is discussed. 
It seems that the claim that the steady-state theory explains the formation of 
the nebulae requires further substantiation. 





1. Introduction.—Jeans’ formula (4, 5) gives the maximum size consistent 
with stability of a cloud of gas whose density p, and equation of state p(p), are 
known. It states that a cloud with linear dimensions greater than 


= dp 1/2 
(=F) (1.1) 


(where G is the constant of gravitation) is unstable to slight fluctuations in the 
density. ‘This formula is often used in astrophysical work, and important 
deductions are made from it. 

In (1.1) the density and pressure are supposed to be uniform throughout the 
gas. In fact, however, if gravitation is ta!.cii into account, the equation of hydro- 
static equilibrium has no solution for a tinite uniform mass. This is disturbing, 
and one’s uneasiness is enhanced if one remembers that certain of the polytropic 
gas spheres (which do satisfy the hydrostatic equation) are mechanically stable, 
so far as is known, whatever their size. On the other hand, as I have recently 
shown (3), an isothermal gas sphere of sufficient size is unstable, and the condition 
for instability is similar to (1.1). 

Since Jeans’ work appeared, uniform distributions of matter have been 
studied in great detail in cosmology. For the most part, this work has been 
based on general relativity, in which uniform distributions are quite unobjection- 
able as there are rigorous solutions of the field equations referring to them. 
However, the Newtonian cosmological models of Milne and McCrea (9, 10) also 
use uniform distributions of matter, and these seem to be a natural vehicle for 
a further study of Jeans’ type of gravitational instability. As a result of recent 
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work of McCrea (8), the assumptions on which these models are based are now 
well understood, and they can be used with confidence as a guide to the models 
of general relativity, and also of the steady-state theory. 

My intention in this paper is to use the Newtonian models to study the 
gravitational instability of a uniform material filling an expanding universe. I 
write with full awareness of the extremely beautiful and powerful work of Lifshitz 
(6) on gravitational instability in relativistic cosmology, and my results, as far as 
they are comparable, agree with his. However, the extra simplicity gained by 
using the Newtonian models, apart from giving formulae which are more tractable 
and readily comprehensible, enables one to extend the investigation to a wider 
class of models—namely, those with cosmological constant, and that representing 
the steady-state universe. 

The result of my work, very briefly, is that something very like Jeans’ 
instability occurs in all the cosmologies, except the steady-state theory. ‘This latter 
conclusion is unfortunate because the steady-state theory is the only one which 
has at present an explanation of the formation of the nebulae (11), and it now 
seems that this explanation is made doubtful by an unjustified use of Jeans’ 
formula. However, it may be possible to retrieve the situation, as will be ex- 
plained later. 

The plan of the paper is as follows: in Section 2, | obtain Jeans’ formula in its 
original form, and I follow this in Section 3 by a derivation of the corresponding 

results in the expanding universe ; this is followed in Sections 4-6 by a discussion 
of the results for certain special models, and in Section 7 by the application of 
the theory to the steady-state universe; in the Conclusion, Section 8, I give a 
brief summary and general discussion of the results of the earlier sections. 

Jeans’ formula.—'Vo derive Jeans’ formula we use the theory of sound in a 
large mass of gas, taking gravitation into account. Let us suppose that the gas 
satisfies an equation of state p _g(p), and is at rest except for a small velocity 
u(r, t) due to the wave motion. Let p(r, ft) and p(r,t) be the density and pres- 
sure; the space derivatives of p and p are not at the moment supposed to be small, 
but we shall assume that the deviations of p and p from their equilibrium values 
are of the same order of smallness as u. Squares and products of small quantities 
(and of their derivatives) are to be neglected. 

On these assumptions the hydrodynamical equations 


+4 (u. V)u=F-~ vp, (2.1) 
Fy .(p)=0, (2.2) 
become = ~ = : Vp, (2.3) 
4 V .(pu)=o0, (2-4) 


where F is the body force per unit mass, consisting solely of gravitation. Since 
the gas is in equilibrium except for the small perturbations, we have 


I 
F,= Po V Po (2.5) 


where a suffix o refers to the equilibrium value. Put 
F=F,+F,, P=Poto, 
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where F, and w are small; then Poisson’s equation gives 


ur F, =— 47Gpo, Vi F, = 47Gw. (2.6) 

Equations (2.3) and (2.5) give 
eu | w dp " 
‘acme 3 on 


where the equilibrium value of dp/dp is to be taken. From (2.6) and (2.7) we 
have 


c wa 7} 
=v ~u=V B- VET 








py dp 
fw d 
_ 47Gw _ 72 = dp\ : (2.3) 
po ap} 
Equation (2.4) may be written, after differentiation and neglect of second order 
terms, 
Cw Cc cu 
0} + py = (V.Uu)+(V po). = =0O. (2.9) 
ct ct ct 
From (2.8) and (2.9) we find 
as V cu 1 
no + (V po) ~~ =47Gspyo+ V7 oP ° (2.10) 
ct Po ct dp 
where we have introduced the condensation s defined by 
S=W Po- 
If now we suppose that py, is constant, (2.10) reduces to the formula given by 
Jeans’ 
cs 5 eae dp 
ro = 47Gpys + ; *(: 7) ‘ (2.11) 


Jeans (4) is not quite explicit on this point, and his method is somewhat different, 
but careful study of his work shows that in order to get (2.11) he also is required 
to take py as constant. If p, (and therefore py) is constant, (2.5) gives F, zero 
throughout the gas, and one is at first inclined to think that the gravitation has 
been removed from the problem. However, one can retrieve the situation, 
after a fashion, by supposing that one has an infinite mass of gas so that, as there 
can be no preferred direction for F, it must be zero. (Strictly, Fy is in this case 
not zero, but undefined.) As, however, Jeans goes on to apply the investigation 
to a finite mass of gas, the argument is, to say the least, precarious. 

The deduction of Jeans’ formula from (2.11) is simple. Consider a wave of 
length A propagated along Ox and suppose a solution of the form 


s=h(t) cos 27x/A. 
Since dp/dp is constant, we find 





d*h ; 47" dp 
r TE = | +769 a dp h. (2.12) 
Hence the disturbance increases exponentially with time if 
Ci 7) 
A> =}. 2.1 
(e dp (2.13) 


A mass of gas with linear dimensions greater than this will allow such disturb- 
ances, and will therefore be unstable. 
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As it will be convenient in the following sections to study spherical waves, | 
shall show that (2.13) holds also for such waves. If we suppose that the disturb- 
ance has spherical symmetry, we may write 

s=A(t)r™ sin 2zr/A. 
Equation (2.11) then gives (2.12), and (2.13) follows as before. 

3. Gravitational stability of Newtonian world-models.—Let us consider smail 
perturbations in a Newtonian world-model. Suppose that a definite observer O 
sees the motion of the world-fluid as purely radial, and that the velocity, density 
and pressure have spherical symmetry about O. ‘The imposition of this sym- 
metry on the perturbations is not essential to the analysis, which can be followed 
through in a similar way if small non-radial motion and deviations from spherical 
symmetry are allowed. ‘These latter perturbations make little difference to the 
final result, and are not important from the point of view of stability, so for 
simplicity I confine the argument to the use of radial perturbations with spherical 
symmetry. 

The hydrodynamical equations (2.1) and (2.2) give 


Cu Cu . ep 
= +ue I =—, (3.1) 
ct cr por ' 
Cp op Ou = 2u 
~~ Se = +e SS 7 THs, (3.2) 
ct cr cr r : 


where wu is the radial velocity measured by O, and /p, p and F are the pressure, 
density and radial force per unit mass. As the motion is supposed to be that of 
a slightly perturbed Newtonian world-model, we may write 

u=nrf(t)+ (7, t), } 

! 

P= Ppolt) + (r, ¢), ; ( 

P=Polt)+ (7, 4), 
where 7, w and g are small functions whose squares and products (and those of 
their derivatives) may be neglected, and where f(t), po(t), Po(t) are the functions 
which occur in the unperturbed model. ‘The body force per unit mass F 
satisfies 


oy) 
ws) 
— 


># 
V.F- — + — =—4nGp+a, (3.4) 


where A is the cosmological constant. In the homogeneous model 7, w and q 
are zero, and the functions p, and f satisfy 


Po — = 3h 3f+ 3f° ? 47Gpy T A, (3-5) 
where * denotes 0/8t. Later we shall need the function R(t) given in terms of f by 
f=RIR. (3.6) 


From (3.5) and (3.6) we find 

ep =aR-, 
where is an arbitrary constant. Equations (3.5) follow from (3.1)-(3-4) 
if 7, w and gare zero. ‘lhe pressure does not appear in the theory of the homo- 
geneous model since it enters (3.1) only through its gradient. We shall suppose 
that there is an equation of state of the form 


Pp a(p) 


at hz lg 
fae: yr, th=alr, OF. (3-7) 





108 W’. B. Bonnor Vol. 117 


Here and elsewhere it is to be understood that dg/dp is to be calculated at p= py. 
If we substitute (3.3) and (3.7) into (3.1) and (3.2), and neglect small quantities 
of the second order, we find 
Bis Olas . pp w dg 
rf+o+nf*+rfe'+fo=F- —=, 8 
f f?+rfe' +f i (3-8) 
py + @ +1fe’ + 3pof + 2povr + por’ + 3uf =0, (3-9) 
where ‘ means ¢/er. If from (3.8) we form 0F/dr + 2F/r, and use (3.4) and (3.5), 
we have 


t 2fv 1d ain. ae 
t' +rfo" + 4fe’ + 2- + afv = —47Gw — E(u + =) ' (3.10) 


where once again products of small quantities have been neglected. Equations 
(3.9) and (3.5) together give 
a +rfw' + 2pevr! + pv’ + 3wf =0, 


which may be written 





C - é 
> (r?2v) = — —(w+rfw' + 3wf) = X, say. (3.11} 
0 
Equation (3.10) may be written 
a f & P 1 dg ” 2w’ 
aumbitth+ + ost ha ~saGe~— 2 3-4 t 12 
r Ara | " r or? (r'v) —e Po dp \* hee 5-88) 


then (3.11) and (3.12) together give 


: y2 d: 2 py’ 
X + rfX’ —— 47Guwr* _- a % (w" + ; ) 
p a 


} 





Introducing for X its value in (3.11), we find, after a rather long calculation in 
which use is made of (3.5), 
£2 dg ~ 2u)’ 
& + 2rfo' + r*f*w" + 8fw + rw’ (9f? +f) + w(18f? + 6f— jo a we + ; (3.13) 
p 3-13 


r 





This is a differential equation for the perturbation in the density, w, in terms of 
the functions f and dg/dp, supposed given. 
To simplify (3.13) we change to a new independent variable x by the trans- 
formation 
r=xR(t) 
where R is the function introduced in (3.6). A long calculation shows that 
(3.13) becomes 


Ow 2 0w dg o-w R R 
553° > ie -(¢ a |S sy Fy to(12% +68 -A)]. (3-14) 


Finally, we introduce the condensation s a 
S=w/p,= «a wR’, 
which reduces (3.14) to 


Os 20s _(de\1 [aR . R 
73 + 5 -(F) re 4 2 s+0(3R-A)]. (3-15) 


Equation (3.14) or (3.15) is the equation representing the propagation of 


spherically symmetric sound waves through the Newtonian world-model. Let 
us study solutions of the form 
sin — 


= h(t), (3.16) 
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where & is a real constant and where A(t) is a function to be determined which 
will show how the amplitude of the wave changes with time. Substituting 
(3.16) into (3.15) we find 


shi (3h op2 48 
hvagh+n(S —A+R*R *) =O, 


which with the help of (3.5) and (3.6) gives 


. u 
h+aRh+h(R*E — 4nGp) =o (3-17) 
In (3.16) sin kx/kx is a bounded oscillating function of r with order of magni- 
tude 1 near the origin for all time, so the growth of the condensation is measured 
by A(t). ‘Thus A(t) is strictly comparable with the function h in Section 2, and 
formula (3.17) corresponds to equation (2.12). ‘The similarity between the 
two equations is obviously close, but in (3.17) there is a secular change in the 
coefficients owing to the expansion. 
It can be shown from (3.17) that condensations can begin to form in favourable 
circumstances. Suppose that there is a perturbation at time ¢, such that 
dh 


h=h,>o0, > temas o at t=t,, (3.18) 


where h,, h, are assumed given; suppose also that 


k2R 3 —4nGpy<0o for t)<t<ty. (3.19) 
P 
Immediately after ¢,, his positive so that h increases; and A will stop increasing 
only if it reaches a maximum, for which the conditions are 
h=o, h-<o. 

Now from (3.17) and (3.19) it follows that, when 4=o0 and h-- 0, h>o so that h 
does not achieve a maximum in ¢, <t ty, and it must go on increasing. Hence, 
given the initial conditions (3.18) and also the condition (3.19), a condensation 
begins to form. 

We may express (3.19) in a form closer to Jeans’ formula if we introduce the 
wave-length A of the disturbance by 





A=27R/k. (3.20) 
Substituting (3.20) into (3.19) we find 
2: 7.@ 
Gp ap” ii 


Here A is a function of the time as expressed by (3.20), and so are py and dg/dp. 
The condensation process continues so long as (3.21) is satisfied. Equation 
(3.21) is formally the same as (2.13), 80, with the appropriate interpretation of the 
quantities involved, Jeans’ formula is true in the expanding universe. In terms of 
the initial wave-length A,, (3.21) may be written 
a dg R? 


A2>=— s+ >> 
fuss Gp, dp R* 


(3-22) 


where R, is the value of R at time ¢,. Supposing this inequality satisfied at ¢=1,, 
the condensation process will continue until such time f, as it ceases to hold. 
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Let us consider the special case in which the equation of state is 


p = g(p)=Kp*”. (3-23) 

‘The formulae (3.21) and (3.22) then simplify, and give as condition for instability 
2 47K 

A,? (3-24) 


; 3Gp,23’ 

2 1 
where p, means the unperturbed density at time ¢,; if this condition is satisfied 
at t,, condensation will proceed indefinitely. Equation (3.24) ts precisely Jeans’ 
condition for a gas with equation of state (3.23). ‘The equation of state for diatomic 
hydrogen (p=xp"*) is not very different from (3.23), and in Section 5 I shall use 
the latter equation to obtain an estimate of the behaviour of an expanding mode! 
filled with hydrogen. 

If we had not confined attention to spherically symmetric perturbations w« 
should have found an equation similar to (3.15) but containing angular terms on 
the left-hand side, making up the full Laplacian. ‘The solution (3.16) would 
have contained a surface harmonic, and a spherical Bessel function would have 
replaced sin kx/kx. The results of this and the subsequent sections would not 
have been basically different. 

The foregoing results will now be applied to certain special models. 

4. The Einstein universe.—In the Einstein universe we have 

R=const. A=47Gpp, 
so that (3.17) becomes 


. 2 do 
h+h (f ~~ +7Gpn) = ©, 


A® dp 
where A=27R/k. 
Thus h depends exponentially on the time if 
_ (_7_ dg\"”. 
(ea) 


This is exactly the same as (2.13) and gives the maximum size of a stable sphere. 
It follows that feans’ formula is true without modification for the Einstein universe. 
It should be noted, however, that the constant density of the unstable gas cloud 
must be that appropriate to the Einstein universe, i.e. A/4mG. The theory 
would not apply to a uniform mass of gas of different density. 

5. The expanding world-models with A=o.—The main interest in the present 
theory is in what it predicts about the formation of the nebulae in an initially 
homogeneous universe. Except for the steady-state theory, which will be dealt 
with later, most of the plausible models start from a singular state of infinite 
density. For some time after the singular state the pressure would be large, 
and for this period the Newtonian models do not give a satisfactory representation 
because Newtonian theory does not allow for the gravitational effects of the pres- 
sure. ‘The theory constructed here will therefore apply only to the later stages 
in the history of the universe when the pressure is small. From Lifshitz’s 
work it seems that the later period is the most favourable for the formation of 
the nebulae, so the inapplicability of the theory to the earlier stages is not very 
serious. 

Equations (3.5) and (3.6) for the unperturbed model may be integrated to 


give , x - 87Ga 


+ ~AR% +e, (5-1) 
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where « and ¢ are constants of integration. Let us now consider a definite 
model, obtained by taking 
A=o0, e= — B?, 


where f is a real constant. ‘This model corresponds to one of general relativity 
which is spatially closed, and has zero cosmological constant and zero pressure. 
Of the pressure-free models with zero cosmological constant, it seems that this 
one is the most favourable to the formation of condensation (2). We can see 
from the Newtonian analogue why this should be so. It is easily shown (9) 
that the total energy of any one of the fundamental particles of the Newtonian 
model has the sign of e, so if € is negative all these particles have velocities less 
than that required to escape from the matter within the spheres on whose 
boundaries they lie. ‘This makes the model more favourable to the growth of 
condensations. 
lhe second of equations (5.1) gives on integration 





87Ga . , 
= 7 sin? $6, (5.2) 
. 47Gau ; 
t+t,= “ae (? — $in @), ‘ (5-3) 
3 


where fy is a Constant of integration which we may take to be zero by choosing 
the origin of 6 at t=o. ‘The observations which the model is required to fit 
are the present values of Hubble’s constant, R/R, and the average density. ‘Taking 
these as 


p=2x 10% g/cm’, R/R=2:8x 10” (yrs.) 4, 
; ; - ‘ G 
we find (2) K - =9°3 x 10° yrs., 
3 


and the value of @ referring to the present time is 
6=1-15 rad. 
We may now return to (3°17) which governs the growth of a condensation. 
If we transform the independent variable from ¢ to @ by (5.3) and express py as 
a function of 6 we find, using (5.1) and (5.2), that (3.17) becomes 


d*h dh 3 k? dg 
—. + cot 46— —h| ~cosec” $@ — = = | =o. (5.4) 
de 2" do 2 2" ~ BE dp a4 
lo proceed further we must make some assumptions about the equation 
of state. For the reasons explained at the end of Section 3, let us suppose that 
p = g(p)=«p"*; 

then proceeding as in Section 3 we find that (5.4) becomes 
d*h 
dé? 


3 27K 
— 2/3 


7A a3) 3 x | =O, ( 


2 Gp; 5) 


wn 


se dh “osec? 16 
+ cot 207 h cosec* } 


where a perturbation of wave-length A, originated the disturbance when the 
density was p,. Suppose that this was at time ¢=¢, when 6=6,. As in 
Section 3, it follows that the disturbance can continue to grow if (3.24) is satisfied. 

Let us get an idea of the behaviour of the solution of (5.5) by using the approxt- 


mations 


cos @~1, sin 6~H 
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‘The equation then becomes 


d*h dh 
ee ae i ae 5.6 
O° aR 2075 4n*h=o, (5.6) 
where ‘iis 3. 2m 
~ 2 Gp*,,*° 


The solution of (5.6) is 


h= Ae + Bo”, 
where A and B are arbitrary constants, and where 
m,, m= —$+3V 14+ 16n?. 
As expected, the disturbance will grow if n?>o. 

From this approximate solution it seems that the most favourable case for 
condensation is that when n® has its maximum value, i.e. 3/2. ‘This is obtained 
by taking p=o, and supposing the model filled with dust. Condensations will 
then form for perturbations of any wave-length. If m*= 3/2, (5.5) can be solved 
exactly in terms of elementary functions. A particular solution is 

h= cosec? 36 cot 30, 
and the general one can be found by substituting in (5.5) 
h=v(8) cosec? $4 cot $8. 
This eventually leads to the solution 
h= — $A6 cosec? $6 cot 36+ A(2 +3 cot? }6)+ B cosec? 36 cot $8, (5-7) 
where A and B are arbitrary constants. 

We now wish to find out how quickly a small perturbation can grow. Sup- 

pose that the time ¢, is early enough for us to put 
cos #, ~ 1, sin 6, ~ 6). 

Then (5.7) gives 

h, ~ 8B0, 3+ ,A6,?. (5.8) 
If we assume that /, is initially small, we must have B very small since @, is small ; 
however, A need not be as small as f,.* At the present time, say when 6~1, 
we have from (5.7) 

h~ A; 

therefore from (5.8) h ~ 10h,0,°?. (5.9) 

It is now easy to see that the condensation process does not take place quickly 
enough to account for the formation of the nebulae. For a collection of V 
molecules of an ideal gas the fluctuation in the density is 

p 
and for a nebula of hydrogen NV ~ 3 x 10%. Hence we may take 


’ 


h,~ 10°, 
Supposing the perturbation to have occurred, say, 1000 years after the singular 
state, we have 
#,~ 10°, 
and from (5.9) h~ 10°, 
Thus although the disturbance has grown in magnitude, it is still quite negligible, 


* One might at first think that, since A and B are independent constants, one could satisfy (5.8) 
with /, small and yet have A and B large. However, it must be remembered that d//dt must 
also be small at time ¢,, and this leads to the conclusion given 
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and there seems no hope whatever of accounting for the formation of the nebulae 
by this mechanism if one assumes fluctuations of the values to be expected on 
ordinary statistical theory. ‘The difficulty cannot be overcome by taking a 
sufficiently small value for 6, because at such times the gravitational effect of the 
pressure could not be ignored and the model (5.2)-(5.3) would not apply. (For 
a fuller discussion, see (2)). 

The foregoing analysis applies to the “ closed”’ model with A=o. The 
‘‘ open ”’ models with A = 0, obtained by putting € >o in (5.1), give similar results ; 
that is to say, small perturbations cannot grow into nebulae in the time available.* 

6. Models with A+o.—The foregoing analysis is not easy to apply to models 
with non-zero cosmological constant, because, in general, elliptic functions are 
needed to integrate the second of equations (5.1). I shall not deal with this 
problem in detail here. I have previously discussed the formation of condensa- 
tions in Lemaitre’s model (2), and shown that this model seems hardly more 
favourable to the condensation process than those with A=o. Considered in 
the light of equation (3.17), the main difference between Lemaitre’s model and 
those with A =o is that in the former R/R is nearly zero for a certain period, and 
at this time the model resembles the Einstein universe. If this period were long 
enough it would certainly allow the formation of the nebulae provided dg/dp were 
small. However, if one supposes, as Lemaitre did, that the model started about 
5 x 10° years ago, this state of affairs cannot go on for very long. It might be 
worth while to investigate this problem further, but I have little doubt that 
detailed calculations would show that the condensations form too slowly to 
account for the formation of the nebulae. 

It is clear that if one allows A to be non-zero, models exist in which the nebulae 
could have formed. If the pressure is small enough to be neglected, and if R/R 
is small for a long period, (3.17) shows that / will increase very rapidly provided 
the density ic not too small. A model with these characteristics is the one 
investigated by Eddington and Lemaitre which expands from the Einstein state, 
though this is not the only one. 

7. Condensations in the steady-state theory.—The only plausible cosmology 
which has a theory of the formation of nebulae is the steady-state theory of Bondi, 
Gold and Hoyle. Sciama (11) has worked out the mechanism in some detail, 
and in his work Jeans’ formula plays a fundamental part. It is therefore important 
to see whether the use of Jeans’ formula in the steady-state theory can be justified. 

Here, however, one meets a difficulty. ‘The steady-state theory, though an 
interesting and imaginative cosmology, has not yet developed either a dynamics 
or a gravitation theory which carries any conviction. There is, though, a 
Newtonian model of the theory, developed with much ingenuity by McCrea (7). 
I shall try to adapt McCrea’s work to the problem in hand. 

According to McCrea, the Newtonian world models will accommodate the 
steady-state universe if one adds two postulates: (a) that mass and energy are 
convertible; (6) that the density of gravitational mass, instead of being p, is 

o=p+ pc, (7-1) 





* Attempts have been made to overcome the difficulty by postulating turbulence in the primordial 
medium; this gives rise to large fluctuations from which the nebulae could have formed. However, 
this proposal can hardly be regarded as a satisfactory solution of the problem unless it can be 
shown first how the turbulence might have arisen, and secondly that it would have survived the 
period of high pressure. Neither of these questions has yet been answered. 


8 
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c being the velocity of light. ‘The equations of the unperturbed model, after 
allowing for McCrea’s two extra assumptions, are 
3h +3f?= — 47G(po + 3Poe), (7.2) 
Po + 3pPof + 3Pofe? =o. (7-3) 
These differ from the equation (3.5) of the standard theory in the presence of 
the terms involving py. The solution of (7.2) and (7.3) which gives the Newtonian 
analogue of the steady-state theory is 
f=RI/R, R= Ae, (7.4) 
Po = 3¢°/82Ga?, Po= — 3¢4/37Ga’, (7-5) 
where A and a are arbitrary constants. 
‘To extend this work to the perturbed Newtonian model considered in Section 
2, let us make two further assumptions : 
(c) equation (7.1) holds even in the non-uniform model ; 
(d) the rate of creation of matter is independent of the density, and is 
gc?/87Ga*_g/cm® sec. 
The latter is the rate in the model (7.4) and (7.5). Although the theory is not 
quite explicit about the dependence of the creation rate on the local density, 
postulate (d) seems a reasonable interpretation of the authors’ intentions (1). 
The last term on the left in (7.3) represents the creation of matter, which, of 
course, causes an alteration in the equation of continuity. 
With these assumptions the equations for the model subject to spherically 
symmetric perturbations are 


eu ou 1 0p 
te OB (7.6) 
0 } ou 
= + ue tp, + 320 f= o, (7-7) 
nt ee ai sf). (7.8) 


These correspond to (3.1), sais and (3.4), and in the unperturbed model are 
equivalent to (7.2)-(7.5).* 

If, following McCrea, we regard the negative pressure p, of the homogeneous 
model as a “ zero-point” pressure, we may suppose that the equation of state 
of the gas of the model takes the form 


P=po+se), 
where g(py)=0. Using the notation of (3.3), we find that for small perturbations 
dg 
ides” * , (7.9) 
op dq 
5, =u is ; (7.10) 


where, as usual, dg/dp is to be calculated at p= pp. 


* A difficulty should be pointed out here. In (7.6) F is the force per unit inertial mass, whereas 
in (7.8) Fis the force on unit “ passive”’ gravitational mass (i.e. mass on which gravitation acts). 
One cannot be sure that these F’s are really the same, especially as we are taking the density of 

‘active ’’ (gravitation-producing) mass as p+ 3pc-*, whereas the density of inertial mass is p. 
My justification for equating the two F’s lies in McCrea’s demonstration of the self-consistency 
of his use of gravitational mass, which is similar to that here. 
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We may now treat (7.6)-(7.8) as we treated the corresponding equations in 
Section 3, using (7.2), (7. 3), (7-9) and (7.10) where necessary. ‘The result, 
corresponding to (3.13) is 

& + 2rfas’ + r7f*a" + wf (8 + 3poc py) 


1£2 f 3Po . 6f OPo 127G 0 
+ rw’ f* (o+4 + Pe) + aft (18+ 7 i cp, 5 a) 


dgf{ , 2’ w d, O7Gpy 
~ plot F)- SE (ors of SP). 
We can simplify this by using (7.4) and (7.5), which reduces it to 


e c? 5c c gc? 1 dg 
@+2r-w'+r?— ow” + i +6r—se’ + fl —wlI- s 
a a* a a* 2a* dp 








dg ~ 2w’ 
=> + Fs 7.11 
dp (« r ) (7 ) 
Now change to a new variable x by the transformation 
r=xR(t); 
this changes (7.11) to 
ee, tf, Bake, 2). 
ot? aot 2a cdp}] dp Ox® x Ox 
If we assume a solution of the form 
sin kx 
= h(t), 
an” (t) 
we find the following — for * : 
I dg k? dg 
ht I-s=> =0. 7.12 
h+2 +2, a = R: dp (7-12) 


Evidently we may a the term mm cMeldp in the bracket, which in the case of a 
gas at ordinary temperatures will be small compared with unity, and it is then 
clear that # cannot increase with time. In the most favourable case, in 
which dg/dp =o, the solution of the equation (7.12) then is 
h= Ae™ + Be, 

—_ my, m= ~<[-5+V7]. 

a 
Since the condensation 8p/p, is proportional to w (because p, is constant) it follows 
that in the steady-state theory condensations cannot form as a result of small perturba- 
tions in the average density. 

In Sciama’s work the density p of the condensing gas-cloud is taken to be 
equal to the average density py, and from the foregoing analysis it seems that no 
condensation would in fact take place if p and p, were equal, or differed by a 
small quantity. Strictly, therefore, Sciama’s theory of the formation of the 
nebulae is not valid. However, in the steady-state theory , large perturbations—of 
the order of several times the average density—are possible as a result of 
clouds of gas breaking away from the accretions of already-existing nebulae. 
It may be that with such large perturbations new nebulae can form in the time 
allowed; but this is not obvious, and until it has been proved one cannot say 
with any confidence that the steady-state theory explains the formation of the 
nebulae. 
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8. Conclusion.—In this paper I have been studying the applicability of Jeans’ 
formula to the expanding universe. Only small perturbations of the average 
cosmic density py have been investigated: this is the only case to which Jeans’ 
work strictly applies because in no other circumstances is one allowed to treat 
the unperturbed gas as uniform. 

The result of this work is that, with a suitable interpretation of the quantities 
involved, Jeans’ formula applies in the present phase of the expansion (except 
in the steady-state theory). Thus disturbances of sufficiently large wave-lengths 
result in instability of the cosmic gas, which causes it to start condensing. If 
one supposes that independent perturbations take place in different regions of the 
universe, then matter will condense into separate conglomerations. The size 
and mass of the condensing regions are connected by (3.21), in which the epoch 
of the perturbation implicitly enters through pp, which, of course, depends on the 
epoch. 

Where Jeans’ original work is misleading (though through no fault of Jeans) 
is in the suggestion (from the solution of (2.12)) that the condensations begin 
to form at an exponential rate. It is only in the Einstein universe that this 
actually happens. When the expansion is taken into account one finds that the 
condensation process takes place very slowly—too slowly to account for the 
formation of the nebulae unless one supposes either that the perturbations are 
much larger than those to be expected on ordinary statistical theory, or that the 
universe has a longer time-scale than that provided by the point-source models. 

Since so much of the matter in the universe is now condensed into nebulae, 
clouds of gas with density very nearly equal to the average density p, must be 
rare at the present time. ‘Thus Jeans’ formula as originally derived and as 
studied in this paper, though it may prove to be important in explaining the 
growth of the nebulae, has little practical application to contemporary astro- 
physics. 

Naturally one would like to know whether Jeans’ type of instability operates 
for clouds of gas with densities larger than py. This represents an unsolved 
problem, even in the spherically-symmetric case, but the stability will presumably 
depend markedly on the equation of state; this follows from the fact, mentioned 
in the introduction, that an isothermal gas sphere of sufficient size is unstable 
whereas adiabatic spheres with y > 4/3 are stable so far as is known. I hope to 
deal with this question in another paper. 
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Note added by D. W. Sciama on“ Jeans’ Formula for Gravitational Instability’. 
—Bonnor’s paper demonstrates in a very clear way what has been generally 
believed for some time, namely that infinitesimal density fluctuations will not 
give rise to galaxies in the times available. ‘This difficulty led Gamow to suggest 
that large density fluctuations would be produced by the high Mach number 
(~10) compressible turbulence that he invoked. Unfortunately the analogue 
of the Jeans’ criterion is not known for this type of motion. 

At first sight the situation is even worse for the steady-state theory. For the 
material density must increase by a factor three in order that self-gravitation 
should exactly compensate the cosmical repulsive forces that the steady-state 
theory must invoke. However, the situation is saved by the fact that these large 
density fluctuations do not have to form in an initially homogeneous universe. 
Indeed, all we have to show is that galaxies are self-propagating. We can hence 
use the fact that a galaxy will produce a large density fluctuation in its wake in 
a systematic way, so that the region of large density will not be highly turbulent, 
and Jeans’ criterion can be applied to it. 

In my paper I gave some sketchy calculations of the process of galaxy 
formation in a steady-state universe. Further details are given in my thesis 
(On the Origin of Inertia, Chapter V, Cambridge University 1953). In my paper 
I glossed over the cosmological aspects of the problem in order to concentrate 
on the astrophysical ones. ‘That is why I assumed that the system was Newtonian 
in the region of interest. In fact one must consider only that part of the wake 
in which the density exceeds 3p. If we go further and require the density to 
exceed 4p, then the net gravitating density to be used in the Jeans formula is 
of order p (because of the cosmical repulsion). | From the formula for p,, given 
on page 5 of my paper we see that p,, >4p when 


: 1 (3M,\'3 
r, sin 6,,: (32) 


where p, is the density at the point p, r, is the distance of p from the parent 
galaxy P, @,, is the angle Pp makes with the peculiar velocity vector of P, and 
M, is the mass of P. Gravitational instability only occurs in the part of the 
wake which satisfies this equation. ‘This restriction does not alter the order 
of magnitude estimates made in my paper. 

I should add that various other points are considered in the thesis, e.g. the 
effect of the Hubble velocities initially present in the inter-galactic gas, and the 
stress produced by newly-created matter which is presumably under the 
gravitational influence of the parent galaxy only after it appears. Again these 
effects do not alter the order of magnitude estimates. I conclude then that 
Bonnor’s calculations do not affect my theory of galaxy formation. 
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8. Conclusion.—In this paper I have been studying the applicability of Jeans’ 
formula to the expanding universe. Only small perturbations of the average 
cosmic density py have been investigated: this is the only case to which Jeans’ 
work strictly applies because in no other circumstances is one allowed to treat 
the unperturbed gas as uniform. 

The result of this work is that, with a suitable interpretation of the quantities 
involved, Jeans’ formula applies in the present phase of the expansion (except 
in the steady-state theory). ‘Thus disturbances of sufficiently large wave-lengths 
result in instability of the cosmic gas, which causes it to start condensing. If 
one supposes that independent perturbations take place in different regions of the 
universe, then matter will condense into separate conglomerations. The size 
and mass of the condensing regions are connected by (3.21), in which the epoch 
of the perturbation implicitly enters through pp, which, of course, depends on the 
epoch. 

Where Jeans’ original work is misleading (though through no fault of Jeans) 
is in the suggestion (from the solution of (2.12)) that the condensations begin 
to form at an exponential rate. It is only in the Einstein universe that this 
actually happens. When the expansion is taken into account one finds that the 
condensation process takes place very slowly—too slowly to account for the 
formation of the nebulae unless one supposes either that the perturbations are 
much larger than those to be expected on ordinary statistical theory, or that the 
universe has a longer time-scale than that provided by the point-source models. 

Since so much of the matter in the universe is now condensed into nebulae, 
clouds of gas with density very nearly equal to the average density p, must be 
rare at the present time. ‘Thus Jeans’ formula as originally derived and as 
studied in this paper, though it may prove to be important in explaining the 
growth of the nebulae, has little practical application to contemporary astro- 
physics. 

Naturally one would like to know whether Jeans’ type of instability operates 
for clouds of gas with densities larger than py). ‘This represents an unsolved 
problem, even in the spherically-symmetric case, but the stability will presumably 
depend markedly on the equation of state; this follows from the fact, mentioned 
in the introduction, that an isothermal gas sphere of sufficient size is unstable 
whereas adiabatic spheres with y > 4/3 are stable so far as is known. I hope to 
deal with this question in another paper. 


The University, 
Liverpool, 3 : 
1956 December. 
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Note added by D. W. Sciama on“ Jeans’ Formula for Gravitational Instability’. 
—Bonnor’s paper demonstrates in a very clear way what has been generally 
believed for some time, namely that infinitesimal density fluctuations will not 
give rise to galaxies in the times available. This difficulty led Gamow to suggest 
that large density fluctuations would be produced by the high Mach number 
(~10) compressible turbulence that he invoked. Unfortunately the analogue 
of the Jeans’ criterion is not known for this type of motion. 

At first sight the situation is even worse for the steady-state theory. For the 
material density must increase by a factor three in order that self-gravitation 
should exactly compensate the cosmical repulsive forces that the steady-state 
theory must invoke. However, the situation is saved by the fact that these large 
density fluctuations do not have to form in an initially homogeneous universe. 
Indeed, all we have to show is that galaxies are self-propagating. We can hence 
use the fact that a galaxy will produce a large density fluctuation in its wake in 
a systematic way, so that the region of large density will not be highly turbulent, 
and Jeans’ criterion can be applied to it. 

In my paper I gave some sketchy calculations of the process of galaxy 
formation in a steady-state universe. Further details are given in my thesis 
(On the Origin of Inertia, Chapter V, Cambridge University 1953). In my paper 
I glossed over the cosmological aspects of the problem in order to concentrate 
on the astrophysical ones. That is why I assumed that the system was Newtonian 
in the region of interest. In fact one must consider only that part of the wake 
in which the density exceeds 3p. If we go further and require the density to 
exceed 4p, then the net gravitating density to be used in the Jeans formula is 
of order p (because of the cosmical repulsion). From the formula for p,, given 
on page 5 of my paper we see that p,, >4p when 


«nr 
r, sin 6, ($5) 


where p, is the density at the point p, 7, is the distance of p from the parent 
galaxy P, @, is the angle Pp makes with the peculiar velocity vector of P, and 
M, is the mass of P. Gravitational instability only occurs in the part of the 
wake which satisfies this equation. ‘This restriction does not alter the order 
of magnitude estimates made in my paper. 

I should add that various other points are considered in the thesis, e.g. the 
effect of the Hubble velocities initially present in the inter-galactic gas, and the 
stress produced by newly-created matter which is presumably under the 
gravitational influence of the parent galaxy only after it appears. Again these 
effects do not alter the order of magnitude estimates. I conclude then that 
Bonnor’s calculations do not affect my theory of galaxy formation. 
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